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2.  Technical  Objectives: 

Wide-bandwidth  microwave  modulation  of  semiconductor  strained  quantum-well  lasers 
using  distributed-feedback  structures  plays  an  important  role  in  high  bit-rate  optical 
communication  systems.  It  is  therefore  critical  to  conduct  research  on  the  device  physics  of 
strained  semiconductor  quantum-well  lasers  under  high-frequency  modulation  and  on  the 
creation  of  new  techniques  to  characterize  high-speed  quantum-well  lasers.  The  goals  of 
this  project  are  (1)  to  design  and  fabricate  high  bandwidth  distributed-feedback  quantum- 
well  lasers  using  strain  effects  and  (2)  to  combine  microwave  and  optical  measurement 
techniques  to  investigate  the  physics  and  device  performance  of  quantum-well  lasers  under 
high-speed  modulation  conditions.  The  project  is  interdisciplinary  in  nature  since 
optoelectronic  device  technology,  microwave  and  optical  measurement  techniques  will  be 
introduced  to  the  study  of  high-frequency  modulation  of  quantum-well  lasers  for  optical 
communication  systems. 

3.  Approach 

a.  Use  electrical  current  injection  to  measure  the  transport  and  intrinsic  high-speed 
parameters  of  semiconductor  laser  and  study  the  temperature  sensitivity  of  these 
parameters. 

b.  Use  optical  injection  by  pumping  in  the  absorption  region  (cross-absorption 
modulation)  to  study  the  intrinsic  small-signal  modulation  response  of 
semiconductor  lasers. 

c.  Use  optical  injection  by  pumping  in  the  gain  region  (cross-gain  modulation)  to  study 
high-speed  gain  competition  phenomena  of  semiconductor  lasers  and  realize 
wavelength  conversion. 

d.  Demonstrate  a  novel  wavelength  conversion  experiment  by  cross-absorption 
modulation  in  an  integrated  electro-absoiption  modulator/laser. 

e.  Demonstrate  a  novel  wavelength  conversion  experiment  by  cross-gain  modulation  in 
a  multisection  integrated  laser/modulator  device. 

f.  Conduct  experiments  on  four-wave  mixing  in  conventional  and  novel  three-electrode 

DFB  lasers  and  study  their  polarization  dependence. 

4.  Accomplishments 

A.  Broaden  Bandwidth  by  Optical  Pumping  in  the  Absorption  Region. 


1 .  Contract  Title: 

Number: 

PI: 


1 


We  study  both  theoretically  and  experimentally  high-speed  optical  injection  in  the 
absorption  region  of  a  quantum-well  laser,  and  compare  the  results  with  those  from 
electrical  injection  including  the  carrier  transport  effect.  We  show  that  the  main  difference 
between  the  two  responses  is  the  low  frequency  roll-off.  By  using  both  injection  methods, 
we  obtain  more  accurate  and  consistent  measurements  of  many  important  dynamic  laser 
parameters  which  include  the  differential  gain,  carrier  lifetime,  K  factor,  and  gain 
compression  factor.  Temperature-dependent  data  of  the  test  laser  are  presented  which 
show  that  the  most  dominant  effect  is  the  linear  degradation  of  differential  gain  and 
injection  efficiency  with  increasing  temperature.  While  the  K-factor  is  insensitive  to 
temperature  variation  for  multiple-quantum-well  lasers,  we  find  that  the  carrier  capture  time 
and  nonlinear  gain  suppression  coefficient  decreases  as  temperature  increases. 

B.  Small-Signal  Cross-Gain  Modulation  in  Semiconductor  Lasers 

We  derive  new  rate  equations  for  the  high-speed  frequency  response  of  a  test  quantum-well 
laser  in  the  presence  of  an  external  microwave  modulated  optical  pump  in  the  gain  region. 
The  model  includes  the  effects  of  pump-induced  stimulated  recombination  and  cross  gain 
saturation.  Expressions  for  the  small-signal  modulation  response  of  the  test  laser  under 
gain  modulation  are  derived.  We  use  a  multiple-quantum-well  InGaAlAs  Fabry-Perot  (FP) 

laser  at  1.552  (im  as  the  test  laser  and  a  DFB  laser  at  1 .542(j.m  as  the  external  pump.  Our 
data  shows  a  reduction  of  carrier  lifetime  with  increasing  optical  pumping,  a  shift  of  the 
probe  laser  threshold  current,  a  change  in  the  K  factor,  and  a  variation  of  the  relaxation 
frequency  with  different  pump  powers.  The  experimental  results  agree  very  well  with  the 
theoretical  results. 

C.  Wavelength  Conversion  by  Cross-Absorption  Modulation  Using  an 
Integrated  Electroabsorption  Laser/Modulator 

Wavelength  conversion  by  cross-absorption  modulation  has  been  demonstrated  using  an 
integrated  electroabsorption  modulator/distributed-feedback  laser  which  acts  as  the 
conversion  medium  and  converted  wavelength  light  source.  Dependencies  on  pump 
wavelength  and  modulator  voltage  have  been  investigated. 

D.  Dual  Pump  Four-wave  Mixing  in  a  Double-moded  Distributed-Feedback 
Laser 

We  present  a  comprehensive  theoretical  model  with  experimental  results  on  four-wave 
mixing  in  a  DFB  laser  with  two  pump  modes.  The  twin  lasing  modes  serve  as  pump 
waves  for  the  four-wave  mixing  which  takes  place  within  the  laser  cavity.  We  show  new 
experimental  measurements  of  the  two-pump  four-wave  mixing  conversion  efficiency 
spectrum  and  show  that  to  first  order  there  is  no  fundamental  dependence  of  the  conversion 
efficiency  on  pump  separation.  The  indicates  significant  improvement  in  the  probe  and 
conjugate  wavelength  separation  and  conversion  efficiency.  Our  theoretical  model  agrees 
quite  well  with  the  measured  data  and  predicts  the  key  experimentally  observed  phenomena 
including  the  relative  invariance  of  the  conversion  efficiency  on  pump  separation  and  the 
strong  enhancement  of  certain  conjugate  waves  when  the  probe  detuning  is  close  to  the 
pump  separation. 


4.  Significance  of  this  work: 

The  fundamental  research  conducted  here,  which  includes  electrical,  optical  absorption  and 
optical  gain  high-speed  responses,  allows  for  a  better  understanding  of  current  device 
limitations  on  broad-bandwidth  semiconductor  laser  performance,  and  will  enhance  high¬ 
speed  device  design.  Our  investigation  of  applications  for  strained  quantum-well 
optoelectronics  such  as  wavelength  conversion  by  cross-gain  and  cross-absorption 


2 


modulation  enable  wavelength  switching  in  dense  WDM  networks.  The  high-speed 
modulation  in  distributed-feedback  (DFB)  lasers  is  essential  for  enabling  the  growth  of 
future  optical  communication  systems.  Study  of  high-speed  modulation  in  integrated 
laser/modulator  devices  is  critical  to  develop  broad-bandwidth  transmitters. 


6.  Papers.  Patents,  and  Awards; 

(Publications) 

(i)  Journal  Papers 

1.  W.  Fang,  M.  Hattendorf,  S.  L.  Chuang,  J.  Minch,  C.  S.  Chang,  C.  G.  Bethea, 
and  Y.  K.  Chen,  "Analysis  of  temperature  sensitivity  in  semiconductor  lasers 
using  gain  and  spontaneous  emission  measurements,"  Appl.  Phys.  Lett.,  vol.  70, 
pp.  796-798,  1997. 

2.  J.  Minch,  C.  S.  Chang,  and  S.  L.  Chuang,  "Four-wave  mixing  in  a  distributed- 
feedback  laser,"  Appl.  Phys.  Lett.,  vol.  70,  pp.  1360-1362,  1997. 

3.  J.  Minch,  S.  L.  Chuang,  C.  S.  Chang,  W.  Fang,  Y.  K.  Chen,  and  T.  Tanbun-Ek, 
"Theory  and  experiment  on  the  amplified  spontaneous  emission  from  distributed- 
feedback  lasers,"  IEEE  J.  Quantum  Electron.,  vol.  33,  pp.  815-823,  1997. 

4.  J.  Minch,  C.  S.  Chang,  and  S.  L.  Chuang,  "Wavelength  conversion  in  distributed- 
feedback  lasers,"  TREE  J.  Selected  Topics  Quantum  Electronics,  vol.  3,  pp.  569- 
576,  1997. 

5.  W.  Fang,  A.  Hsu,  S.  L.  Chuang,  T.  Tanbun-Ek,  and  A.  M.  Sergent, 
"Measurement  and  modeling  of  distributed-feedback  lasers  with  spatial  hole 
burning,"  IEEE  J.  Selected  Topics  Quantum  Electronics,  vol.  3,  pp.  547-554, 
1997. 

6.  T.  Keating,  J.  Minch,  C.  S.  Chang,  P.  Enders,  W.  Fang,  S.  L.  Chuang,  T. 
Tanbun-Ek,  T.  K.  Chen,  and  M.  Sergent,  "Optical  gain  and  refractive  index  of  a 
laser  amplifier  in  the  presence  of  pump  light  for  cross-gain  and  cross-phase 
modulation,"  IEEE  Photon.  Technol.  Lett.,  vol.  9,  pp.  1358-1360, 1997. 

7.  J.  Minch,  S.  H.  Park,  T.  Keating  and  S.  L.  Chuang,  "Theory  and  Experiment  of 
InGaAsP  and  InGaAlAs  long-wavelength  strained  quantum-well  lasers,"  IEEE  J. 
Quant.  Electron.,  vol.  35,  pp.  771-782,  1999. 

8.  T.  Keating,  S.  H.  Park,  J.  Minch,  X.  Jin,  and  S.  L.  Chuang,  "Optical  gain 
measurements  based  on  fundamental  properties  and  comparison  with  many-body 
theory,"  J.  Appl.  Phys.,  vol.  86,  pp.  2945-2952,  1999. 

9.  A.  Hsu,  S.  L.  Chuang,  W.  Fang,  L.  Adams,  G.  Nykolak,  and  T.  Tanbun-Ek,  "A 
Wavelength-Tunable  Curved  Waveguide  DFB  Laser  with  Integrated  Modulator," 
IEEE  J.  Quant.  Electron.,  vol.  35,  pp.  961-969,  1999. 

10.  T.  Keating,  X.  Jin,  S.  L.  Chuang,  and  K.  Hess,  "Temperature  dependence  of 
electrical  and  optical  modulation  responses  of  quantum-well  lasers,"  IEEE  J. 
Quantum  Electron.,  vol.  35,  pp.  1526-1534,  1999. 

11.  J.  Minch  and  S.  L.  Chuang,  "Dual-pump  four-wave  mixing  in  a  double-mode 
distributed  feedback  laser,"  J.  Opt.  Soc.  Am.  B,  vol.  17,  pp.  53-62,  2000. 

(ii)  Presentations 

12.  J.  Minch,  C.  S.  Chang,  W.  Fang,  S.  L.  Chuang,  Y.  K.  Chen,  and  T.  Tanbun-Ek, 
"Theory  and  experiment  on  the  amplified  spontaneous  emission  from  DFB  lasers," 
Conference  on  Lasers  and  Electro-Optics,  Anaheim,  CA,  June  2-7, 1996. 

13.  W.  Fang,  C.  G.  Bethea,  Y.  K.  Chen,  T.  Tanbun-Ek,  and  S.  L.  Chuang, 
"Measurements  of  spatially  dependent  carrier  lifetimes  in  1.55  pm  DFB  quantum- 


3 


well  lasers  using  microwave  modulation,"  Conference  on  Lasers  and  Electro- 
Optics,  Anaheim,  CA,  June  2-7,  1996. 

14.  *S.  L.  Chuang,  C.  S.  Chang,  J.  Minch,  and  W.  Fang,  "Amplified  spontaneous 
emission  spectroscopy  of  strained  quantum- well  lasers-Theory  and  experiment," 
Optoelectronics’97  Symposium,  San  Jose,  Feb.  1997.  (Invited). 

15.  W.  Fang,  S.  L.  Chuang,  T.  Tanbun-Ek,  and  Y.  K.  Chen,  "Modeling  and 
experiment  of  1.55|im  integrated  electroabsorption  modulator  with  distributed- 
feedback  lasers,"  SPIE,  vol.  3006,  Edited  by  Y.  S.  Park  and  R.  V.  Ramaswamy, 
Photonics  West,  Optoelectronics'97  Symposium,  San  Jose,  Feb.  1997. 

16.  J.  Minch,  C.  S.  Chang,  and  S.  L.  Chuang,  "Wavelength  conversion  using 
distributed-feedback  lasers,"  Conference  on  Lasers  and  Electro-Optics,  Baltimore, 
MD,  May,  1997. 

17.  J.  Minch,  C.  S.  Chang,  and  S.  L.  Chuang,  "Wavelength  conversion  using  two- 
pump  four-wave  mixing  in  a  double-moded  distributed-feedback  laser,"  Laser  and 
Electro-optics  Society  Annual  Meeting,  MN3,  San  Francisco,  CA,  Nov.  10-13, 
1997. 

18.  T.  Keating,  J.  Minch,  C.  S.  Chang,  and  S.  L.  Chuang,  "Optimal  refractive  index 
changes  for  cross-gain  modulation  and  cross-phase  modulation,"  Photonics  West, 
Physics  and  Simulation  of  Optoelectronic  Devices  VI,  San  Jose,  CA,  1998. 

19.  A.  Hsu,  W.  Fang,  S.  L.  Chuang,  T.  Tanbun-Ek,  C.  Bethea,  and  R. 
People"Integrated  tunable  laser  with  mode  selection  modulator,"  Photonics  West, 
Optoelectronic  Integrated  Circuits  n,  SPIE,  vol.  3290,  Edited  by  S.  Y.  Wang  and 
Y.  S.  Park,  San  Jose,  CA,  1998. 

20.  A.  Hsu,  W.  Fang,  and  S.  L.  Chuang,  "Modeling  of  normal  and  backward 
integrated  electroabsorption  modulator  and  laser,"  Integrated  Photonic 
Research., Victoria,  Canada,  April  1998. 

21.  *S.  L.  Chuang,  S.  H.  Park,  J.  Minch,  and  T.  Keating,  "Strained  zinc-blende  and 
wurtzite  quantum- well  lasers:  theory  and  comparison  with  experiment," 
Semiconductor  Science  and  Technology'98,  La  Jolla,  CA,  Sept.  1998. 

22.  E.  Young,  J.  Minch,  and  S.  L.  Chuang,  "Polarization  resolved  four-wave  mixing 
in  a  distributed-feedback  laser,"  Photonics  West,  Physics  and  Simulation  of 
Optoelectronic  Devices  VII,  San  Jose,  CA,  Jan.  1999. 

23.  A.  Hsu  and  S.  L.  Chuang,  "Wavelength  conversion  by  cross-absorption 
modulation  using  an  integrated  electroabsorption  modulator  /  laser,"  Conference  on 
Lasers  and  Electro-Optics,  Baltimore,  MD,  May  1999. 

24.  A.  Hsu,  S.  L.  Chuang,  and  T.  Tanbun-Ek,  "Multi-wavelength  conversion  in  an 
integrated  wavelength-tunable  laser-modulator,"  Photonics  West,  OEIC  IV,  Proc. 
SPIE,  Ed.,  Y.  S.  Park  and  R.  T.  Chen,  vol.  3950,  2000  (in  press). 

*  invited  talk 


4 


ATTACHMENTS 


Five  major  representative  publications 


IEEE  JOURNAL  OF  QUANTUM  ELECTRONICS,  VOL.  35,  NO.  5,  MAY  1999 


771 


Theory  and  Experiment  of  Ini-^Ga^AsyPi-y 
and  Ini-^-yGa^AlyAs  Long-Wavelength 
Strained  Quantum- Well  Lasers 

J.  Minch,  S.  H.  Park,  T.  Keating,  and  S.  L.  Chuang,  Fellow,  IEEE 


Abstract — We  present  a  comprehensive  model  for  the  calcu¬ 
lation  of  the  bandedge  profile  of  both  the  Ini-xGaxAsvPi-y 
and  Ini_x_vGarAlyAs  quantum-well  systems  with  an  arbitrary 
composition.  Using  a  many-body  optical  gain  model,  we  compare 
the  measured  net  modal  gain  for  both  material  systems  with 
calculations  from  the  realistic  band  structure  including  valence 
band  mixing  effects.  Calibrated  measurements  of  the  side  light 
spontaneous  emission  spectrum  based  on  its  fundamental  relation 
to  the  optical  gain  spectrum  give  values  for  the  radiative  current 
density.  These  measurements  allow  us  to  extract  the  relationship 
between  total  current  density  and  carrier  density.  A  fit  of  this 
relation  yields  values  for  the  Auger  coefficient  for  each  material 
system. 

Index  Terms — Auger  recombination,  gain,  InGaAsP,  InGaAlAs, 
semiconductor  lasers,  strained  quantum-well  lasers. 


L  Introduction 

THERE  are  two  main  material  systems  used  to  fab¬ 
ricate  long-wavelength  semiconductor  lasers.  While 
the  Ini-xGaxAsyPi_y  system  has  been  firmly  established 
as  a  candidate  for  high-performance  lasers  [1],  [2],  the 
Ini-x-yGasAlyAs  system  has  only  recently  begun  to  yield 
quality  results  [3]-[5].  The  interest  in  the  Ini_x_yGaxAlyAs 
system  stems  from  its  larger  conduction  band  offset  of 
AEc/AEg  =  0.7  compared  to  a  value  of  AEc/AEg  =  0.4 
for  Ini_xGaxAsyPi_y.  This  larger  conduction  band  offset  had 
been  predicted  to  result  in  better  electron  confinement  in  the 
conduction  band  and,  therefore,  a  higher  temperature  stability. 

Previous  reports  have  been  given  on  the  model  for  the 
band  structure  of  these  two  systems  [6],  theoretical  calcu¬ 
lations  comparing  the  systems  in  terms  of  threshold  current 
density  and  differential  gain  [7],  as  well  as  the  tempera¬ 
ture  dependence  of  the  threshold  current  [8].  There  have 
also  been  measurements  of  the  contributions  of  the  various 
recombination  mechanisms  for  each  system  [9].  However, 
little  has  been  published  directly  comparing  the  gain  and 
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recombination  mechanisms  of  these  two  systems,  or  comparing 
the  existing  many-body  optical  gain  models  with  experimental 
measurements. 

In  this  paper,  we  present  a  comprehensive  model  for  the 
calculation  of  the  band  structure  of  samples  made  of  either 
of  these  systems  with  an  arbitrary  composition.  Then,  using 
a  many-body  model  for  the  optical  gain  and  spontaneous 
emission,  we  compare  the  model  to  experimental  data  for 
each  material  system.  Also,  from  simultaneous  measurements 
of  both  the  spontaneous  emission  spectrum  from  the  side  of 
the  lasers  and  the  gain  spectrum  from  the  facet  emission,  we 
are  able  to  compare  the  relative  contributions  of  the  different 
recombination  mechanisms  in  each  laser  system. 

In  Section  II,  we  present  our  model  for  the  calculation  of 
the  bandedge  of  each  material  system  and  the  quantum-well 
(QW)  band  edge  discontinuities,  followed  by  the  model  for 
the  optical  gain  and  spontaneous  emission.  In  Section  HI,  we 
describe  the  structure  of  the  lasers  as  well  as  our  experimental 
measurements.  Section  IV  describes  the  comparison  of  the 
optical  gain  data  with  calculations  using  our  model,  and  the 
extraction  of  the  relative  contributions  of  the  recombination 
mechanisms  in  each  system.  Finally,  a  brief  conclusion  is 
given  in  Section  V. 

n.  Theoretical  Model 

A .  Calculation  of  Bulk  Bandgap  with  Strain 

To  obtain  most  parameters  for  both  the  Ini_xGaxAsyPi_y 
and  Ini_x..yGaxAlyAs  material  systems,  a  linear  interpolation 
between  the  parameters  of  the  relevant  binary  semiconductors 
is  used.  The  interpolation  formulas  for  all  physical  parameters 
P  used  in  the  calculation  of  the  band  edge,  except  for  the 
bandgap,  are  given  as  [10],  [11] 

P  (Ini — a;Gax  Asy  Pi  -  y ) 

=  P(GaAs)rcy  -I-  P(GaP)x(l  —  y)  +  P(InAs)(l  —  x)y 
4-  P(InP)(l  —  x)(l  -  y)  (1) 

P(Ini-x_yGaxAly  As) 

=  P(InAs)(l  -  x  -  y)  4-  P(GaAs)x  4-  P(ALAs)?/.  (2) 

The  material  parameters  of  the  binary  semiconductors  can  be 
found  in  Table  I. 

The  one  exception  to  the  linear  interpolation  is  the  formula 
for  the  unstrained  bandgap.  For  each  material  system,  this 
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TABLE  I 

Parameters  for  the  Calculation  of  Strain  and  Bandgap  Energy  for  the  Ini_xGax  AsyPi-y  and  Ini_x-yGaz  AlyAs  Material  Systems  [10],  [11] 


Parameter 

Symbol(unit) 

GaAs 

InAs 

InP 

GaP 

AlAs 

Lattice  Constant 

a(A) 

5.6533 

6.0584 

5.8688 

5.4505 

5.660 

Elastic  Stiffness  Constant 

Cu(lOltdynfcrn2) 

11.879 

8.329 

10.11 

14.05 

12.5 

Elastic  Stiffness  Constant 

Ci2(10lidyn/cni2) 

5.376 

4.526 

5.61 

6.203 

5.34 

Hydrostatic  deformation  potential 

for  conduction  band 

«c(eV) 

-7.17 

-5.08 

-5.04 

-7.14 

-5.64 

for  valence  band 

av{eV) 

1.16 

1.00 

1.27 

1.70 

2.47 

Shear  deformation  potential 

for  valence  band 

b{eV) 

-1.7 

-1.8 

-1.7 

-1.8 

-1.5 

Valence  band  parameter 

7i 

6.8 

20.4 

4.95 

4.05 

3.45 

72 

1.9 

8.3 

1.65 

0.49 

0.68 

73 

2.73 

9.1 

2.35 

1.25 

1.29 

Electron  effective  mass 

me/m0 

0.067 

0.023 

0.077 

0.25 

0.15 

Heavy-hole  effective  mass 

0.50 

0.40 

0.60 

0.67 

0.79 

quantity  is  given  as 
Ini_xGaxAsyPi_y 

Eg(x,  y)  =  1.35  +  0.642a:  -  l.lOly  +  0.758a:2  +  O.lOlj/2 
-  0.159xy  -  0.28x2y  +  0.109xy2  eV  (3) 

Eg(x,  y)  =  0.36  +  2.093j/  +  0.629s;  4-  0.577 y2  4-  0.436s;2 

4*  1.013xy  —  2.0xy(l  -  x  -  y)  eV.  (4) 


where  ac  and  av  are  the  conduction-band  and  valence-band 
hydrostatic  deformation  potentials,  and  b  is  the  valence-band 
shear  deformation  potential. 

The  strained  bandgaps  can  then  be  expressed  as 

Ec—hhfay  y )  —  Eg 2/)  “b  6Ec(x,  y)  6EhH{x,  y)  (14) 
and 

Ec-ih(x ,  y)  =  Eg(x,  y)  +  8Ec(x ,  y)  -  6Eih(x,  y).  (15) 


The  effects  of  strain  are  calculated  in  the  following  way. 
First,  the  strain  in  the  plane  of  the  epitaxial  growth  is 

_  aG  —  a 
evv 


6  —  £xx  —  €yv  — 


(5) 


where  a  is  the  lattice  constant  of  the  quaternary  epitaxial  layer 
and  aQ  is  the  lattice  constant  of  the  substrate  which  is  assumed 
to  be  InP.  The  strain  in  the  perpendicular  direction  can  be 
expressed  as 


e“-  2Cn€ 


(6) 


where  Cu  and  C\2  are  the  elastic  stiffness  constants.  The 
conduction  band  is  shifted  by  the  energy 

6Ec(x,  y )  =  ac(exx  +  eyy  +  ezz)  =  2 ac ^1  - 

and  the  valence  bands  are  shifted  by 

6Ehh(x,  y)=  -  Pe-Qe 
6Eth(x,  y)=  -P(  +  Q€ 


:r 


(7) 


(8) 

(9) 


where 


P€  -  U„(exx  +  tyy  4“  t-zz) 


=  -2^ 


ii-k 


X 


Qt-  -  2^xz  eyy 

~<1+2§7> 


(10) 

(ID 

(12) 

(13) 


Fig.  1(a)  shows  a  contour  plot  relating  x  and  y  to  the  strain 
and  energy  gap  of  the  Ini.-xGaxAsyPi-_y  system  grown  on 
an  InP  substrate.  Fig.  1(b)  shows  the  same  information,  but 
inverted,  so  that,  given  a  strain  and  a  bandgap,  the  Ga  and 
As  mole  fractions  can  be  obtained.  Fig.  2(a)  and  (b)  shows 
the  results  for  the  Ini_x_yGaxAlyAs  system  with  an  InP 
substrate.  In  all  of  these  figures,  the  bandgap  with  the  smallest 
energy  is  presented  (i.e.,  c-hh  for  compressive  strain,  c-lh  for 
tensile  strain)  and  the  bandgap  shift  due  to  the  strain  effects  is 
taken  into  account.  Note  that  these  results  are  different  from 
a  conventional  plot  where  the  “unstrained”  bandgap  energy  is 
usually  presented. 


B.  Calculation  of  Band  Offset 

When  modeling  the  bandstructure  for  a  QW,  the  relative 
alignment  of  the  band  edges  of  the  well  and  barrier  material 
is  very  important.  Unfortunately,  there  are  few  data  for  arbi¬ 
trary  compositions  of  these  two  quaternary  material  systems, 
especially  for  the  strained  cases.  Here  we  present  and  compare 
two  models  for  the  band  alignment  of  these  systems. 

1)  Model-Solid  Theory  [10],  [12]:  The  valence  band  posi¬ 
tion  of  a  quaternary  is  given  by 


Ev(x ,  y)  — 


EVyav(x>  y)  4-  . 4-  8Ehh(x,  y), 

for  hh  (compressive  strain) 

Ev,av(xy  y)  H - -g—  +  8Eih(x,  y ), 

,  for  Ih  (tensile  strain) 


(16) 
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(b) 

Fig.  1.  Strain  and  strained  bandgap  energy  versus  Ga  (x)  and  As  (y)  mole 
fractions  for  the  Ini_xGaxAsyPi_y  material  system  with  (a)  contours  of 
constant  strain  and  bandgap  energy  and  (b)  contours  of  constant  x  and  y. 


where  Ev^av{x^  y)  is  the  average  valence  subband  energy  and 
A  is  the  spin-orbit  split-off  band  energy.  These  values  are 
obtained  by  a  linear  interpolation  [see  (1)  or  (2)]  of  the  binary 
values  listed  in  Table  II.  The  conduction  band  position  may 
be  calculated  by  simply  adding  the  strained  bandgap  energy 
to  the  valence  band  position 


Ec{x,  y)  = 


Ev{x,  y)  -f-  Ec—hh{Xi  y")i 
for  hh  (compressive  strain) 
Ev(x,  y)  +  Er-ih(x ,  y), 

.  for  Ih  (tensile  strain). 


(17) 


The  conduction  band  offset  is  given  by 

A  Ec  =1  K-_Et 
A  Eg  Ebg  -  E™ 


(18) 


where  Eg  and  Eb  are  the  valence  band  positions  in  the  well 
and  barrier  materials,  respectively,  and  Eg  and  Eg  are  the 
strain  adjusted  band  gaps  (E^hh  for  compressive  strain  and 
E^-ih  for  tensile  strain)  for  the  well  and  barrier  materials. 

.  2)  Harrison' s  Model  [13]:  The  position  of  both  the  con¬ 
duction  and  valence  bands  are  determined  by 


Ev{x ,  y)  =  < 


Ey(x,  y)  +  SEhh(x,  y), 
for  hh  (compressive  strain) 
E?(x,  y)  +  SElh(x,  y), 
for  Ih  (tensile  strain) 


Ec(x,  y)  =  E?  +  6Ec(x,  y) 


(19) 

(20) 


where  E^(x,  y)  and  E^(x,  y)  are  obtained  by  a  linear 
interpolation  of  the  binary  parameters  found  in  Table  H,  and 
SEhh(x,  y),  6Eik(x,  y),  and  SEc(x,  y)  are  the  strain-induced 


Fig.  2.  Strain  and  strained  bandgap  energy  versus  Ga  (x)  and  A1  (y)  mole 
fractions  for  the  Ini_x_yGaxAlyAs  material  system  with  (a)  contours  of 
constant  strain  and  bandgap  energy  and  (b)  contours  of  constant  x  and  y. 


energy  shifts  given  in  (7M9).  The  superscript  “H"  refers 
to  Harrison’s  model.  The  conduction  band-edge  discontinuity 
may  then  be  calculated  as 


A  Ec 
A  Ea 


-Ef'm 


(E^w  -  E“’b)  +  (E?'b  -  E?'w) 


pH,  b 


H,b 


-%H,  w\ 


(21) 


where  the  superscripts  w  and  b  indicate  the  well  and  barrier 
materials,  respectively. 

It  should  be  noted  that  this  method  simply  is  meant  to  give 
the  parameter  AEc/AEg  which  may  be  used  to  determine 
the  alignment  of  the  well  and  barrier  materials.  The  difference 
between  Ev(x ,  y)  and  Ec(x ,  y)  from  (19)  and  (20)  should  not 
be  used  to  calculate  the  bandgap  of  the  quaternary  material. 
Rather,  (3)  and  (4)  should  be  used. 

Fig.  3  shows  a  comparison  of  the  model-solid  theory 
and  Harrison’s  model  for  the  calculation  of  the  conduction 
band  offset  for  both  lattice-matched  Inx^xGa^ASyPi^y  and 
Ini_x_yGaxAlyAs.  Also  shown  are  empirical  curves  (solid) 
based  on  experimental  measurements  of  these  quantities 
[14],  [15].  For  the  Ini_xGaxAsyPi_y  system,  Harrison’s 
model  is  excellent  while  the  model-solid  theory  gives  the 
opposite  trend,  only  being  good  near  the  InGaAs  boundary. 
For  Ini_x_yGaxAlyAs,  both  models  have  some  error,  but 
Harrison’s  model  is  again  closer  to  the  empirical  curve.  For 
the  calculations  of  the  band  structure  in  this  paper,  we  therefore 
have  chosen  to  use  Harrison’s  model. 


C.  Calculation  of  Band  Structure  and  Gain 

The  optical  gain  with  many-body  effects  is  calculated  by  a 
non-Markovian  model  [16].  The  plasma  screening,  bandgap 
renormalization  (BGR),  and  the  excitonic  or  the  Coulomb 
enhancement  (CE)  of  the  interband  transition  probability  are 
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TABLE  D 

Parameters  for  the  Calculation  of  Band  Alignment  of  the  Ini_xGarAsyPi-y  and  Ini_x-yGarAlyAs 
Material  Systems  Using  toe  Model-Solid  Theory  [12]  and  Harrison-s  Model  [13] 


Parameter 

Symbol  (unit) 

GaAs 

InAs 

InP 

GaP 

AlAs 

Model-solid  theory: 

Average  valence  band  position 

Ev,m(eV) 

-6.92 

-6.67 

-7.04 

-7.40 

-7.49 

Spin-orbit  split-off  energy 

A(eV) 

0.34 

0.38 

0.11 

0.08 

0.28 

Harrison’s  model: 

Conduction  band  position 

E»(eV) 

1.53 

0.801 

1.35 

2.352 

2.5255 

Valence  band  position 

E»(eV) 

0.111 

0.441 

0.00 

-0.388 

-0.4245 

1.0 
0.8 

£  °’6 

3°  °*4 
0.2 

0.0 

0.0  0.2  0.4  0.6  0.8  1.0 

y 

(a) 


■  j  j  r:’-*  t  ■  » ■np-'^nrr 

- Empirical 

- Him  son's  model 


Model-Solid  Theory 


In,  GaAsP,  1 
i-x  *  y  i-y  j 


Fig.  3.  The  conduction  band  offset  for  (a)  Ini_xGax  AsyPi_y  and  (b) 
Ini_r_yGaxAlyAs  lattice-matched  to  InP.  The  solid  lines  show  empiri¬ 
cal  results  [14],  [15]  obtained  from  experimental  data.  Dashed  lines  show 
calculations  using  the  model-solid  theory  and  Harrison’s  model. 


included  in  the  model.  The  valence  band  structure  is  calculated 
by  using  the  block-diagonalized  3x3  Hamiltonian  based  on  the 
kp  method  [17].  Aside  from  the  conduction  band  effective 
mass  for  the  Ini_xGaxAsyPi_y  system,  parameters  used  in 
the  calculation  are  interpolated  using  (1)  and  (2)  and  the 
binary  parameters  summarized  in  Table  I.  The  conduction 
band  effective  mass  for  Ini_xGaxAsyPi_y  is  given  as 

Trig.  (Ini_xGaxAsyP  \  — y ) 

=  0.08  —  0.116a;  -I-  0.026?/  -  0.059 xy 
+  (0.064  -  0.02x)y2  +  (0.06  +  0.032j/)rr2 .  (22) 

The  optical  gain  spectra  are  related  to  the  spontaneous 
emission  spectra  from  the  detailed  balance  [18],  [19]  between 
absorption  and  emission  of  photons.  This  means  that  there  is 
a  transparency  point  in  the  gain  spectra  that  is  determined  by 
the  Fermi  level  separation  A F  that  suggests  the  carriers  are 
in  a  quasi-equilibrium  distribution.  The  following  fundamental 


A 


relationship  [18],  [19]  is  then  used  to  calculate  the  optical  gain: 

g{w)  =  [l  -  exp  (^~TAF)]gsPM  (23) 

9sM  =  (24) 

where  A  F  is  the  quasi-Fermi  level  separation,  k  is  Boltz¬ 
mann’s  constant,  T  is  the  device  temperature,  c  is  the  speed 
of  light,  n  is  the  effective  refractive  index  of  the  optical 
waveguide  mode,  and  rspon(ui)  is  the  spontaneous  emission 
rate.  The  spontaneous  emission  coefficient  gs p(u)  taking  into 
account  the  non-Markovian  relaxation  and  the  many-body 
effects  is  given  by  [16] 

$sPm = l  dk^wasm 

.  x(l-C(fc|,))(l-Regfcl|) 


ReL(^m(fe||))  -  Imgfc||Im£(E<m(fc||)) 
(l-Regfc|1)2  +  (Imgfc||)2 

(25) 


where  u  is  the  angular  frequency,  n0  is  the  vacuum  permeabil¬ 
ity,  e  is  the  dielectric  constant,  k\\  is  the  in-plane  wavevector, 
Lz  is  the  well  thickness,  |Mjm|2  is  the  momentum  matrix 
element  in  the  strained  QW,  ff  and  are  the  Fermi  function 
for  the  conduction  band  states  and  the  valence  band  states, 
and  ft  is  the  reduced  Planck’s  constant.  The  indices  l  and  m 
denote  the  electron  states  in  the  conduction  band  and  the  heavy 
hole  (light  hole)  subband  states  in  the  valence  band.  Also, 
Eim{k\\)  =  Ef(kn)  -  E^(kn)  +  Eg  +  AEsx  +  A£Ch  -  ^ 
is  the  renormalized  transition  energy  between  electrons  and 
holes,  where  Eg  is  the  bandgap  of  the  material,  and  AEsx 
and  A2?ch  are  the  screened  exchange  and  Coulomb-hole 
contributions  [21]  to  the  bandgap  renormalization.  The  factor 
qka  accounts  for  the  excitonic  or  Coulomb  enhancement  of 
the  interband  transition  probability  [20],  [21],  The  lineshape 
function  is  Gaussian  for  the  simplest  non-Matkovian  quantum 
kinetics  and  is  given  by  [16] 


ReL(£,„(fc,))  =  ^W%xp(.z; 


t 


(26) 


ft  *  (not  out)  - 
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TABLE  m 

A  Summary  of  the  Structure  of  the  Three  Laser  Systems  Studied 


Sample  A 

Well: 

Number  of  wells 

7 

SCH  width: 

700A 

Material 

InGaAsP 

Stripe  width: 

1.2pm 

Strain 

-0.9%  compressive 

Cavity  length: 

350pm 

Width 

80A 

PL  wavelength 

1.58pm 

Barrier: 

Material 

InGaAsP 

Strain 

lattice  matched 

Width 

80A 

PL  wavelength 

1.255 pm 

Sample  B 

Well: 

Number  of  wells 

5 

SCH  width: 

ioooA 

Material 

InGaAlAs 

Stripe  width: 

1.5pm 

Strain 

—0.78%  compressive 

Cavity  length: 

890pm 

Width 

54A 

PL  wavelength 

1.56pm 

Barrier: 

Material 

InGaAlAs 

Strain 

lattice  matched 

Width 

57A 

PL  wavelength 

1.21  /im 

Sample  C 

Well: 

Number  of  wells 

5 

SCH  width: 

600A 

Material 

InGaAlAs 

Stripe  width: 

1.46pm 

Strain 

lattice  matched 

Cavity  length: 

638pm 

Width 

86A 

PL  wavelength 

1.56pm 

Barrier: 

Material 

InGaAlAs 

Strain 

lattice  matched 

Width 

50A 

PL  wavelength 

1.21pm 

and 

Im  L(Eim(k\\))  =  J^°  exp^ 

h 

The  intraband  relaxation  time  Tm  and  correlation  time  rc  are 
assumed  to  be  constant  and  used  as  fitting  parameters. 

in.  Experiments 

A.  Laser  Structures 

In  our  experiments,  we  perform  measurements  using  lasers 
with  three  types  of  active  regions.  First,  we  tested  a  buried 
heterostructure  laser  with  five  lattice-matched  QW’s  made  of 
the  Ini_a.-yGaxAlyAs  material  system,  as  well  as  a  buried 
heterostructure  laser  with  five  -0.78%  compressively  strained 
QW’s  of  the  same  system  [5].  Also  tested  was  a  buried 
heterostructure  laser  with  seven  —0.9%  compressively  strained 
QW’s  made  of  the  Ini-xGa^AsyPx^y  material  system. 


Table  HI  contains  detailed  information  for  the  structure  of 
each  active  region. 

B.  Experimental  Procedure 

The  facet  emission  spectrum  of  each  laser  was  measured  for 
various  currents  below  lasing  threshold.  The  device  under  test 
is  mounted  on  a  heat  sink  and  is  held  at  a  constant  temperature 
of  25  °C  using  a  thermoelectric  cooling  unit.  The  current 
to  the  device  was  supplied  by  an  ILX  Lightwave  LDX-3412 
current  source  and  monitored  to  the  hundredth  of  a  milliamp 
using  a  Fluke  87  in  digital  multimeter.  The  facet  output 
was  collected  using  an  E-TEK  Dynamics  laser  optical  fiber 
interface  with  a  built-in  polarizer  to  control  the  polarization 
of  the  measured  emission.  An  optical  isolator  was  used  to 
prevent  optical  feedback  effects  and  the  light  was  monitored 
using  an  HP  7095  IB  optical  spectrum  analyzer  with  a  spectral 
resolution  of  0.08  nm.  The  net  modal  gain  for  these  devices 
was  obtained  using  the  well  known  Hakki-Paoli  method  [22] 
in  which  the  ratio  of  the  maximum  Jmax  and  minimum  Imin 
of  the  amplified  spontaneous  emission  spectrum  yields  the  net 


2nn(*||) 


^jdt. 


(27) 
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modal  gain  Gnet  by  the  following  relations: 


Gnet  —  Tg  —  Cti  —  —  In 


(28) 

(29) 


where  g  is  the  material  gain  for  one  well,  T(=NWTW)  is  the 
optical  confinement  factor  per  well  (I1*)  multiplied  by  NWJ 
the  number  of  wells,  a*  is  the  intrinsic  loss,  L  is  the  cavity 
length,  and  R\  and  R2  are  the  facet  power  reflectivities. 

The  side  light  spontaneous  emission  [24]  of  the  same  device 
was  also  monitored  for  currents  up  to  lasing  threshold.  In  this 
measurement,  a  cleaved  multimode  fiber  is  brought  up  close 
to  the  side  of  the  device,  usually  near  the  middle  in  order  to 
avoid  scattered  light  from  the  laser  facets.  First,  using  the  same 
equipment  as  for  the  gain,  the  spontaneous  emission  spectrum 
is  taken  using  a  resolution  bandwidth  of  10  nm  for  each  of 
the  currents  at  which  the  facet  emission  was  measured.  Next, 
the  current  is  controlled  using  an  HP  4145  semiconductor 
parameter  analyzer  in  order  to  achieve  currents  as  low  as 
20  /iA.  The  current  is  swept  and  the  light  captured  by  the 
fiber  is  then  fed  into  an  HP  8 153 A  Lightwave  Multimeter  to 
obtain  the  integrated  spontaneous  emission  power.  For  these 
measurements,  the  width  of  the  active  region  is  sufficiently 
small  (1.2-1.46  fim)  so  that  amplification  of  the  SE  spectra 
is  negligible. 


IV.  Analysis  and  Discussion 

The  measured  net  modal  gain  obtained  from  the  amplified 
spontaneous  emission  (ASE)  spectrum  using  (28)  for  each  of 
the  three  material  systems  studied  is  presented  in  Fig.  4  for 
several  currents  up  to  threshold.  Typical  mode  spacings  in 
the  ASE  spectrum  for  samples  A,  B,  and  C  are  about  0.93, 
0.38,  and  0.51  nm,  respectively.  To  check  for  the  effects  of 
spectrometer  resolution  on  the  measured  net  modal  gain,  the 
min/sum  method  [23],  which  is  less  sensitive  to  the  resolution 
limits  of  the  measurement  (although  more  sensitive  to  noise), 
was  also  used  to  extract  the  gain  spectra.  The  resulting  spectra 
were  nearly  identical  except  for  values  closest  to  the  threshold 
gain  value  where  the  Hakki-Paoli  method  underestimates  the 
gain,  but  these  differences  give  less  than  a  5%  error  on  the 
value  of  the  modal  gain. 

Also  presented  in  Fig.  4  are  calculated  spectra  using  the 
model  presented  above  for  the  material  gain  of  the  three 
systems.  The  intraband  relaxation  time  rin  and  the  correlation 
time  tc  used  for  the  fits  are  given  in  Table  IV.  For  the  screened 
exchange  AUsx  and  Coulomb-hole  A2?ch  contributions  to 
the  bandgap  renormalization,  [21,  eqs.  (27)  and  (28)]  were 
used  with  a  C  parameter  of  1.1,  1.7,  and  2.0  for  Samples 
A,  B,  and  C,  respectively.  Each  of  these  parameters  affect  the 
peak  gain  position  and  overall  shape  of  the  spectra.  They  were 
selected  to  yield  the  best  overall  fit,  balancing  these  criteria. 
The  net  modal  gain  is  obtained  using  an  intrinsic  loss  which 
is  independent  of  wavelength,  and  a  confinement  factor  with 
a  linear  dependence  on  wavelength,  with  a  slope  equal  to 
1.66xl0~4  nm-1  (estimated  from  the  wavelength-dependent 
refractive  index  and  waveguide  theory  for  the  fundamental 
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Fig.  4.  The  measured  modal  gain  (solid  lines)  and  calculated  (dashed  lines) 
values  from  the  electronic  band  structure  for  (a)  Ini_zGazAsyPi-y  with 
—0.9%  compressive  strain  (Sample  A),  (b)  Ini_z-yGaxAlyAs  with  —0.78% 
compressive  strain  (Sample  B),  and  (c)  Im_z_yGazAlyAs  lattice  matched 
(Sample  C).  The  device  parameters  are  listed  in  Table  HI. 


mode)  and  the  vertical  intercept  treated  as  a  parameter  to  fit 
the  data. 

Fig.  5(a)  shows  a  plot  of  the  fitted  values  of  the  intrinsic 
loss  as  a  function  of  the  injected  current  for  each  device. 
The  observed  increase  in  the  loss  is  consistent  with  previous 
reports  [25]  and  can  be  explained  by  the  increase  in  free- 
carrier  absorption  with  an  increasing  carrier  density.  The  fitted 
optical  confinement  factor  used  for  the  three  devices  is  shown 
in  Fig.  5(b).  Values  for  the  optical  confinement  factor  per  well 
(Tw)  and  the  overall  confinement  factor  (T  =  NWTW)  are 
listed  in  Table  IV.  These  fitted  values  are  in  good  agreement 
with  typical  values  for  InP-based  lasers  [26]  when  scaling  due 
to  well  width  is  taken  into  account.  Finally,  Fig.  5(c)  shows  a 
plot  of  peak  net  modal  gain  versus  the  total  current  density  for 
each  laser  system.  The  fits  of  the  gain  spectra  of  these  three 
material  systems  allow  us  to  obtain  the  density  of  carriers  N 
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TABLE  IV 

Parameters  Used  in  and  Obtained  from  the  Theoretical  Tits 


Parameter 

Symbol 

Sample  A 

- ^ - 

Sample  B 

Sample  C 

Fit  of  the  modal  gain: 

Interband  relaxation  time 

n.(/*) 

40 

60 

65 

Correlation  time 

»;(/') 

20 

20 

20 

Quasi-Fermi  level  separation 

AF(eV) 

0.801  (4mA) 

0.804  (5mA) 

0.806  (6mA) 

0.813  (8mA) 

0.815  (7mA) 

0.816  (8mA) 

0.822  (8mA) 

0.823  (9mA) 

0.823  (10mA) 

0.828  (10mA) 

0.829  (12mA) 

0.827  (11.7mA) 

Transformation  of  spontaneous  emission: 

Scaling  coefficient 

6.11  x  10”T  dk  0.22  x  10“T 

1.61  x  10-»  dk  0.08  x  10-T 

6.2x10“*  ±0.02x10-* 

Effective  index  of  refraction 

n 

3.3  ±  0.15 

3.3  ±0.15 

3.3  ±0.15 

Optical  confinement  factor 

r(«  NJtm) 

0.15  A  0.02 

0.06  ±0.02 

0.11  ±0.02 

Optical  confinement  factor  per  well 

r„ 

0.021  ±  0.0029 

0.012  ±0.004 

0.022  ±0.004 

Fit  of  the  J,N  relationship: 

Coupling  coefficient 

1.58  x  10-“  ±  0.01  xl0’M 

1.64x10-“  ±0.05x10-” 

1.21  x  10-“  ±  0.01  x  10‘“ 

Monomolecular  coefficient 

Mr') 

1.1  x  10*  A  0.2  x  10* 

7  x  10T  ±  2  x  10r 

6.7  x  10r  ±  3  x  10r 

Auger  coefficient 

C(em**-') 

1.1x10-”  ±0.15x10'” 

1.4x10-”  ±0.3x10-" 

3x10-”  ±2x10-” 

Injection  efficiency 

*•> 

0.92  ±0.17 

0.95  ±  0.15 

0.45  ±0.15 

Characteristic  temperature  (measured] 

T.(‘K) 

38 

53 

51.4 

in  the  QW’s  corresponding  to  the  measured  injection  currents. 
For  each  system,  other  samples  from  the  same  wafer  exhibit 
similar  gain  spectra. 

The  same  model  used  to  create  the  gain  spectrum  was  also 
used  to  calculate  the  spontaneous  emission  rate  for  each  of  the 
material  systems  studied.  The  radiative  current  density  can  be 
related  to  the  spontaneous  emission  rate  as 


roo 

J rad  —  I  ^spon(^)  ^  —  eiV wLwBefiN 

Jo 


(30) 


where  e  is  the  charge  on  an  electron,  Nw  is  the  number  of 
wells,  Lw  is  the  width  of  a  QW,  rspon(A)  is  the  spontaneous 
emission  rate,  Beff  is  the  effective  radiative  recombination 
coefficient,  and  N  is  the  carrier  density.  Jraa  can  be  obtained 
from  our  calculated  spontaneous  emission  spectrum  using  (30) 
with 

fspon(A)  =  ^f-ffsp(A)  (31) 


where  gsp( A)  is  given  in  (25),  A  =  2kc/uj,  c  is  the  speed  of 
light,  and  n  is  the  refractive  index  of  the  QW.  The  calculated 
curves  for  Jrad  of  each  sample  are  shown  in  Fig.  6. 

In  order  to  validate  these  calculations  from  the  band  struc¬ 
ture,  it  is  important  to  be  able  to  obtain  the  radiative  current 
density  experimentally.  As  described  above,  we  measure  the 
spontaneous  emission  power  spectra  from  the  side  of  each 
laser  device.  The  spontaneous  emission  power  captured  by  the 
optical  fiber  can  be  related  to  the  spontaneous  emission  rate  by 

he 

P(X)  =  Ki—rBpon(X)  d\V  (32) 


where  P(  A)  is  the  measured  power  in  the  fiber,  Ki  is 
a  coupling  coefficient  which  includes  the  fraction  of  the 
spontaneous  emission  captured  by  the  optical  fiber,  h  is 


Planck’s  constant,  c  is  the  speed  of  light,  d\  is  the  spectrometer 
resolution,  and  V  is  the  active  region  volume.  Since  K\  is  an 
unknown  quantity,  it  is  not  straightforward  to  obtain  a  value 
for  the  spontaneous  emission  rate  so  that  (30)  can  be  used  to 
extract  Jrac i.  A  calibration  can  be  accomplished  based  on  the 
fundamental  relationship  between  the  material  gain  and  the 
spontaneous  emission  rate  [(23)  and  (24)],  rewritten  as 


3(A)  = 


A4 

87 rcn2 


1  —  exp 


/  hc/X  -  AF 

V  kBT 


(33) 


where  n  is  the  effective  refractive  index,  A F  is  the  quasi- 
Fermi  level  separation  energy,  kB  is  Boltzmann’s  constant, 
and  T  is  the  temperature.  The  material  gain  per  well  can  be 
related  to  the  measured  net  modal  gain  as 


G«t  =  rfl(A)-Oi  (34) 


where  T  is  the  optical  confinement  factor  and  a*  is  the  intrinsic 
loss  for  the  optical  waveguide  mode.  Combining  (32M34), 
we  obtain 

Gnet(A)  =  QA5  (l  -  exp  />(  A)  _  a,  (35) 

where 

r 

Q  =  8xKlhc2n2  dXV'  (36) 

This  equation  relates  the  measured  net  modal  gain  to  the 
measured  spontaneous  emission  power  using  the  three  pa¬ 
rameters  Q,  A  F,  and  or<.  A  F  and  a,-  are  not  independent, 
however,  since  the  net  modal  gain  value  measured  at  the 
wavelength  corresponding  to  A  F  is  exactly  the  transparency 
level  and  equal  to  cti.  (We  assume  here  that  the  intrinsic  loss  is 
independent  of  wavelength.)  Hence,  the  short-wavelength  side 
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Wavelength  (nm) 
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(c) 

Fig.  5.  (a)  The  optical  confinement  factor  and  (b)  intrinsic  loss  used  in  the 
fits  of  the  modal  gain  for  Samples  A,  B,  and  C  in  Fig.  4.  (c)  The  peak  net 
modal  gain  versus  the  total  current  density. 


Wavelength  (nm) 
(a) 


Wavelength  (nm) 

(b) 

Fig.  7.  (a)  The  measured  spontaneous  emission  spectra  for  the  —0.78%  com- 
pressively  strained  Ini_x-yGar  AlyAs  laser  system  (Sample  B  in  Table  HI) 
and  (b)  the  measured  net  modal  gain  spectra  (solid  lines)  compared  to  the 
transformed  spontaneous  emission  spectra  (dashed  lines)  using  (35). 

into  the  net  modal  gain  spectrum.  By  choosing  the  best 
fit  between  the  measured  and  transformed  net  modal  gains, 
an  accurate  value  for  Q  can  be  obtained.  As  an  example 
of  this  procedure.  Fig.  7(a)  shows  the  spontaneous  emis¬ 
sion  power  spectra  measured  for  the  compressively  strained 
Ini_x_yGaxAlyAs  system.  These  spectra  were  measured  for 
each  of  the  currents  that  the  gain  was  measured  at.  The  best  fits 
between  the  transformed  spontaneous  emission  and  actually 
measured  net  modal  gain  spectra  are  shown  in  Fig.  7(b).  The 
extracted  values  of  Q  are  about  1.8xl0~7  and  vary  by  less 
than  4%  among  all  of  the  currents  measured,  showing  excellent 
uniformity  in  the  measurement  of  the  SE  spectra.  For  AF, 
the  extracted  values  for  each  system  deviate  from  the  values 
obtained  for  the  fitting  of  the  experimental  gain  curves  by  less 
than  0.1%  at  low  currents  where  the  difference  between  two 
adjacent  curves  in  Fig.  4  represent  a  shift  of  about  1%  in  A F. 
At  higher  currents,  the  deviation  gets  as  large  as  0.45%  while 
the  separation  between  curves  is  about  0.5%. 

Once  accurate  values  of  Q  have  been  extracted,  the  mea¬ 
sured  spontaneous  emission  power  spectra  can  be  calibrated 
in  the  following  way  to  get  the  actual  spontaneous  emission 
rate.  By  solving  (32),  rspon(A)  can  be  expressed  as 


Fig.  6.  The  calculated  total  radiative  current  density  through  all  wells  [see 
(30)]  as  a  function  of  carrier  density  (lines)  compared  to  the  extracted  total 
radiative  current  density  from  the  calibrated  spontaneous  emission  spectrum 
measurements. 

of  the  measured  net  modal  gain  spectrum  gives  the  relationship 
between  these  two  quantities. 

Using  Q  and  A F  as  fitting  parameters,  (35)  may  then  be 
used  to  transform  the  measured  spontaneous  emission  power 


rspon(A)  =  TAP(A)  (37) 

where 

(38) 

Note  that  this  procedure  has  eliminated  the  need  to  derive 
the  coupling  coefficient  K\  explicitly.  The  data  points  in  Fig.  6 
were  obtained  by  substituting  (37)  and  (38)  into  (30)  and 
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Fig.  8.  The  radiative  recombination  coefficient  as  a  function  of  carrier 
density  extracted  from  the  calculated  total  radiative  current  density. 

numerically  integrating  the  calibrated  data.  The  value  for  n 
was  assumed  to  be  3.3  with  an  error  of  0.15  for  each  of  the 
material  systems  measured.  The  optical  confinement  factor  was 
assumed  to  be  constant  with  a  value  equal  to  the  average  of  the 
range  used  for  the  gain  fits,  as  shown  in  Fig.  5(b).  The  error 
range  is  taken  from  the  values  of  the  confinement  factor  in 
Fig.  5(b)  at  the  edges  of  the  transformed  spontaneous  emission 
spectra.  This  amounts  to  an  error  of  about  ±12%  for  the  value 
of  the  confinement  factor,  leading  to  a  total  error  of  slightly 
over  ±20%  for  the  radiative  current  density. 

From  the  calculated  curves  of  the  radiative  current  den¬ 
sity  Jra d  as  a  function  of  carrier  density  N  in  Fig.  6,  the 
radiative  recombination  coefficient  Beff  for  each  material 
system  is  obtained  from  (30)  and  is  presented  in  Fig.  8.  This 
parameter  is  independent  of  the  well  width/number  product, 
so  it  serves  as  a  fair  comparison  of  the  relative  probabil¬ 
ity  of  the  radiative  recombination  of  an  electron-hole  pair 
for  each  structure.  We  see  that,  at  low  carrier  densities, 
the  compressively  strained  Ini_xGaxAsyPi_y  system  has  the 
largest  radiative  recombination  coefficient  followed  by  the 
compressively  strained  Ini_x_yGaxAlyAs  system  and  then  the 
lattice-matched  Ini_x_yGaxAlyAs  system.  It  is  not  surprising 
that  the  strained  samples  are  more  efficient  since  strain  has 
been  shown  to  decrease  the  in-plane  effective  mass  in  the 
top  (hh)  valence  band,  leading  to  a  more  symmetric  injection 
of  electrons  and  holes  in  the  conduction  and  valence  bands 
[27].  Also,  the  heavy-hole  and  light-hole  subband  separation 
becomes  larger  with  strain,  leading  to  less  wasted  carriers. 
At  higher  carrier  densities,  all  three  systems  show  a  decrease 
in  the  radiative  recombination  coefficient.  This  is  especially 
pronounced  in  the  compressive  Ini_xGaxAsyPi_y  system. 
This  decrease  is  probably  due  to  the  significant  population 
of  higher  energy  valence  subbands  while  the  higher  energy 
conduction  subbands  are  virtually  empty.  For  example,  in  all  of 
these  systems,  the  second  conduction  subband  has  less  than  3% 
of  the  total  carrier  population.  Meanwhile,  for  the  compressive 
Ini_xGaxASyPi_y  system,  the  first  two  heavy-hole  subbands 
(which  are  the  lowest  subbands)  have  a  spacing  of  36  meV 
and  a  distribution  of  approximately  60.1%  of  carriers  in  the 
first  heavy-hole  subband  and  29.6%  in  the  second  heavy- 
hole  subband.  Also,  for  the  compressive  Ini_x_yGaxAlyAs 
sample,  there  is  a  spacing  of  65  meV  between  the  first  two 
heavy-hole  subbands  and  a  carrier  distribution  of  about  80% 
and  19.4%  for  the  first  and  second  heavy-hole  subbands, 
respectively.  This  difference  is  mostly  due  to  the  thinner 


Fig.  9.  The  square  root  of  the  measured  integrated  spontaneous  emission 
power  >/P  (plotted  along  the  x  axis)  versus  the  measured  injection  current 
along  the  vertical  axis.  The  x  axis  is  roughly  proportional  to  carrier  density, 
and  the  slope  of  the  lines  on  the  log-log  plot  give  the  power  law  for  the 
injection  current  and  carrier  density  relation. 

well  of  the  compressive  Ini_x_yGaxAlyAs  sample.  Holes 
in  the  second  heavy-hole  subband  are,  therefore,  unable  to 
significantly  contribute  to  radiative  recombination  processes 
since  they  can  only  recombine  with  electrons  in  the  second 
conduction  subband,  which  is  sparsely  populated  for  each 
laser  system  studied.  The  decrease  seen  in  the  compressive 
Ini_xGaxAsyPi_y  system  is  more  significant  due  to  the 
smaller  spacing  between  valence  subbands.  At  the  higher 
carrier  densities,  this  spacing  leads  to  an  even  larger  increase 
in  the  population  of  higher  energy  subbands,  and  thus  more 
wasted  carriers. 

To  study  the  major  recombination  mechanisms  in  each  of 
these  lasers,  it  is  convenient  to  use  the  measured  spontaneous 
power  which  has  been  monitored  by  an  optical  power  meter. 
In  this  case,  the  data  are  easily  measured  over  a  large  current 
range,  and  the  photodetector  performs  the  integration  of  the 
spontaneous  emission  power  spectrum.  Fig.  9  shows  a  log-log 
plot  of  the  square  root  of  spontaneous  emission  power  versus 
the  total  injected  current  for  each  device  studied.  This  format 
was  chosen  [28]  to  give  an  indication  of  the  relative  con¬ 
tributions  of  each  recombination  mechanism  (monomolecular, 
radiative,  and  Auger).  The  value  of  the  x  axis  can  be  expressed 
as 

\fP  =  K2y/KsN  (39) 

where  K2  is  another  unknown  coefficient  which  includes 
the  fraction  of  spontaneous  emission  power  captured  by  the 
optical  fiber.  We  see  that  this  is  almost  directly  proportional 
to  Ny  but  not  perfectly,  due  to  the  carrier  dependence  of 
Beff,  as  discussed  above.  This  deviation  is  relatively  small, 
however,  since  we  can  see  clear  regions  of  slope  one  and  slope 
three  which  correspond,  respectively,  to  regimes  where  I  is 
proportional  to  N  (monomolecular  mechanisms  dominate)  and 
Nz  (where  Auger  mechanisms  dominate).  The  region  of  slope 
two  could  indicate  that  the  radiative  processes  are  dominant, 
but  this  is  not  clear  since  there  would  be  a  point  at  which 
the  slope  is  two  due  to  the  transition  between  the  N  and  N 3 
regimes,  even  if  radiative  processes  were  missing.  These  issues 
will  be  clarified  in  the  following  analysis. 

It  is  necessary  to  transform  the  measured  data  into  a 
relationship  between  total  current  density  and  carrier  density. 
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(a) 


(b) 

Fig.  10.  The  extracted  current  density/carrier  density  relation  (dashed 
lines)  as  well  as  the  best  fits  from  the  theory  (solid  lines)  are  shown 
for  (a)  the  Ini_xGarAsyPi_y  laser  system  (Sample  A)  and  (b)  the  two 
Ini_z-vGazAlyAs  systems  (Samples  B  and  C). 


Current  density  is  easily  obtained  from  the  current  as 


(40) 


where  w  is  the  stripe  width  and  L  is  the  cavity  length.  To 
obtain  N ,  however,  the  unknown  coefficient  K2  must  be 
determined.  The  fits  of  the  net  modal  gain  with  our  theory 
give  a  relation  between  I  and  N  over  a  limited  range  of 
the  measured  data  presently  under  consideration.  Using  the 
calculated  dependence  of  £eff  on  the  carrier  density  for  each 
system,  and  the  measured  \[P  value  corresponding  to  the 
currents  for  the  fitted  gain  data,  (39)  yields  a  value  of  K2 
for  each  material  system.  The  average  between  the  values 
obtained  for  each  of  the  I-N  pairs  are  used  for  our  calculations 
and  are  listed  in  Table  IV.  All  values  are  within  2%  of  the 
average  value.  Since  K2  should  be  independent  (for  the  most 
part)  of  injection  current,  it  is  fair  to  extrapolate  the  value 
of  K2  from  this  small  range  of  currents  to  the  entire  range 
of  data  under  consideration.  Using  these  K2  values  and  the 
calculated  Beff(N),  it  is  straightforward  to  obtain  a  value  for 
N  corresponding  to  the  measured  y/P  for  each  device. 

The  extracted  relationships  between  Jto t  and  N  are  shown 
as  dashed  lines  in  Fig.  10(a)  for  the  Ini_xGaxAsyPi_y  system 
and  in  Fig.  10(b)  for  the  two  Ini-x-yGaxAlyAs  systems.  The 
solid  lines  in  these  plots  correspond  to  the  best  fit  results  using 
the  relation 

Jtot  =  ^^{AN  +  Beff(N)N2  +  CN 3]  (41) 

where  r)inj  is  the  injection  efficiency  accounting  for  the 
fraction  of  carriers  which  make  it  into  the  QW’s,  Nw  is 
the  number  of  wells,  A  is  the  monomolecular  recombination 
coefficient,  Beft{N)  is  the  calculated  radiative  recombination 


coefficient,  and  C  is  the  Auger  recombination  coefficient.  A  is 
a  fitting  parameter,  but  is  only  important  for  very  low  current 
densities.  Its  value  is  determined  by  fitting  the  behavior  of 
the  curve  at  these  low  current  densities.  7?inj  and  C  are  the 
major  parameters  which  control  the  behavior  of  the  curve.  For 
a  set  value  of  7?inj,  the  penalty  function  for  the  fit  is  very 
sensitive  to  the  C  coefficient,  so  the  value  can  be  very  well 
determined.  However,  a  range  of  7fcnj  will  give  decent  fits,  so 
the  extracted  C  value  has  some  uncertainty.  The  best  fit  values 
of  7jinj,  Ay  and  C  are  presented  in  Table  IV.  The  uncertainties 
on  rjinj  and  C  are  determined  by  selecting  those  values  which 
raise  the  penalty  function  20%  over  the  best  fit  value.  The 
injection  efficiency  for  the  unstrained  Sample  C  appears  to 
have  a  relatively  low  injection  efficiency. 

It  should  be  noted  that  previous  measurements  of  the  Auger 
coefficient  in  the  Ini-.xGaxAsyPi_y  material  system  have 
given  values  between  2xl0~29cm6s~1  and  lxl0~28  cm6s^1 
[29]— [31].  As  discussed  previously,  the  fact  that  a  signifi¬ 
cant  amount  of  holes  occupy  states  where  they  may  only 
recombine  by  nonradiative  processes  could  explain  why  the 
extracted  Auger  coefficient  for  our  sample  is  significantly 
higher  than  this  range.  This  indicates  that  the  Auger  coefficient 
can  be  minimized  by  selecting  the  strain  and  well  width 
to  increase  the  separation  of  the  valence  subbands.  For  the 
Ini_x_yGaxAlyAs  material  system,  a  reported  value  of  the 
Auger  coefficient  [9]  is  3.6xl0~29  cm6s“1  which  is  smaller 
than  our  value  for  the  compressive  sample.  Our  procedure, 
based  on  a  comprehensive  gain  model  and  a  fitting  of  the  gain 
spectra  to  determine  N  and  Beff,  should  provide  a  reliable 
estimate  of  the  C  coefficient. 

These  extracted  values  are  used  in  Fig.  11  to  com¬ 
pare  the  relative  contributions  of  Jmo no*  Aa.  and  Aug 
for  the  compressive  Ini_x_yGaxAlyAs  and  compressive 
Ini_xGaxAsyPi_y  systems.  Each  device  is  dominated  by  the 
Auger  recombination  well  before  the  laser  reaches  threshold. 
However,  due  to  an  Auger  coefficient  about  10  times  larger, 
the  Auger  recombination  dominates  at  a  much  lower  point 
and  becomes  much  larger  in  the  Ini_xGaxAsyPi_y  system. 
This  leads  to  a  much  larger  threshold  current  density.  The 
To  for  each  of  these  lasers,  measured  from  the  ( shift  in 
threshold  over  a  temperature  range  from  15°  to  60°,  is  also 
presented  in  Table  IV.  These  relative  values  agree  with  the 
trend  in  the  extracted  Auger  coefficients,  indicating  that 
Auger  recombination  is  a  major  factor  reducing  To  for 
these  laser  systems  [28],  [32].  Finally,  it  is  interesting  to 
note  that  the  radiative  current  density  never  becomes  the 
dominant  recombination  mechanism.  At  the  injection  level 
where  Auger  recombination  becomes  significant,  the  radiative 
current  density  is  about  the  same  as  the  monomolecular  current 
density  for  each  system  studied. 

V.  Conclusion 

In  this  paper,  we  have  presented  a  comprehensive  model  for 
the  calculation  of  the  band  edges  of  both  Ini_xGaxAsyPi_y 
and  Ini-x~yGaxAlyAs  QW  lasers  with  an  arbitrary  com¬ 
position.  Using  a  many-body  optical  gain  model,  we  have 
compared  the  net  modal  gain  for  both  material  systems  with 
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Fig.  11.  The  relative  contributions  of  the  monomolecular,  radiative, 
and  Auger  recombination  currents  for  (a)  —0.9%  compressively  strained 
Ini_xGarAsyPi_y  (Sample  A),  (b)  —0.78%  compressively  strained 
Ini_x-yGaxAlyAs  (Sample  B),  and  (c)  lattice-matched  Ini_x_yGaxAlyAs 
(Sample  C). 


calculations  from  the  realistic  band  structure.  These  compar¬ 
isons  allow  us  to  obtain  carrier  densities  for  each  of  the 
measured  currents.  From  a  calibrated  measurement  of  the 
spontaneous  emission  spectrum,  we  were  able  to  obtain  the 
radiative  current  density.  These  results  are  all  in  excellent 
agreement  with  our  calculations.  From  a  fit  of  the  extracted 
relationship  between  carrier  density  and  total  current  density, 
we  were  able  to  extract  values  for  the  Auger  coefficient  of 
each  material  system  studied.  The  relative  values  of  these  pa¬ 
rameters  correspond  very  well  with  the  observed  temperature 
performance  of  each  device. 

From  this  study,  it  is  difficult  to  make  a  strong  statement  as 
to  which  of  the  two  material  systems  is  superior.  Obviously, 
more  studies  on  a  variety  of  active  region  designs  would  be 
needed  to  gauge  all  of  the  relevant  factors  and  draw  such 
a  conclusion.  Our  measurements  seem  to  indicate  that  the 


InGaAlAs  system  is  at  least  as  good  as  and  probably  better 
than  the  InGaAsP  systems  in  terms  of  Auger  recombination 
and  temperature  stability.  Also,  care  should  be  taken  to  design 
the  QW’s  such  that  there  is  a  large  separation  between  the 
lowest  and  highest  level  subbands,  in  order  to  reduce  the 
carriers  wasted  by  nonradiative  recombination  and  to  reduce 
the  temperature  sensitivity  of  the  lasers. 
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Abstract — A  distributed-feedback  (DFB)  laser  with  a  curved 
waveguide  in  the  active  region  integrated  with  an  electroabsorp¬ 
tion  modulator  is  studied  experimentally  and  theoretically.  The 
modulator  controls  the  lasing  wavelength  of  the  DFB  laser  by 
acting  as  an  optical  phase  shifter.  Our  model,  which  is  based  on 
the  transfer  matrix  method,  is  used  to  simulate  this  multisection 
device  with  the  curved  waveguide,  self-consistently  including  the 
effects  of  spatial  hole  burning  (SHB).  The  model  explains  the 
features  and  wavelength-tuning  behavior  of  the  spectrum  and 
shows  good  agreement  with  experimentally  measured  spectra.  We 
also  show  theoretically  that  the  curved  waveguide  suppresses  the 
longitudinal  photon  density  profile  of  the  device  compared  with 
a  straight  waveguide  case,  which  implies  reduced  SHB  effects. 

Index  Terms — Curved  waveguide,  distributed  feedback,  mod¬ 
ulator,  spatial  hole  burning,  wavelength-tunable. 


I.  Introduction 

THE  STUDY  OF  wavelength-tunable  laser  sources  is  an 
important  topic  of  research,  emphasized  by  the  rapidly 
expanding  development  of  wavelength-division-multiplexed 
(WDM)  optical  communications  systems,  which  can  presently 
transmit  data  at  rates  of  over  1  Tb/s  over  a  single  optical  fiber. 

Research  in  multisection  wavelength-tunable  lasers  has  been 
investigated  extensively  for  use  in  WDM  and  coherent  optical 
communications  systems,  and  some  of  the  progress  in  this 
field  is  well  summarized  in  [1]  and  [2].  Originally,  the  goal 
of  developing  distributed-feedback  (DFB)  lasers  integrated 
with  electroabsorption  modulators  was  to  improve  high-speed 
performance  by  reducing  frequency  chirping  of  DFB  lasers  by 
external  modulation  [3]-[9].  Recently,  these  devices  have  also 
been  used  as  wavelength-tunable  sources  [10],  [11].  The  use 
of  a  chirped  grating  pitch  [12]— [14]  has  been  investigated  as  a 
way  to  achieve  a  distributed  phase  shift  in  the  grating.  Several 
groups  have  implemented  this  chirped  grating  pitch  using 
a  bent  or  curved  waveguide  to  achieve  wavelength  tuning 
or  high-power  operation  and  have  shown  promising  results 
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[15]— [17].  In  previous  work,  we  have  modeled  a  DFB  laser  and 
an  integrated  electroabsorption  modulator-laser  (EML)  using 
the  transfer  matrix  method  [18],  [19].  Multisection  DFB  lasers 
have  been  modeled  using  a  transmission-line  model  as  well 
[20]. 

The  wavelength-tunable  laser  studied  here,  a  possible  candi¬ 
date  for  WDM  systems,  consists  of  a  DFB  laser  with  a  curved 
waveguide,  integrated  with  an  electroabsorption  modulator 
which  acts  as  a  wavelength-tuning  section.  The  structure  for 
this  device  is  identical  to  that  of  an  integrated  EML  with  the 
exception  that  the  facet  coatings  on  the  laser  and  modulator 
are  opposite  to  those  of  the  conventional  EML.  A  voltage 
bias  applied  to  the  modulator  effectively  controls  the  lasing 
wavelength  of  the  DFB  laser.  Previously  reported  results  for 
this  device  have  shown  a  wavelength-tuning  range  of  3.5  nm 
and  high-speed  optical  packet  switching  over  four  wavelength 
channels  at  a  rate  of  2.5  Gb/s  [11]. 

The  purpose  of  this  work  is  to  model  this  wavelength- 
tunable  multisection  device,  including  spatial  hole  burning 
(SHB)  effects  self-consistently,  and  to  explain  the  charac¬ 
teristics  of  the  spectrum  and  the  wavelength-tuning  behavior 
of  this  device  by  showing  good  agreement  between  theory 
and  experiment.  We  first  describe  the  device  structure  in 
Section  II.  In  Section  III,  we  explain  our  model  based  on  the 
transfer  matrix  method  that  includes  SHB  effects.  In  Section 
IV,  we  compare  measured  spectrum  with  theory  and  explain 
the  spectral  features  of  this  wavelength-tunable  laser.  This 
includes  explanations  for  the  wavelength-tuning  behavior  of 
this  device  as  a  function  of  modulator  voltage  and  theoretical 
results  which  show  SHB  suppression  due  to  the  use  of  a  curved 
waveguide.  In  Section  V,  we  summarize  our  work  and  present 
our  conclusions. 

n.  Device  Structure 

A  schematic  of  the  wavelength-tunable  laser,  which  consists 
of  a  DFB  laser  with  a  curved  waveguide  integrated  with  an 
electroabsorption  modulator  separated  by  an  isolation  region, 
is  shown  in  Fig.  1.  The  length  of  the  DFB  section  is  300  jum, 
the  isolation  region  is  about  80  /zm,  and  the  modulator  is  250 
fj,m.  The  facet  reflectivity  at  the  high-reflection  (HR)-coated 
facet  on  the  modulator  side  is  99%,  and  the  antireflection  (AR)- 
coated  facet  reflectivity  on  the  DFB  laser  side  is  estimated  to 
be  1%  using  a  50-jzm  window  structure  on  the  AR  side  of  the 
DFB  section.  The  shape  of  this  curved  waveguide  has  been 
modeled  previously  as  a  raised  sine  shape  [11].  However,  this 
waveguide  can  also  be  modeled  as  an  S-bend  shape  which 
has  more  flexible  parameters  and  provides  a  more  general 
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(b) 

Fig.  1.  (a)  A  side-view  schematic  of  the  two-section  wavelength -tun able 

laser  which  consists  of  a  DFB  laser  section  at  the  AR  side  and  a  modulator 
at  the  HR  side  of  the  device  which  tunes  the  lasing  wavelength.  Ltota \  is 
the  length  of  the  entire  device,  Ldfb»  £iso*  and  Lrnod  are  the  lengths  of 
the  DFB  laser,  isolation,  and  modulator  sections,  respectively,  Z]  =  Lj) fb< 
and  h  =  Ldfb  +  ^iso-  (b)  A  top-view  schematic  illustrating  the  curved 
waveguide  of  the  DFB  laser  section. 

waveguide  shape.  The  shape  of  the  S-bend  curved  waveguide 
is  given  by  the  following  equations  [16]: 


The  laser,  isolation,  and  modulator  sections  can  each  be 
represented  by  the  transfer  matrices  Fdfb,  F “o1,  and  Fmod, 
respectively,  which  gives  the  slowly  varying  amplitudes  of 
the  forward  and  backward  propagating  waves  A(z)  and  B(z ), 
respectively,  in  each  section.  The  elements  of  these  matrices 


are  given  by  [18],  [19] 

Fn(l)  =  V^r(ciql-rprme-^) 

J-  “  '  p1  m 

(2) 

1  ~  "  v"  m 

(3) 

W)=rJfr- 

J-  “  f  pf  m 

(4) 

F22Q)  —  1  (  rprme  +e  ) 

J-  •  p*  m 

(5) 

where 

q—\/  (A/3)2  +  KabKba 

(6) 

and 

Kba 

rp~A0  +  q 

K ab 

rm~A(3  +  q 

(7) 

(8) 

y(z)=2W,  for  0  <  2  <  Z/i 

W  W  (  z-L\\ 

y{z)  =2 W  +  -(L1-z)  +  -  sia(n  —  j, 

for  Li  <  z  <  Li  4-  2Lb 

y(z)  =  0,  for  Lx  +  2LB  <  z  <  L  (1) 

where  L  is  the  total  length,  Lx  is  the  length  of  the  straight 
section  before  the  S-bend,  Lb  is  equal  to  half  of  the  length 
of  the  S-bend  region,  and  W  is  the  width  of  half  of  the  S- 
bend.  For  this  device,  we  estimated  the  S-bend  parameters  to 
be  L  =  360  fjm,  Lx  =  47  ym,  Lb  =  133  /um,  and  W  =  10 

fim. 

In  both  devices,  the  active  region  in  the  DFB  and  the 
modulator  sections  consist  of  seven  layers  of  compressively 
strained  InGaAsP  QW’s  and  were  fabricated  using  a  selective- 
area  metal-organic  vapor  phase  epitaxy  (MOVPE)  technique. 
The  DFB  region  has  an  estimated  compressive  strain  of  at 
least  —1%  while  the  compressive  strain  in  the  modulator  was 
measured  to  be  about  -0.5%.  The  peak  wavelength  of  the 
photoluminescence  was  measured  at  1.4891  fjm.  The  isolation 
region  between  the  DFB  and  modulator  sections  was  created 
by  shallow  trenching  a  depth  of  about  0.7  ym  to  yield  an 
electrical  isolation  with  a  typical  value  greater  than  50  kO. 

in.  Theory 

A.  Longitudinal  Fields  in  a  Multisection  Device 

We  model  the  longitudinal  properties  of  the  optical  fields 
in  the  curved  waveguide  DFB  laser  with  integrated  modulator 
and  isolation  sections  using  the  transfer  matrix  method.  Most 
of  the  development  of  the  transfer  matrix  method  used  in  our 
model  can  be  found  in  [18]  and  [19]. 


where  [3  =  27rne(A)/A  is  the  propagation  constant  where  ne  is 
the  effective  index  of  the  fundamental  guided  mode,  (30  =  7r/A 
where  A  is  the  grating  pitch,  and  the  detuning  parameter 
A fi  =  p-  0o. 

The  interface  between  sections  which  includes  differences 
in  refractive  index  between  each  section  is  also  represented 
by  a  transfer  matrix  T 1x11  and  is  derived  from  the  boundary 
conditions  that  tangential  electrical  and  magnetic  fields  are 
continuous  at  the  interface  at  a  given  2  =  lx  or  I2,  given  by 


pw+r 

M+). 


=  T'nt(l) 


(9) 


The  forward  and  backward  propagating  waves  in  the  whole 
device  can  be  modeled  by  multiplying  the  matrices  for  each 
section  and  interface  together  as  given  by 


b{l“T1  =  [r“M 

.  p«(il)]  ^(LDFB)]  [£<“>]  (10) 


where  Ltotai  is  *he  length  of  the  entire  device,  Ldfb>  ^iso,  and 
Lmod  are  the  lengths  of  the  DFB  laser,  isolation,  and  modulator 
sections,  respectively,  and  k  is  the  location  of  section  i,  where 
h  =  LDFb,  h  =  LDfb  +  £iso>  as  shown  in  Fig.  1. 

We  model  the  amplified  spontaneous  emission  spectrum  for 
the  integrated  device  by  including  an  equivalent  spontaneous 
emission  source  [21]  into  the  transfer  matrix  method.  We  place 
an  equivalent  spontaneous-emission  surface  current  density  Js 
at  a  position  in  the  (active)  DFB  section  of  the  cavity  and 
solve  the  boundary  conditions  for  the  electric  and  magnetic 
fields  on  either  side  of  the  current  sheet.  From  this,  we  obtain 
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the  relation  between  the  field  amplitudes  near  the  source 
position 


'A+(zs) ' 
B+(zs)_ 


'(*«) 

'(*•) 


=  M 


-1 


0 

U)pJ8 


521 

522 


(ID 


where 

3?  = 


eifioZs 
\j3ei0OZs  +  KbaQ-ipoz 


o-iflaz. 


Kabeif3oZ’  -  \ 


(12) 


where  nab  and  Kia  are  the  coupling  coefficients  used  in 
the  coupled-mode  equations  [19].  Next,  we  match  boundary 
conditions  at  the  facets  and  obtain 


-A(itotal) 

.  0)  . 


=  D  1F(Ltotal  -  zs) 


821 

S22 


(13) 


where  the  matrix  D  contains  the  boundary  condition  informa¬ 
tion  at  the  facets 


1  -riFii(Ltotal)  “  ^12(£total) 

?2  —T1F21  (itotal)  —  ^(Aotal) 


(14) 


and  and  7*2  are  the  complex  reflectivities  at  the  facets.  The 
matrix  F(Ltotai )  is  the  product  of  all  the  transfer  matrices  in 
(10),  ancHt  describes  the  total  three-section  device  whereas  the 
matrix  i^Ltotai  —  Zs)  contains  only  those  matrices  between 
z3  and  T/totai* 

The  amplified  spontaneous  emission  intensity  from  each 
facet  is  obtained  by  summing  the  square  of  the  electric  field 
envelopes  J3(0)  and  A(L)  due  to  a  single  source  at  z8  over 
all  the  sources  J8  within  the  cavity  (0  <  z8  <  L)  as  given  by 


Jleft  = 


-bright  — 


(1  -  hi2)  fL\B(0)\2dzs 
Jo 

(l-\r2\2)£\A(L)\2dzs. 


(15) 

(16) 


The  ASE  is  calculated  using  (15)  and  (16)  for  each  wavelength 
to  get  the  entire  ASE  spectrum. 

The  photon  density  profile  is  obtained  from  the  electric  field 
envelopes  using 


S(z)  oc  \A(z)\2  +  \B(z)\2.  (17) 


Again,  the  total  photon  density  is  found  by  summing  the 
photon  density  S(z)  due  to  a  single  source  at  over  all  the 
sources  J8  within  the  cavity. 


B .  Spatial  Hole  Burning 

In  the  first  iteration  of  the  photon  density  calculation,  the 
gain  and  carrier  density  profiles  along  the  device  are  consid¬ 
ered  uniform.  If  the  applied  injection  current  is  above  thresh¬ 
old,  the  calculated  photon  density  profile  will  be  perturbed  as 
a  function  of  longitudinal  cavity  position  and  therefore  not 
constant.  The  carrier  density  and  gain  are  recalculated  using 
the  carrier  density  rate  equation  at  steady  state  given  by 

^  =  jy-AN-  BN 2  -  CN 3  -  vgg(N)S  (18) 

where  I  is  the  injected  current,  e  is  the  electron  unit  charge,  V 
is  the  volume  of  the  laser  active  region,  A  is  the  monocarrier 


lifetime  coefficient,  B  is  the  bimolecular  carrier  recombination 
coefficient,  C  is  the  Auger  recombination  coefficient,  vg  is 
the  group  velocity,  g(N )  is  the  material  gain  as  a  function 
of  carrier  density  N ,  and  S  is  the  photon  density.  In  the 
model,  we  consider  a  linearized  material  gain  given  by  g(N)  = 
g*( N  —  No),  where  gf  is  the  differential  gain  and  N0  is  the 
carrier  density  at  transparency.  This  results  in  a  nonuniform 
carrier  density  and  profile.  The  refractive  index  profile  also 
becomes  nonuniform  due  to  the  perturbed  carrier  density 
profile  through  the  relation 

A ne  =  aat  (19) 

47 r 

where  A ne  is  the  change  in  effective  index,  A  is  the  wave¬ 
length,  ae  is  the  linewidth  enhancement  factor,  and  gr  is 
the  differential  gain.  This  change  in  refractive  index  to  a 
nonuniform  carrier  density  profile  changes  the  threshold  gain 
of  potential  lasing  modes  and  accounts  for  the  decrease  in  the 
sidemode  suppression  ratio  and  mode  hopping  associated  with 
SHB.  The  photon  density  and  carrier  density  are  continually 
recalculated  using  (18)  until  a  self-consistent  solution  is  found. 
Only  a  single  longitudinal  mode  is  taken  into  account.  Our 
model  has  been  applied  to  the  study  of  SHB  in  conventional 
DFB  lasers  with  excellent  agreement  with  experimental  data 
[19].  It  should  be  noted  that  each  section  of  the  device  is 
further  subdivided  into  multiple  sections,  and  the  coupled¬ 
mode  equations  are  used  to  derive  the  transfer  matrices  to 
account  for  the  spatially  dependent  refractive  index  and  gain 
profiles  and  the  grating  pitch  parameter. 

C.  Curved  Waveguide  and  Dispersion 

The  curved  waveguide  in  the  DFB  laser  section  induces 
an  effective  variation  in  grating  pitch  Aeff  as  a  function  of 
longitudinal  position  and  is  included  by  [11] 


a  Ao 

cos  <(>(z) 

(20) 

«.)  =  tan’1 

z 

(21) 

where  y  and  2  are  defined  in  (1),  and  Ao  is  the  grating  pitch  for 
the  straight  waveguide  case.  Aeff  is  taken  into  account  in  the 
model  through  the  propagation  constant  fJc  =  7r/Aeff,  whereas 
in  our  previous  work  [19]  f)0  =  7r/A0.  In  addition,  the  use  of  a 
curved  waveguide  slightly  increases  the  effective  length  of  the 
waveguide  in  the  DFB  laser  section  through  which  the  fields 
propagate.  This  is  taken  into  account  in  the  model. 

The  effects  of  dispersion  are  also  considered  in  the  model. 
Over  a  narrow  range  of  A  near  a  reference  wavelength  Ao, 
the  effective  index  ne(A)  can  be  approximated  by  a  linear 
function  of  A 


ne(A)  =  ne(A0)  —  s(A  —  A0)  (22) 

where  the  slope  s  is 

s  =  (ng  ~ne(Ao))  (23) 

Ao 

and  Tig  is  the  group  index,  which  can  be  approximated  as  a 
constant  within  the  wavelength  range  of  interest.  We  estimate 
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the  group  refractive  index  from  the  mode  spacing  outside 
the  stopband  in  the  experimental  spectrum,  and  the  value  of 
s  is  estimated  using  previous  experimental  results  [21]  for 
the  InGaAsP  system.  We  then  calculate  ne(A0)  and  solve  for 
ne( A)  using  (22)  over  the  wavelengths  of  interest. 


D.  Voltage-Dependent  Refractive  Index  and  Gain 
in  the  Electoabsorption  Modulator 


Typically  in  EML’s,  the  bandgap  of  the  modulator  is  slightly 
larger  than  that  in  the  DFB  laser  section  with  a  detuning  of 
about  30-60  nm  [8],  [9],  so  that  for  zero  or  forward  modulator 
voltage  bias  there  is  no  attenuation  except  for  insertion  loss. 
For  a  positive  voltage  bias,  carriers  are  injected  so  that  light 
propagating  through  the  modulator  experiences  gain,  and  the 
refractive  index  of  the  modulator  decreases.  If  a  negative 
voltage  bias  is  applied  to  the  modulator,  absorption  is  increased 
due  to  the  quantum-confined  Stark  effects  or  Franz-Keldysh 
effects  which,  in  turn,  increase  the  refractive  index  of  the 
modulator.  We  model  the  voltage-dependent  gain  (negative 
absorption)  and  refractive  index  in  the  modulator  through  a 
simple  linear  additive  term  as  given  by 


9mod  —  Po  + 


^mod  —  Ti0 


dg 

dVm  od 
dn 

dVm  od 


^mod 

(24) 

^mod 

(25) 

where  gmod  and  nmod  are  the  gain  coefficient  (1/cm)  and 
effective  index  in  the  modulator,  respectively,  gQ  and  n0  are 
the  gain  and  effective  index  in  the  modulator  with  no  voltage 
bias,  dgldVmod  and  dnldV1  nod  are  empirical  constants,  and 
Vmod  is  the  voltage  bias  at  the  modulator.  Note  that  both 
dgidVm od  and  dn/dVmod  are  positive.  Because  the  change  in 
refractive  index  in  the  forward  bias  case  is  caused  by  carrier 
injection  while  the  change  in  refractive  index  in  the  reverse 
bias  case  is  caused  by  the  electroabsorption  effect,  the  value  of 
dnldVmod  is  different  for  each  case.  In  both  cases,  the  change 
in  refractive  index  in  the  modulator  due  to  the  voltage  bias 
is  the  mechanism  responsible  for  the  phase-shifting  effect,  as 
will  be  explained  irr  Section  IV. 


E.  Leakage  Current  from  the  Modulator 

The  function  of  the  fabricated  isolation  region  is  to  provide 
electrical  isolation  between  the  modulator  and  laser  sections  so 
that  the  biasing  of  one  section  does  not  affect  the  performance 
of  the  other.  However,  the  isolation  region  in  this  particular 
device  does  not  provide  perfect  electrical  isolation.  When 
the  modulator  is  forward-biased,  carriers  are  injected  into 
the  modulator  and  a  small  leakage  current  flows  from  the 
modulator  into  the  laser  section.  This  leakage  current  will 
decrease  the  refractive  index  of  the  laser  section  and  is 
included  in  the  model  using  a  simple  empirical  relation 


^DFB  —  no 


drijeak 
dV mod 


^mod 


(26) 


where  tidfb  is  the  modified  effective  index  of  the  laser 
section,  n0  is  the  original  effective  index  and  dnieak/dVmod 
is  the  change  in  effective  index  due  to  the  leakage  current 


and  is  determined  empirically  by  the  observed  shift  in  lasing 
wavelength  in  the  measured  spectrum  as  a  function  of  Vmod. 

IV.  Results  and  Discussions 

A.  Wavelength-Tuning  Mechanisms — Phase-Shifting 
and  Current  Tuning 

Although  our  transfer  matrix  method  provides  the  full 
numerical  solution  taking  into  account  all  multiple  reflection, 
we  present  a  simple  physical  model  to  illustrate  the  physics 
of  the  tuning  mechanisms  below.  Since  the  isolation  region 
and  modulator  are  homogeneous,  meaning  that  there  is  no 
distributed  feedback  present,  the  round  trip  through  these 
sections  can  also  be  perceived  as  an  effective  reflectivity  at  the 
DFB-isolation  region  interface.  Varying  the  modulator  voltage 
will  change  the  magnitude  and  phase  of  the  reflectivity  at  this 
interface;  hence,  the  modulator  acts  as  a  phase-shifting  section. 
The  voltage-dependent  magnitude  and  phase  of  the  interface 
reflectivity  is  approximately 

Tint  =  rH  Ret2f)i‘°  Liso  gi2/3mo<l  Lmod  (27) 

where  rint  and  run  are  the  magnitudes  of  the  field  reflectivity 
at  the  laser-isolation  region  and  HR  facet,  respectively,  L ^ 
and  Lmod  are  the  lengths  of  the  isolation  region  and  modulator, 
respectively,  /?is0  is  the  propagation  constant  of  the  isola¬ 
tion  region,  and  (3mo a  is  the  voltage-dependent  propagation 
constant  of  the  modulator  given  by 

n  27T  .gm  od 

Pmod  —  ^  Tlmod  *  2  \^®) 

where  gmod  and  nmod  are  functions  of  modulator  voltage  and 
are  given  by  (24)  and  (25),  respectively.  Therefore,  the  change 
in  phase  at  the  interface  due  to  a  voltage-dependent  change  in 
refractive  index  is  given  by 

A<p  =  47rAneLmoJ  (29) 

A 

where  A ne  is  given  by  the  voltage-dependent  term  in  (25). 

It  is  well  known  that  in  an  HR-/AR-coated  DFB  laser 
the  lasing  wavelength  can  change  dramatically  for  different 
phases  at  the  HR  facet  [22].  By  integrating  a  phase-tuning 
section  into  the  HR-coated  side  of  the  DFB  laser,  we  use  this 
idea  to  achieve  voltage-controlled  wavelength  tuning  over  the 
stopband  of  the  DFB  laser. 

The  other  mechanism  behind  wavelength  tuning  in  this 
integrated  device  is  due  to  leakage  current  from  the  modulator. 
An  increasing  positive  modulator  voltage  will  inject  carriers 
into  the  laser  section,  decrease  the  refractive  index  of  the  laser, 
and  shift  the  spectrum  toward  longer  wavelengths.  For  larger 
negative  voltage  biases,  carriers  are  absorbed  from  the  laser 
and  the  spectrum  shifts  toward  shorter  wavelengths.  Although 
leakage  current  is  not  usually  viewed  favorably,  this  effect, 
when  combined  with  the  phase-shifting  effect,  actually  allows 
this  device  to  achieve  a  larger  wavelength-tuning  range  than 
phase  shifting  alone. 

The  wavelength-tuning  behavior  of  this  device  naturally 
leads  to  a  discussion  on  the  advantages  and  disadvantages 
of  current  tuning  and  phase  shifting  as  wavelength-tuning 
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mechanisms  in  lasers.  The  idea  of  using  different  injection 
currents  and  operating  temperatures  to  produce  a  change  in  the 
wavelength  of  the  lasing  mode  has  long  been  known  and  has 
been  extended  to  a  variety  of  other  types  of  devices,  including 
multisection  lasers.  The  advantage  of  current  tuning  is  that  one 
can  achieve  a  continuous  tuning  range.  The  main  disadvantage 
is  that,  for  wavelength  switching,  as  used  in  traffic  routing 
applications,  current  tuning  is  a  relatively  slow  process,  lim¬ 
ited  by  the  carrier  lifetime,  and  cannot  support  wavelength 
switching  at  rates  in  the  multigigahertz  range.  In  order  to 
increase  wavelength  switching  rates,  several  groups  have  taken 
advantage  of  the  faster  Franz-Keldysh  effects  in  a  reverse- 
biased  electroabsorption  modulator  to  produce  wavelength 
switching.  Nakajima  et  al  demonstrated  wavelength  switching 
at  10  Gb/s  over  a  0.2-nm  range  [10].  This  wavelength  tuning 
was  achieved  by  absorption  of  carriers  by  the  modulator  and 
is  essentially  voltage  tuning  using  a  faster  process. 

Using  the  same  type  of  device  presented  here,  it  has  been 
shown  that  wavelength  switching  at  2,5  Gb/s  between  four 
wavelength  channels  spanning  a  2.1-nm  range  can  be  achieved 
[11].  The  wider  wavelength  range  achieved  here  is  due  to  the 
fact  that  both  phase  shifting  and  voltage  tuning  from  quantum- 
confined  Stark  effects  were  used  to  produce  the  wavelength 
tuning. 

B.  Spectrum  of  the  Integrated  DFB  Laser-Modulator 

Our  model  was  first  used  to  explain  some  of  the  features 
present  in  the  output  spectrum  of  this  curved  waveguide 
DFB  laser  with  an  integrated  modulator.  The  spectrum  of  a 
conventional  HR/AR-coated  DFB  laser  was  measured  above 
threshold  as  shown  in  Fig.  2(a).  The  DFB  laser  spectrum 
typically  has  a  lasing  mode  with  a  sidemode  suppression  ratio 
(SMSR)  greater  than  30  dB  on  one  side  of  the  stopband 
and  side  modes  with  amplitudes  which  decrease  gradually 
for  shorter  and  longer  wavelengths  away  from  the  stopband. 
The  measured  spectrum  of  the  integrated  device  with  no 
modulator  voltage  bias  is  shown  in  Fig.  2(b)  and  shows  several 
interesting  features  which  differ  from  the  conventional  DFB 
laser.  The  first  feature  is  the  absence  of  a  visible  stopband. 
The  second  feature  is  that  the  modes  alternate  in  amplitude, 
whereas  the  mode  amplitudes  of  the  conventional  DFB  laser 
gradually  decrease  away  from  the  stopband.  The  third  feature 
is  that  competing  side  modes  can  be  seen  on  both  sides  of  the 
lasing  mode,  whereas  in  the  conventional  DFB  laser,  there  is 
usually  only  one  competing  mode  located  on  the  other  side 
of  the  stopband.  The  theoretically  calculated  spectrum  shown 
in  Fig.  2(c)  exhibits  features  which  are  quite  similar  to  the 
measured  spectrum  shown  in  Fig.  2(b).  The  parameters  used 
in  simulations  of  the  integrated  wavelength-tunable  device  are 
shown  in  Table  I. 

The  spectrum  of  the  integrated  device  can  be  explained  by 
looking  at  the  threshold  gain  of  the  modes  as  a  function  of 
wavelength  for  several  cases.  First  we  consider  the  effect  of 
using  a  two-section  DFB  laser-modulator  device  compared  to 
the  single-section  DFB  laser  case.  Shown  in  Fig.  3(a)  are  the 
modes  for  a  single-section  DFB  laser  with  a  straight  waveguide 
case.  Shown  in  Fig.  3(b)  are  the  modes  for  the  two-section 
integrated  laser-modulator.  For  the  single-section  DFB  laser. 
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Fig.  2.  (a)  Measured  spectrum  of  a  conventional  HR/AR-coated  DFB  laser. 

The  spectrum  typically  has  a  lasing  mode  with  good  SMSR  on  one  side  of  the 
stopband  and  side  modes  with  amplitudes  which  decrease  gradually  for  shorter 
and  longer  wavelengths  away  from  the  stopband,  (b)  Measured  spectrum 
for  the  integrated  wavelength-tunable  laser.  This  spectrum  differs  from  the 
conventional  DFB  laser  in  that  the  stopband  is  less  visible,  the  modes  alternate 
amplitudes  away  from  the  stopband,  and  competing  side  modes  are  present 
on  either  side  of  the  lasing  mode,  (c)  Theoretically  calculated  spectrum  of  the 
integrated  wavelength-tunable  laser.  This  spectrum  shows  similar  features  to 
the  experimentally  measured  spectrum  in  (b). 


the  modes  show  a  well-defined  stopband.  However,  for  the 
two-section  DFB  laser-modulator  device,  the  modes  seem  to 
separate  into  two  sets  near  the  stopband,  where  the  set  of 
modes  with  higher  threshold  gain  corresponds  to  the  coupled- 
cavity  effects  from  the  modulator  region,  and  the  modes  with 
lower  threshold  gain  correspond  to  the  DFB  laser  region.  This 
accounts  for  the  alternating  mode  amplitude  behavior  away 
from  the  stopband  toward  shorter  and  longer  wavelengths  as 
seen  in  the  spectrum. 

In  a  straight  waveguide,  the  grating  pitch  is  constant  as 
a  function  of  longitudinal  cavity  position  which  produces 
two  well-defined  Bragg  frequency  modes  located  opposite 
each  other  across  the  stopband.  In  the  curved  waveguide,  the 
grating  pitch  is  chirped  as  a  function  of  cavity  position,  which 
produces  a  mode-flattening  effect  around  the  stopband  which 
increases  as  a  function  of  S-bend  width  W.  This  effect  is 
illustrated  by  simulating  a  single-section  AR/AR-coated  DFB 
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TABLE  I 

Parameter  List  Used  for  Simulations 


Power  reflectivity  at  HR  facet  J?hr  0.99 

Power  reflectivity  at  AR  facet  Rar  0.01 

Phase  at  HR  facet  <£hr  30° 

Phase  at  AR  facet  <Par  0° 

Coupling  coefficient-length  product  kLdfb  3.0 

Length  of  laser  section  I/DFB(Aim)  360 

Length  of  isolator  section  Liso  (/*m)  80 

Length  of  modulator  section  Z/MOD(Mm)  260 

Grating  pitch  A0  (jt m)  0.2308 

S-bend  parameter  L\  Qtm)  47 

S-bend  parameter  Lq  (/tm)  133 

S-bend  parameter  W  Qtm)  10 

Group  index  n(J  3.6 

Linewidth  enhancement  factor  aK  4 

Optical  confinement  factor  T  0.03 

Transparent  carrier  density  N0  (cm-3)  1.5x  10“ 18 

Differential  gain  g*  (cm2)  2.7xl0-18 

dn^k/VJnod  (V-1)  0.002169 

dg/Vx mod  (cm-1  V~l)  8.0 

dn/V ;nod  for  ymot,  >  0  (V-1)  0.013135 

dn/Vmful  for  Vinod  <  0  (V-1 )  0.01 129 

Group  velocity  vu  (m/s)  8.33 XlO8 

Dispersion  parameter  .s  (1/^m)  0.21930 
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Fig.  3.  (a)  Modes  for  a  one-section  DFB  laser,  where  the  modes 
show  a  well-defined  stopband,  (b)  Modes  for  the  two-section  integrated 
laser- modulator:  the  modes  show  a  splitting  near  the  stopband,  the  upper 
modes  corresponding  to  the  modulator  region  and  the  lower  modes 
corresponding  to  the  DFB  laser  region. 

laser  in  order  to  isolate  the  effects  of  using  the  curved  wave¬ 
guide  and  to  more  clearly  illustrate  the  effect.  In  Fig.  4(a)-(c), 
the  threshold  gain  of  the  theoretically  calculated  modes  are 
plotted  versus  wavelength  for  increasing  values  of  W.  The 
mode-flattening  effect  comes  from  the  fact  that,  using  a  curved 
waveguide  or  chirped  grating  pitch,  the  Bragg  wavelength  is  no 
longer  defined  as  a  single  frequency  but  a  range  of  frequencies. 
Therefore,  the  distributed  feedback  is  weaker  for  one  particular 


Fig.  4.  The  theoretically  calculated  threshold  gain  of  the  modes  for  an 
AR/AR-coated  DFB  laser  as  a  function  of  wavelength  for  (a)  W  =  0  jtm, 
(b)  W  =  5  fi m,  and  (c)  W  =  15  f.im.  In  general,  increasing  the  width  of 
curved  waveguide  W  produces  a  mode-flattening  effect  around  the  stopband, 
which  can  enhance  mode  tunability. 

frequency,  and  the  SMSR  as  determined  by  the  threshold 
gain  decreases  as  a  result.  This  mode-flattening  effect  around 
the  stopband  accounts  for  the  presence  of  multiple  competing 
side  modes  in  the  spectrum  of  the  integrated  device  as  well 
as  the  absence  of  a  stopband.  Although  a  reduced  SMSR 
is  not  desirable,  this  mode-flattening  effect  may  extend  the 
wavelength-tuning  range  of  this  type  of  device  by  providing 
more  competing  modes  for  lasing,  which  may  be  selected  by  * 
the  modulator  section.  Far  outside  the  stopband,  the  mode 
behavior  was  observed  to  be  determined  by  the  Fabry-Perot 
cavity  in  both  experiment  and  theory  as  expected.  Simulations 
of  the  modes  for  an  HR/AR-coated  DFB  laser  show  a  similar 
mode-flattening  effect  due  to  the  curved  waveguide. 

C.  Wavelength  Tuning  of  the  Lasing  Mode 

The  model  was  used  to  simulate  the  wavelength-tuning 
behavior  in  the  integrated  DFB  laser-modulator  for  several 
values  of  modulator  voltage.  Fig.  5(a)  shows  the  measured 
spectrum  of  the  integrated  device  for  modulator  voltage  val¬ 
ues  of  0,  +0.47,  +0.78,  +0.87,  +0.92,  and  +1.03  V.  As 
explained  earlier,  the  wavelength  tuning  occurs  due  to  a 
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Fig.  5.  (a)  Measured  spectrum  of  the  integrated  device  at  modulator  voltage 

values  of  0,  +0.47,  +0.78,  +0.87,  +0.92,  and  +1.03.  (b)  Theoretically 
calculated  spectrum  of  the  integrated  devices  at  modulator  voltage  values 
of  0,  +0.65,  +0.75,  +0.85,  +0.95,  and  +1.10. 

combination  of  phase  shifting  and  current  tuning.  In  general, 
the  wavelength-tuning  behavior  of  this  device  can  be  described 
in  the  following  way.  As  the  modulator  voltage  is  increased 
from  0  V,  the  lasing  mode  remains  at  one  wavelength  until 
+0.47  V,  where  the  lasing  mode  switches  abruptly  to  a  second 
wavelength.  As  the  modulator  voltage  is  increased  further,  the 
lasing  mode  remains  at  this  second  wavelength  until  +0.78 
V,  where  the  lasing  mode  switches  to  another  wavelength, 
and  this  behavior  continues  until  +1.03  V.  Fig.  5(b)  shows 
the  theoretically  calculated  spectrum  for  modulator  voltages 
of  0,  +0.65,  +0.75,  +0.85,  +0.95,  and  +1.10  V,  which 
shows  reasonably  good  agreement  using  the  linear  model  for 
voltage-dependent  refractive  index. 

For  a  positive  voltage  bias,  the  modulator  gain  and  refractive 
index  could  be  modeled  more  accurately  by  making  them 
functions  of  carrier  density  in  the  modulator,  which  is  a 
function  of  injected  current  in  the  modulator.  Also,  wavelength 
tuning  of  the  lasing  mode  for  positive  bias  voltages  is  also 
associated  with  changes  in  threshold  current,  linewidth,  and 
output  power.  This  was  not  included  here,  but  we  refer  to  [24] 
for  a  more  detailed  treatment. 

Wavelength  tuning  of  the  lasing  mode  was  also  observed 
for  negative  modulator  bias  voltages,  although  for  a  more 
limited  range  for  this  particular  device.  While  in  the  positive 
voltage  bias  case  the  refractive  index  decreases  due  to  the 
injection  of  carriers,  for  this  case,  the  refractive  index  increases 
due  to  absorption  of  carriers  from  quantum-confined  Stark 
effects.  Fig.  6(a)  shows  the  measured  spectra  for  modulator 
voltage  values  of  0,  -0.45,  and  -0.78  V.  Fig.  6(b)  shows 
the  theoretically  calculated  spectra  for  modulator  values  of 
0,  -0.50,  and  -0.90  V,  which  show  good  agreement  with 
the  measured  spectra.  The  lasing  wavelength  as  a  function  of 
positive  and  negative  modulator  voltages  is  shown  for  both  the 


Fig.  6.  (a)  Measured  spectra  of  the  integrated  device  at  modulator  voltage 

values  of  0,  —0.45,  and  -0.78.  (b)  Theoretically  calculated  spectra  of  the 
integrated  devices  at  modulator  voltage  values  of  0,  —0.50,  and  —0.90. 


Fig.  7.  Measured  and  theoretically  calculated  values  of  the  lasing  wavelength 
as  a  function  of  modulator  voltage.  The  closed  circles  represent  measured 
lasing  wavelengths.  The  lines  represent  observed  lasing  wavelengths  and 
are  accurate  to  within  ±1  A.  The  open  triangles  represent  the  theoretically 
calculated  values. 


experimentally  measured  and  theoretically  calculated  spectra 
in  Fig.  7.  The  closed  circles  are  measured  data  points.  The 
lasing  wavelengths  between  data  points,  shown  by  the  lines, 
were  observed  but  not  measured  and  are  accurate  to  within  ±1 
A.  The  open  triangles  are  the  theoretically  calculated  values 
and  show  reasonably  good  agreement  with  the  measured 
spectrum  values. 


D.  Reduced  SHB  Due  to  the  Curved  Waveguide 

SHB  is  an  important  phenomena  which  has  been  studied 
extensively  [2,5] — [32],  particularly  in  DFB  lasers,  because  it 
can  severely  limit  the  maximum  output  power  and  cause 
mode  hopping  or  multimode  behavior.  SHB  is  caused  by 
a  nonuniform  distribution  of  the  photon  density  along  the 
longitudinal  cavity  which  can  be  caused  by  the  distributed 
feedback  in  DFB  lasers  as  well  as  asymmetries  in  HR/AR- 
coated  DFB  lasers.  A  greater  photon  density  in  one  region 
of  a  DFB  laser  will  cause  greater  carrier  depletion  in  that 
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Fig.  8.  The  theoretically  calculated  photon  density  as  a  function  of  longi¬ 
tudinal  cavity  position  for  the  straight  waveguide  (W  ~  0  fim)  and  curved 
waveguide  (W  =  10  /urn)  cases.  The  solid  lines  are  for  Vmod  =  0,  and  the 
dashed  lines  are  for  a  different  lasing  mode  at  Knod  =  +0.65.  The  curved 
waveguide  shows  less  perturbation  in  the  photon  density  profile  and  therefore 
less  SHB  effects  compared  to  the  straight  waveguide  case. 


region  so  the  carrier  density  profile  also  becomes  nonuniform 
along  the  cavity.  This  in  turn  produces  a  nonuniform  refractive 
index  profile  which  affects  the  threshold  gain  conditions  of  the 
modes  and  can  result  in  mode  hopping  or  multimode  lasing  at 
moderate  injection  currents  above  threshold.  The  model  was 
used  to  calculate  the  longitudinal  photon  density  profile  for  the 
straight  and  curved  waveguide  cases,  as  shown  in  Fig.  8.  The 
corresponding  facet  coatings  and  device  sections  are  labeled 
with  cavity  position.  The  solid  lines  are  for  V^noa  =  0  V,  and 
the  dashed  lines  are  for  Vmofx  =  -1-0.65  V.  Considering  only 
the  laser  section,  the  photon  density  profile  shows  a  larger 
perturbation  in  the  straight  waveguide  case  than  that  for  the 
curved  waveguide  case.  One  way  to  evaluate  the  amount  of 
SHB  present  in  a  DFB  laser  is  to  calculate  a  flatness  parameter 
[33].  We  use  a  modified  flatness  parameter  given  by 


1)  The  present  device  mode  shifts  between  only  two  modes. 
By  optimizing  the  curved  waveguide  shape  and  increas¬ 
ing  the  mode-flattening  effect,  we  should  be  able  to 
obtain  mode  shifting  between  three  or  more  modes, 
thereby  increasing  the  wavelength-tuning  range.  How¬ 
ever,  the  tradeoff  would  be  a  decrease  in  SMSR. 

2)  The  concept  of  an  amplitude-modulated  grating  [37]  can 
be  used  to  extend  the  phase-shifting  range.  Although 
originally  designed  to  flatten  the  photon  profile  and 
reduce  SHB  effects,  this  idea  may  also  be  used  to 
enhance  the  phase-shifting  range. 

V.  Conclusions 

We  have  theoretically  modeled  a  DFB  laser  with  a  curved 
waveguide  integrated  with  an  electroabsorption  modulator 
which  acts  as  a  wavelength-tuning  section.  The  characteristics 
of  the  output  spectrum  have  been  identified,  and  the  mecha¬ 
nisms  behind  the  wavelength  tuning,  namely,  phase  shifting 
and  voltage  tuning,  have  been  explained  using  simulation 
results  which  show  good  agreement  with  experiment.  The 
model  also  shows  that  the  curved  waveguide  suppresses  the 
photon  density  profile  and  reduces  SHB  compared  to  the 
straight  waveguide  case.  The  development  of  wavelength- 
tunable  lasers  with  greater  tuning  range  and  higher  power 
operation  should  be  achievable  using  a  curved  waveguide. 
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F=  ~  J  (/(*> n0rm  ~  4v g)2  dz  (30) 
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We  present  high  accuracy  measurements  of  gain,  loss,  and  transparency  energy  in  long-wavelength 
semiconductors  based  on  a  hybrid  approach  using  the  fundamental  relationship  between  the  gain 
and  the  spontaneous  emission  spectra.  Independent  measurements  of  optical  gain,  transparency 
energy,  and  loss  show  the  accuracy  and  validity  of  this  technique.  These  results  are  compared  with 
those  obtained  by  the  non-Markovian  gain  model  with  many-body  effects  under  the  spontaneous 
emission  transformation  method.  It  is  found  that  the  hybrid  approach  for  the  gain  spectrum 
alleviates  many  of  the  problems  related  to  the  poor  signal  to  noise  ratio  in  the  amplified-spontaneous 
emission  near  and  below  the  band  edge.  The  theoretical  spectra  compare  well  with  the  measured 
spectra  for  both  the  transverse  electric  and  transverse  magnetic  polarizations.  ©  1999  American 
Institute  of  Physics.  [S002 1-8979(99)010 18-X] 


I.  INTRODUCTION 

Optical  gain  and  loss  are  of  fundamental  importance  to 
the  study  of  lasers,  since  the  gain  determines  many  device 
performance  characteristics.  Therefore,  accurate  and  reliable 
gain  measurement  techniques  are  essential  tools  for  the  ex¬ 
perimental  study  of  lasers.  Likewise,  highly  accurate  theoret¬ 
ical  gain  models  incorporating  the  most  realistic  physical 
effects  are  also  essential  to  advance  the  state  of  the  art.  Such 
models  are  most  stringently  tested  against  the  best  obtainable 
data,  since  the  usefulness  of  theoretical  predictions  of  gain 
and  carrier  density  depends  both  on  the  model  and  on  close 
agreement  with  accurately  measured  data. 

On  the  experimental  side,  several  excellent  and  widely 
used  techniques  already  exist  for  measurements  of  gain  and 
loss,  yet  there  are  instances  where  improvements  are  still 
needed.  Measuring  the  gain  in  a  Fabry-Perot  or  distributed- 
feedback  laser  typically  first  involves  measuring  the  ampli¬ 
fied  spontaneous  emission  (ASE)  of  the  laser.  Since  the  ASE 
signals  are  very  weak  below  the  band  edge,  however,  the 
evaluation  of  gain  may  degrade  strongly  due  to  poor  signal 
to  noise  ratios  in  the  ASE.  As  a  consequence,  commonly 
used  techniques  for  measuring  the  modal  gain,  such  as  the 
Hakki-Paoli  method1  or  Cassidy’s2  method,  show  a  signifi¬ 
cant  amount  of  noise  in  the  region  of  the  gain  tail  because  of 
reduced  signal  intensity.  While  some  reduction  in  the  noise 
susceptibility  of  both  methods  can  be  achieved  through  an 
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averaging  procedure,3  the  results  thus  obtained  are  less  accu¬ 
rate. 

Difficulties  in  measuring  the  gain  tail  can  compound  to 
create  further  problems,  since  estimates  of  the  intrinsic  loss 
are  often  obtained  by  examining  the  limit  of  the  gain  tail.  For 
instance,  assuming  the  intrinsic  losses  are  constant  with  re¬ 
spect  to  wavelength,  the  measured  loss  may  be  used  to  find 
the  transparency  energy.  But  it  is  impossible  to  do  so  if  the 
losses  are  inaccurate.  In  turn  each  of  these  problems  impacts 
the  comparison  between  theory  and  experiment.  Accurate 
values  of  the  loss  are  needed  to  properly  renormalize  theo¬ 
retical  spectra.  Moreover,  an  incorrect  renormalization  will, 
in  effect,  shift  the  location  of  the  transparency  energy.  Thus, 
clear  data  in  the  tail  region  is  necessary  for  comparison  to 
theoretical  gain  models  in  terms  of  predicted  line  broadening 
behavior. 

On  the  theoretical  side,  conventional  optical  gain  calcu¬ 
lations  are  usually  based  on  the  density-matrix  theory  with  a 
phenomenological  damping  term  which  gives  the  Lorentzian 
line  shape  function.4  A  few  groups5”8  have  argued  that  the 
optical  gain  spectra  calculated  with  the  Lorentzian  line  shape 
function  deviate  from  the  experimental  results.  Especially,  an 
anomalous  absorption  region  below  the  band  gap  and  dis¬ 
crepancies  between  the  transparency  points  and  the  Fermi- 
level  separation  appear  in  the  gain  spectra  when  the  Lorent¬ 
zian  line  shape  is  used.  To  overcome  these  artifacts,  the 
intraband  relaxation  times  by  the  carrier-carrier  and  carrier- 
longitudinal  optical  (LO)  phonon  scatterings  have  been  con¬ 
sidered  to  obtain  the  line  shape  function.7,8  Its  low-energy 
tail  decreases  much  faster  than  the  Lorentzian  line  shape 
while  its  high-energy  tail  is  close  to  the  Lorentzian  line 
shape.  This  shows  that  the  steepness  of  the  low-energy  tail  of 
the  line  shape  function  is  most  important.  However,  these 
models  require  a  large  amount  of  computational  time.  As  a 
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simplified  approach,  approximate  formula  considering  only 
the  hole-hole  scattering  has  often  been  used  in  calculating 
the  optical  gain  spectrum.9-11  Another  approach  is  to  replace 
the  Lorentzian  line  shape  function  with  a  sharper  spectral 
function,  i.e.,  with  a  Gaussian  line  shape.  The  effect  of 
Gaussian  line  shape  is  similar  to  that  of  the  line  shape  func¬ 
tion  considering  the  carrier-carrier  and  carrier-LO  phonon 
scattering. 

On  the  other  hand,  the  spontaneous  emission  transforma¬ 
tion  method  using  the  fundamental  relations  between  the 
spontaneous  emission  and  the  absorption  coefficients  has 
been  proposed  as  an  improved  method.5’12-15  It  guarantees 
that  artifacts  mentioned  above  are  circumvented  independent 
of  the  line  shape  functions  (Lorentzian,  Gaussian,  or  others). 
Also,  the  slow  convergence  of  the  Lorentzian  leads  to  a  very 
long  tail  of  the  gain  spectra  into  the  band  gap,  in  contradic¬ 
tion  to  experimental  observations.15  Hence,  a  non-Markovian 
(Gaussian)  gain  model  under  the  spontaneous  emission  trans¬ 
formation  method  is  desirable  to  obtain  the  optical  gain 
spectrum.16,17 

In  this  article,  we  present  high  accuracy  measurements 
of  the  optical  gain,  using  the  fundamental  relationship  be¬ 
tween  the  spontaneous  emission  spectrum  and  the  gain  spec¬ 
trum.  The  values  for  transparency  energy  and  loss  obtained 
from  this  technique  compare  favorably  with  those  obtained 
from  several  independent  measurements.18  Experimentally, 
the  major  techniques  involve  measuring  (1)  gain  and  sponta¬ 
neous  emission,  (2)  polarization  resolved  gain  spectra,  (3) 
optical-pump  induced  modified  gain  spectra,  and  (4)  pump- 
induced  electrical  responses19  for  directly  obtaining  the 
transparency  energy  and  losses.  The  key  information  each 
technique  provides  is  summarized  as  follows.  First,  since 
spontaneous  emission  and  gain  are  related  by  only  three  pa¬ 
rameters,  the  transparency  energy,  the  loss,  and  a  coupling 
constant,  matching  gain  and  spontaneous  emission  provides 
a  way  to  estimate  these  parameters.  Second,  the  quasi-Fermi 
level  separation  is  the  same  at  a  fixed  injection  current  for 
each  polarization,  so  the  polarization-resolved  gain  spectra 
will  cross  at  the  same  quasi-Fermi  level  separation.  Third, 
the  net  absorption  of  an  optical  pump  will  be  zero  at  the 
quasi-Fermi  level  separation,  so  the  change  in  gain  will  be 
zero  for  a  pump  wavelength  located  at  that  energy.  Finally, 
the  change  in  induced  device  voltage  due  to  absorption  of  an 
external  optical  pump  light  will  be  zero  when  the  pump 
wavelength  is  at  the  transparency  energy. 

The  transformed  spontaneous  emission  procedure  yields 
highly  accurate  measurements  of  the  gain  tail,  which  are  dif¬ 
ficult  to  obtain  from  the  other  techniques.  The  experimen¬ 
tally  obtained  spectra  are  presented  along  with  theoretical  fits 
based  on  a  many-body  gain  model.  Theoretical  gain  models 
provide  estimates  of  the  carrier  density  within  the  active  re¬ 
gion  by  matching  the  predicted  gain  spectrum  to  the  mea¬ 
sured  gain.  For  the  best  accuracy  it  is  necessary  to  include  all 
of  the  most  important  physical  effects  which  determine  the 
gain.  A  non-Markovian  optical  gain  model  with  many-body 
effects  is  used  to  calculate  the  spectra.  The  plasma  screening, 
band-gap  renormalization  (BGR),  and  the  excitonic  or  the 
Coulomb  enhancement  (CE)  of  the  interband  transition  prob- 
abilitv  are  included  in  the  model.  Good  agreement  between 


theory  and  experiment  for  both  the  transverse  electric  (TE) 
and  transverse  magnetic  (TM)  polarizations  is  obtained.  Our 
model  also  provides  the  numerical  values  for  the  carrier  den¬ 
sities. 

Experimental  results  for  the  measurement  of  gain,  loss, 
and  transparency  energy  using  the  spontaneous  emission 
spectrum  and  its  fundamental  relationship  with  the  gain  is 
presented  in  Sec.  II  along  with  a  comparison  to  other  meth¬ 
ods  for  obtaining  these  quantities.  A  description  of  the  many- 
body  gain  model  and  comparisons  between  theoretical  and 
experimental  results  are  discussed  in  Secs.  IE  and  IV.  Fi¬ 
nally,  the  summary  and  conclusions  are  presented  in  Sec.  V. 


II.  EXPERIMENTAL  TECHNIQUES 

The  following  experiments  were  performed  using  a 
nominally  lattice-matched  AlGalnAs  multiple  quantum-well 
laser.  An  unstrained  quantum-well  laser  is  chosen  so  that 
there  will  be  a  significant  amount  of  gain  for  both  the  TE  and 
TM  polarizations.  The  laser  is  a  Fabry-Perot  device  with  a 
cavity  length  of  638  jim  and  one  facet  coated  for  high  (85%) 
reflectivity.  Round  trip  mirror  loss  is  11.4  cm”1.  The  device 
was  temperature  controlled  and  all  data  were  taken  at  25  °C. 
For  all  spectra  the  device  current  was  measured  to  within  10 
/jA  accuracy.  Threshold  at  this  temperature  was  11.7  mA. 
Further  structural  details  are  described  in  Sec.  IV. 

A.  Method  I:  Measurement  of  spontaneous  emission 
and  the  use  of  the  fundamental  relation 

As  described  above,  direct  measurements  of  the  gain 
provided  insufficient  data  in  the  low  energy  tail  region.  On 
the  other  hand,  very  accurately  measured  values  of  transpar¬ 
ency,  energy  and  loss  can  be  yielded  by  using  an  indirect 
technique  based  on  the  relationship  between  spontaneous 
emission  and  gain.  The  optical  gain  spectra  are  related  to  the 
spontaneous  emission  spectra  from  the  detailed  balance  be¬ 
tween  absorption  and  emission  of  photons.5,12-15  That  is,  the 
modal  gain  from  spontaneous  emission  is  given  by 

h3c2  I  E~  AE\ 

s(£)=8^1^exriiir)r"(E)  (1) 

in  which  rsp  is  the  spontaneous  emission  rate,  h-hflTr,  h  is 
the  Planck’s  constant,  T  is  the  temperature,  E=fto>,  T  is  the 
optical  confinement  factor,  and  A  F  is  the  transparency  en¬ 
ergy,  which  is  given  by  the  energy  separation  of  the  quasi- 
Fermi  levels  of  the  electrons  and  holes.  In  practice,  the  spon¬ 
taneous  emission  power  Psp(\)  is  measured  within  small 
resolution  intervals  d\  determined  by  the  spectrometer  over 
the  entire  spectrum, 

Psb(V=KiYrsp{x)dKV  (2) 

so  the  relation  in  Eq.  (1)  must  be  modified.  First,  the  energy 
is  expressed  in  terms  of  the  wavelength,  then  the  spontane¬ 
ous  emission  rate  is  substituted  by  the  power  through  Eq.  (2). 
The  result  is 
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Next  the  material  gain  is  written  in  terms  of  the  mea¬ 
sured  net  modal  gain  as 

Sm(M  =  rg(X)-a„  (4) 

where  T  is  the  optical  confinement  factor  and  cl{  is  the  in¬ 
trinsic  loss  for  the  optical  waveguide  mode.  Now  the  modal 
gain  can  be  written  in  terms  of  the  measured  spontaneous 
emission  power: 
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In  the  abovee  expression,  there  are  three  parameters  which 
relate  gain  to  spontaneous  emission:  the  transparency  energy 
A F,'  the  loss  a, ,  and  the  constant  Q ,  which  includes  the 
optical  coupling  constant.  The  losses,  which  actually  vary 
slightly  as  a  function  of  wavelength,  are  assumed  constant  in 
this  technique.  That  is,  the  determination  of  the  transparency 
energy  is  equivalent  to  determining  the  losses,  since  the  ma¬ 
terial  gain  at  that  energy  is  zero.  The  coupling  constant  and 
transparency  energy  are  therefore  used  as  fitting  parameters 
for  each  pair  of  gain  and  spontaneous  emission  spectra. 

In  principle  it  should  be  possible  to  fix  the  value  of  the 
coupling  constant  when  matching  a  pair  of  gain  and  sponta¬ 
neous  emission  spectra,  and  use  that  fixed  value  for  subse¬ 
quent  matches.  Alternatively,  the  best-fit  matches  of  the  gain 
and  spontaneous  emission  spectra  ought  naturally  to  result  in 
a  uniform  value  for  the  coupling  constant.  Recent 
studies20-22  employing  this  technique,  however,  have  permit¬ 
ted  some  variation  in  the  coupling  constant  from  spectrum  to 
spectrum.  But  allowing  the  coupling  constant  to  fluctuate, 
even  by  a  few  percent,  can  lead  to  significant  changes  in  the 
extracted  values  of  transparency  energy  and  loss.  Therefore 
it  is  necessary  to  consider  whether  it  is  valid  to  allow  the 
coupling  constant  to  vary  for  each  pair  of  gain  and  sponta¬ 
neous  emission  spectra. 

In  the  experiments  the  spontaneous  emission  spectra  are 
measured  from  the  side  of  the  laser  using  a  bare  fiber.  It  is 
possible  to  mount  the  fiber  rigidly  and  isolate  it  from  vibra¬ 
tion  and  air  currents  so  that  the  position  of  the  fiber  is  held 
constant.  The  test  laser  is  temperature  controlled  and  biased 
at  a  constant  current,  and  monitored  throughout  the  duration 
of  the  experiments.  When  the  total  spontaneously  emitted 
power  is  measured  and  observed  for  stability,  it  varies  by 
less  than  one  percent  over  times  comparable  to  the  duration 
of  the  experiment.  Thus,  it  does  not  seem  reasonable  to  allow 
the  coupling  coefficient  to  vary  in  the  fitting  procedure. 


FIG.  1.  The  spontaneous  emission  spectra  at  four  injection  currents,  1-6 ,  8, 
10,  and  11.7  mA,  in  increasing  order.  The  wavelengths  corresponding  to  the 
gain  spectra  extending  from  1500  to  1650  nm. 


It  is  also  worthwhile  to  consider  wavelength  dependence 
of  the  coupling  coefficient.  A  multimode  silica  core  fiber  was 
used  to  couple  the  spontaneous  emission.  Figure  1  shows  the 
spontaneous  emission  for  four  injection  currents  up  to  laser 
threshold.  Here  the  portion  of  the  spontaneous  emission 
spectrum  which  is  transferred  to  obtain  gain  spectra  is  at 
wavelengths  of  1500  up  to  1650  nm.  In  coupling  to  the  fiber, 
small  variations  in  the  glass  refractive  index  (from  1.44462 
to  1.44267)  lead  to  variation  in  the  reflectivity,  and  a  simple 
estimate  can  be  made  using  the  Fresnel  formula,  i.e.. 


(«,+  !)* 


in  which  ne  is  the  effective  refractive  index  of  the  fundamen¬ 
tal  mode  of  the  fiber.  The  reflectivity  changes  from  3.31%  to 
3.28%  over  the  wavelength  range  of  1500-1650  nm  used  to 
match  gain  and  spontaneous  emission.  Absorption  within  the 
fiber  is  negligible.  So,  a  0.1%  variation  in  the  coupling  may 
be  reasonable,  though  in  obtaining  these  results  the  coupling 
constant  was  not  varied  with  respect  to  wavelength. 

Polarization  dependent  gain  data  are  extracted  by  mea¬ 
suring  the  ASE  spectra  through  a  polarizing  laser-optical  fi¬ 
ber  interface.  The  gain  is  then  calculated  using  the  Hakki- 
Paoli  method,  in  which  the  ratio  of  intensity  maxima  and 
minima  of  the  Fabry-Perot  modes  is  related  to  the  gain.1  As 
mentioned  earlier,  this  technique  relies  on  a  good  signal  to 
noise  ratio,  which  may  be  insufficient  for  wavelengths  near 
the  band  edge  where  the  ASE  intensity  is  extremely  low.  In 
the  peak  gain  region  for  currents  near  threshold,  the  modes 
are  almost  lasing  and  may  not  be  completely  resolved  by  the 
spectrometer.  The  mode  maxima  become  artificially  lowered 
and  the  Hakki-Paoli  method  underestimates  the  gain.  In  this 
region,  Cassidy’s  method,  employing  the  ratio  of  the  average 
mode  power  to  the  mode  minimum,  is  less  susceptible  to 
these  effects.  For  our  gain  measurements,  the  Hakki-Paoli 
technique  is  used  at  nearly  all  wavelengths,  except  near  the 
peak  gain  at  just  below  threshold,  where  Cassidy’s  method  is 
employed. 

To  obtain  the  best  fit  results  from  Eq.  (6),  we  choose  a 
constant  value  of  Q  to  be  used  for  each  of  the  currents.  The 
quasi-Fermi  level  of  separation  AF  is  then  varied  for  each 
current  to  yield  transformed  spontaneous  emission  spectra 
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FIG.  2.  Measured  net  modal  gain  spectra  (solid-TE  polarization,  dotted-TM 
polarization)  and  transformed  spontaneous  emission  spectra  (dashed-TE  po¬ 
larization)  are  shown  for  injection  current  1=6,  8,  10,  and  1 1.7  mA.  Note  the 
tail  regions  of  the  transformed  spectra  (on  the  long-wavelength  side),  which 
show  closer  agreement  for  TE  and  TM  polarizations. 

with  the  smallest  deviation  from  the  measured  modal  gain. 
The  value  of  Q  which  is  selected  is  that  which  gives  the  best 
simultaneous  fit  for  all  the  spectra.  The  resulting  transformed 
gain  spectra  then  provide  accurate  gain  data  in  the  tail  region 
where  none  was  available  from  the  direct  measurement  using 
the  Hakki-Paoli  or  Cassidy’s  method. 

The  measured  net  modal  gain  spectra,  Tg-a,, 
(solid-TE  polarization,  dotted-TM  polarization)  are  shown  in 
Fig.  2  with  the  transformed  spontaneous  emission  based  on 
the  above  fitting  techniques.  Overall  the  agreement  is  quite 
good;  but  the  improvement  is  most  noticeable  in  the  tail 
region,  where  the  transformed  spectra  (dashed)  show  signifi¬ 
cantly  less  noise  than  the  TE  spectra  (solid).  It  is  interesting 
to  note  the  tail  of  the  TE  gain  spectra,  which  initially  appear 
to  taper  off  at  a  loss  value  different  than  that  of  the  TM 
spectra.  The  transformed  data  show  that  the  actual  tail 
reaches  a  loss  point  in  much  closer  agreement  with  the  mea¬ 
sured  TM  losses.  The  extracted  values  of  transparency  en¬ 
ergy  and  loss  are  shown  in  Figs.  3(a)  and  3(b),  respectively. 
Also  shown  for  comparison  are  the  values  obtained  from  the 
independent  measurements  of  the  loss  and  transparency 
energies18  described  below.  The  extracted  values  for  all 
methods  are  in  close  agreement.  The  SE-transform  method  is 
no  better  or  worse  than  the  other  techniques.  The  advantage 
of  this  method,  however,  is  that  the  entire  gain  spectrum  can 
be  extracted  from  the  spontaneous  emission  spectrum,  yield¬ 
ing  very  accurately  measured  values  of  transparency  energy 
and  loss,  and  giving  a  more  accurate  tail  for  the  gain  spec¬ 
trum. 

B.  Method  II:  Crossing  of  polarization  dependent 
optical  gain  spectra 

Figure  4  shows  the  polarization  dependent  gain  spectra 
obtained  from  the  facet  ASE  spectrum.  The  transparency  en¬ 
ergy  may  be  measured  directly  by  examining  the  crossing 
point  of  the  TE  and  TM  gain  spectra,  and  the  value  of  loss  at 
the  crossing  point  is  equal  to  the  intrinsic  loss.  The  values 
extracted  from  comparison  of  TE  and  TM  gain  spectra  are 


Current  (mA) 

FIG.  3.  A  comparison  of  (a)  the  measured  transparency  energy  based  on  the 
TE-TM  crossing  method  (circles),  the  electrical  method  (triangles),  and  the 
spontaneous  emission  (SE)  transformation  method  (crosses)  and  (b)  the  in¬ 
trinsic  loss  using  the  TE-TM  crossing  method  (circles)  and  the  SE  transfor¬ 
mation  method  (crosses). 

presented  in  Fig.  3.  There  are  certain  flaws  in  measuring  the 
loss  and  transparency  in  this  manner.  Because  of  noise  near 
the  gain  tail  (especially  for  the  TE  spectra),  it  appears  that 
the  losses  may  be  polarization  dependent.  If  that  is  the  case, 
estimates  of  both  loss  and  transparency  energy  based  on  the 
crossing  of  the  TE  and  TM  gain  spectra  are  questionable.  Yet 
even  if  the  losses  were  not  polarization  dependent,  measur¬ 
ing  them  by  examining  the  crossing  point  of  TE  and  TM 
spectra  is  not  always  straightforward.  For  instance,  the  slope 
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FIG.  4.  Polarization  dependent  gain  spectra  for  currents  of  6, 8,  10,  and  1 1 .7 
mA,  from  bottom  to  top  curves,  are  shown.  The  plateau  at  the  long  wave¬ 
length  sides  gives  the  intrinsic  loss.  The  crossing  point  of  the  TE-  and 
TM-polarization  gain  spectra  gives  the  transparent  wavelength  at  each  in¬ 
jection  current.  Note  the  apparent  discrepancy  in  the  location  of  the  gain 
tails,  which  suggests  a  polarization-dependent  intrinsic  loss. 
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FIG.  5.  Gain  spectra  with  (dotted  and  dashed)  and  without  (solid)  an  exter¬ 
nal  optical  pump.  Pump  light  is  absorbed  at  1529  nm  and  amplified  at  1541 
nm,  causing  the  gain  spectrum  to  be  enhanced  (dotted)  or  compressed 
(dashed),  respectively. 


of  the  gain  is  relatively  steep  near  the  crossing  point,  so  that 
while  the  measured  transparency  energy  may  be  relatively 
accurate,  the  exact  value  of  loss  is  much  less  certain.  A  more 
general  problem  is  that  it  is  not  always  possible  to  obtain 
clear  TM  spectra,  especially  in  compressively  strained  lasers. 
Nevertheless,  measuring  the  TE-TM  crossing  is  straightfor¬ 
ward,  and  has  the  advantage  of  simplicity. 


C.  Method  III:  Measurement  of  gain  compression  in 
the  presence  of  an  external  laser  pump 

One  technique  for  measuring  the  transparency  energy 
involves  measuring  the  gain  of  a  semiconductor  laser  with  a 
constant  biased  current  in  the  presence  of  an  external  optical 
pump.  In  contrast  to  measurements  of  the  single-pass  gain  at 
the  pump  wavelength,  changes  in  the  overall  gain  spectrum 
induced  by  the  absorption  or  amplification  of  a  pump  are 
examined.  If  the  pump  wavelength  is  located  within  the  gain 
region,  amplification  of  the  pump  will  consume  carriers  and 
reduce  the  overall  gain.  Conversely,  if  the  pump  wavelength 
is  in  the  absorption  region,  the  pump  will  create  carriers  and 
enhance  the  gain.  In  these  experiments  an  external  distrib¬ 
uted  feedback  (DFB)  laser  is  used  as  the  pump  source,  and  it 
is  filtered  to  prevent  the  ASE  of  the  pump  DFB  laser  from 
entering  the  Fabry-Perot  laser  under  study.  Figure  5  shows 
measured  gain  spectra  at  6  mA  with  and  without  an  external 
pump.  First  the  pump  is  located  in  the  absorption  region  at 
1541  nm,  and  the  gain  is  enhanced.  Then  the  pump  is  tuned 
to  the  gain  region  at  1529  nm,  and  the  gain  is  compressed. 

The  two  pump  wavelengths  in  Fig.  5  are  separated  by  12 
nm  or  6.3  meV,  which  corresponds  to  an  uncertainty  in  the 
transparency  energy  of  less  than  1  %.  This  degree  of  uncer¬ 
tainty  in  the  transparency  energy  is  too  large  to  make  an 
accurate  comparison.  Fortunately,  one  can  interpolate  to  im¬ 
prove  the  estimate,  which  gives  an  estimated  transparency 
wavelength  of  1535  nm,  which  is  in  close  agreement  with  the 
results  obtained  by  TE-TM  crossing  method.  Likewise,  the 
estimated  losses  can  be  interpolated  from  the  pump  loca¬ 
tions,  for  an  intrinsic  loss  of  25  cm-1. 


FIG.  6.  A  plot  of  the  induced  ac  voltage  magnitude  across  the  laser  termi¬ 
nals,  due  to  the  injection  of  the  modulated  external  pump  laser,  as  the  dc 
bias  current  ( I0 )  of  the  test  laser  is  varied.  When  the  magnitude  of  the 
induced  ac  voltage  is  zero,  the  laser  wavelength  is  at  the  transparency  en¬ 
ergy. 


While  improvement  might  seem  possible  by  pumping  at 
wavelengths  closer  to  transparency,  in  practice  very  high 
pump  powers  are  required  since  the  absorption  or  gain  be¬ 
comes  very  small.  Little  gain  change  was  observed  for  pump 
wavelengths  between  1529  and  1541  nm  in  experiments  con¬ 
ducted  at  a  bias  of  6  mA,  though  pump  powers  as  high  as  3 
mW  were  used.  Other,  more  subtle  problems  arise  as  well; 
for  instance,  the  sensitivity  of  this  technique  varies  depend¬ 
ing  on  whether  the  pump  wavelength  overlaps  a  Fabry-Perot 
resonance  or  not.  If  the  pump  does  not  overlap  a  Fabry-Perot 
maximum,  the  gain  or  absorption  is  weaker  and  the  gain 
change  may  be  misleadingly  small.  A  slightly  different  ap¬ 
proach  based  on  an  electrical  characterization  technique,  de¬ 
scribed  in  the  next  section,  overcomes  these  difficulties  and 
provides  superior  accuracy  in  measuring  the  transparency  en¬ 
ergy  and  loss. 


D.  Method  IV:  Measurement  of  electrical 
characteristics  in  the  presence  of  an  external  laser 
pump 

An  excellent  method19  for  determining  the  transparency 
energy  to  a  high  degree  of  accuracy  also  involves  the  use  of 
an  external  pump  laser.  In  contrast  to  the  preceding  method, 
the  pump  signal  is  modulated  and  injected  into  the  laser  un¬ 
der  study  (test  laser)  while  the  ac  signal  appearing  across  the 
test  laser  terminals  is  monitored  by  a  lock-in  amplifier.  A 
modulating  frequency  of  a  few  kHz  is  sufficient  for  high 
signal  to  noise  ratia  When  using  a  pump  filtered  to  remove 
its  ASE,  the  magnitude  of  the  induced  ac  voltage  in  the  test 
laser  will  be  a  minimum  when  the  pump  laser  is  tuned  to  the 
transparency  energy.  Alternatively,  the  bias  of  the  laser  un¬ 
der  test  may  be  varied  while  the  pump  laser  is  fixed.  The 
latter  method  was  used  to  measure  the  transparency  energy 
and  loss  independently.  The  pump  sources  used  were  several 
single-mode  DFB  laser  operating  between  1514  and  1555 
nm.  Figure  6  shows  the  measured  ac  voltage  magnitude  ver- 


2950  J.  Appl.  Phys.,  Vol.  86,  No.  6,  15  September  1999 


Keating  et  al 


sus  current  for  three  of  the  pump  wavelengths  used.  The  ' 
negative  or  positive  sign  in  the  plot  is  used  to  indicate  the 
relative  phase.  This  plot  shows  that  as  the  pump  wavelength 
is  decreased,  the  current  which  has  that  wavelength  as  the 
transparency  point  increases,  which  is  the  correct  trend. 
From  our  measurements,  the  zero  crossings  are  obtained  by 
fitting  a  curve  to  the  measured  data  and  taking  the  fitted  zero 
value.  The  resulting  transparency  energies  are  compared  to 
those  of  the  other  methods  in  Fig.  3(a). 


III.  MANY-BODY  GAIN  THEORY 


The  optical  gain  with  the  Lorentzian  line  shape  function 
gives  an  anomalous  absorption  region  below  the  band  gap 
and  discrepancies  between  the  transparency  points  and  the 
Fermi-level  separation.  The  spontaneous  emission  transfor¬ 
mation  method  using  the  fundamental  relations  between  the 
spontaneous  emission  and  the  absorption  coefficients  guaran¬ 
tees  that  artifacts  related  to  the  Lorentzian  line  shape  are 
circumvented  independent  of  the  line  shape  functions 
(Lorentzian,  Gaussian,  or  others).  However,  the  slow  conver¬ 
gence  of  the  Lorentzian  still  leads  to  a  very  long  tail  of  the 
gain  spectra  into  the  band  gap.15  Here,  a  non-Markovian 
(Gaussian)  gain  model  with  many-body  effects  under  the 
spontaneous  emission  transformation  method  is  used  to  cal¬ 
culate  the  optical  gain,  and  we  will  show  that  this  model 
agrees  very  well  with  the  experimental  data  in  Sec.  IV.16,17 

Under  the  spontaneous  emission  transformation  method, 
the  optical  gain  with  non-Markovian  relaxation  and  many- 
body  effects  is  given  by16 


1  -  exp 


h<o-AF\  [JT0I 

kT  )  V  €  \ 


ReL[E/m(/:|i)]-Im^Jt]iImL[E/m(/:ll)] 
X  (l-Re^  +  am^ 


(8) 


where  <o  is  the  optical  angular  frequency,  fi0  is  the  vacuum 
permeability,  e  is  the  dielectric  constant,  &ft  is  the  in-plane 
.vector,  Lz  is  the  well  thickness,  |A/,J2  is  the  momentum 
matrix  element  in  the  strained  QWy  fj  and  fm  are  the  Fermi 
functions  for  the  conduction  band  states  and  the  valence 
band  states,  and  AF  is  the  difference  in  quasi-Fermi  energy 
levels  of  the  electrons  and  holes.  The  indices  l  and  m  denote 
the  electron  states  in  the  conduction  band  and  the  heavy  hole 
(light  hole)  subband  states  in  the  valence  band.  Also, 
Eim(kO  =  Eei{kO-E°m(ki)  +  Eg  +  AF5*+AFC//  is  the  nor¬ 
malized  transition  energy  between  electrons  and  holes, 
where  Eg  is  the  band  gap  of  the  material,  and  A Esx  and 
A  Ech  are  the  screened  exchange  and  Coulomb-hole  contri¬ 
butions  to  the  band  gap  renormalization.  The  factor  q(k\{)  is 
the  excitonic  or  Coulomb  enhancement  of  the  interband  tran¬ 
sition  probability.23,24  The  line  shape  function  [E(m{k^)]  is 
Gaussian  for  the  simplest  non-Markovian  quantum  kinetics 
and  is  given  by16,17 


FIG.  7.  The  measured  gain  spectra  (solid  curves)  and  theoretical  fits  (dashed 
lines)  using  a  many-body  gain  model  are  shown  for  both  TE  (upper  group  of 
curves)  and  TM  (lower  group)  polarizations. 
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The  intraband  relaxation  time  rin  and  correlation  time  rc  are 
assumed  to  be  85  and  20  fs  in  our  calculations,  respectively, 
for  the  best  fit  with  the  experimental  data.  The  valence  band 
structure  is  calculated  by  using  the  block-diagonalized  3X3 
Hamiltonian  based  on  the  k-p  method.25 


IV.  COMPARISON  OF  THEORETICAL  AND 
EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Figure  7  shows  the  comparison  between  experimental 
(solid  curves)  and  calculated  (dashed  curves)  TE  and  TM 
gain  spectra  for  the  lattice-matched  InGaAlAl  multi-QW  la¬ 
ser.  The  structural  details  are  summarized  in  Table  I.  Experi¬ 
mental  data  are  measured  at  injection  currents,  7=6,  8,  10, 
and  11.7  mA.  The  corresponding  carrier  densities  are  1.38, 
1.58,  1.72,  and  1.81  X  1018  cm“3,  respectively.  The  carrier 


TABLE  I.  Structure  of  the  laser  under  study. 


Well: 

Number  of  wells 
Material 

Strain 

Width 

PL  wavelength 

5 

InGaAlAs 
lattice  matched 

86  A 

w.56  fim 

SCH  width: 
Stripe  width: 
Cavity  length: 

600  A 

1 .46  fim 
638  /xm 

Barrier. 

Material 

Strain 

Width 

PL  wavelength 

InGaAlAs 
lattice  matched 
50  A 

1.21  fim 
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FIG.  8.  The  measured  intrinsic  loss  and  trend  for  the  fits  are  shown  in  (a)  for 
the  TE  polarization.  In  (b)  the  optical  confinement  factors  for  both  polariza¬ 
tions  are  plotted  as  a  function  of  wavelength,  and  in  (c)  the  measured  and 
predicted  net  modal  gain  versus  the  injection  current  density  are  shown. 


ment.  If  the  Lorentzian  line  shape  is  used  instead  of  the 
Gaussian  line  shape,  it  would  be  difficult  to  obtain  the  exact 
optical  losses  because  the  slow  convergence  of  the  Lorentz¬ 
ian  leads  to  a  very  long  tail  of  the  gain  spectra  into  the  band 
gap.  This  is  mainly  due  to  the  steepness  of  its  low-energy  tail 
in  the  Gaussian,  which  is  similar  to  the  line  shape  obtained 
by  considering  the  carrier-carrier  and  carrier-longitudinal  op¬ 
tical  (LO)  phonon  scatterings.7,8 

V.  SUMMARY  AND  CONCLUSIONS 

High  accuracy  measurements  of  gain,  loss  and  transpar¬ 
ency  energy  in  long-wavelength  semiconductors  based  on  a 
hybrid  approach  using  the  fundamental  relationship  between 
the  gain  and  the  spontaneous  emission  spectra  have  been 
demonstrated.  Our  hybrid  approach  for  the  gain  spectrum 
alleviates  many  of  the  problems  related  to  the  poor  signal  to 
noise  ratio  in  the  ASE  near  and  below  the  band  edge,  yield¬ 
ing  more  accurate  gain  data  in  this  region  than  any  of  the 
other  available  techniques.  The  values  of  loss  and  transpar¬ 
ency  energy  obtained  from  the  spontaneous  emission  trans¬ 
formation  method  agree  favorably  with  those  obtained  from 
other  techniques  such  as  the  measurement  of  polarization 
dependent  gain,  and  measurements  of  the  gain  compression 
and  electrical  characteristics  in  the  presence  of  an  optical 
pump  beam.  These  independent  measurements  of  transpar¬ 
ency  energy  and  loss  provide  validation  that  our  correlated 
gain  and  spontaneous  emission  spectra  have  a  high  degree  of 
reliability.  We  have  also  compared  the  measured  gain  spectra 
for  both  polarizations  with  calculations  using  a  non- 
Markovian  gain  model  with  many-body  effects  under  the 
spontaneous  emission  transformation  method.  This  model 
yields  very  good  agreement  to  the  measured  gain  spectra. 


densities  are  selected  to  give  the  best  fit  to  the  experimental 
data.  Recently,  Shtengel  et  al  showed  that  the  absolute  value 
of  the  optical  loss  increases  with  an  increasing  current.26 
Here,  we  obtain  the  optical  losses  for  the  TE  polarization 
equal  to  -21,  -22.3,  —23.9,  and  -24.9  cm"1  for  6,  8,  10, 
and  11.7  mA  injection  currents,  respectively.  These  values 
are  generally  within  1  cm"1  of  the  best  measured  values  and 
are  plotted  in  Fig.  8(a);  for  the  8  mA  spectra  the  losses  are 
within  3  cm"1.  The  theoretical  gain  spectra  show  reasonably 
good  agreement  with  the  experimental  results. 

The  optical  confinement  factor  Y  per  unit  width  of  the 
quantum  well  obtained  from  the  fitting  result  is  shown  in  Fig. 
8(b)  and  is  approximately  1.8X  10“ 3  nm"1  at  1.5  yum.  This 
value  is  in  good  agreement  with  that  reported  for  InP-based 
quantum  well  lasers.27  The  absolute  magnitude  of  the  TE  and 
TM  confinement  factors  differs  by  about  0.01;  the  TM  po¬ 
larization  being  less  confined.  Variation  in  the  confinement 
factor  with  respect  to  wavelength  follows  the  same  trend  for 
each  polarization.  Confinement  decreases  with  increasing 
wavelength,  since  the  ratio  of  the  active  region  length  or 
width  to  the  wavelength  is  decreasing.  In  Fig.  8(c),  the  pre¬ 
dicted  and  measured  values  of  peak  net  modal  gain  versus 
injected  current  density  also  show  excellent  agreement.  From 
these  results,  we  know  that  the  non-Markovian  gain  model 
with  many-body  effects  is  useful  for  comparison  with  experi- 
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Temperature  Dependence  of  Electrical  and  Optical 
Modulation  Responses  of  Quantum- Well  Lasers 
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Abstract — We  present  theory  and  experiment  for  high-speed 
optical  injection  in  the  absorption  region  of  a  quantum-well  laser 
and  compare  the  results  with  those  of  the  electrical  injection  in¬ 
cluding  carrier  transport  effect.  We  show  that  the  main  difference 
between  the  two  responses  is  the  low-frequency  roll-off.  By  using 
both  injection  methods,  we  obtain  more  accurate  and  consistent 
measurements  of  many  important  dynamic  laser  parameters, 
which  include  the  differentia]  gain,  carrier  lifetime,  K  factor,  and 
gain  compression  factor.  Temperature-dependent  data  of  the  test 
laser  are  presented  which  show  that  the  most  dominant  effect  is 
the  linear  degradation  of  differential  gain  and  injection  efficiency 
with  increasing  temperature.  While  the  A'-factor  is  insensitive 
to  temperature  variation  for  multiple-quantum -well  lasers,  we 
find  that  the  carrier  capture  time  and  nonlinear  gain  suppression 
coefficient  decreases  as  temperature  increases. 

Index  Terms — Differential  gain,  electrical  modulation,  optical 
injection,  transport  effect. 


I.  Introduction 

HIGH-SPEED  modulation  of  long-wavelength  semicon¬ 
ductor  lasers  is  of  major  technological  importance.  The 
continuing  development  of  high-bandwidth  fiber  optical  com¬ 
munication  systems  and  the  unceasing  demand  for  greater  data 
transmission  capacity  over  fiber-optic  cables  is  dependent  on 
the  superior  qualities  offered  by  semiconductor  laser  optical 
transmitters.  Long- wavelength  lasers,  as  studied  in  this  paper, 
are  especially  useful  for  long-haul  communications  since  the 
loss  in  optical  fibers  is  minimal  near  wavelengths  of  1.55  /un. 
Practical  communication  systems  based  on  these  lasers  also 
require  robust  temperature  performance  for  application  in 
ambient,  uncooled  environments  where  the  device  may  be 
required  to  operate  between  -40  °C  and  85  °C.  To  im¬ 
prove  these  devices,  it  is  therefore  important  to  study  their 
temperature-dependent  response  and  analyze  which  effects  are 
most  responsible  for  limiting  the  modulation  bandwidth. 

The  excessive  damping  that  limits  the  electrical  modulation 
bandwidth  of  quantum-well  (QW)  lasers  originates  from  a 
number  of  possible  mechanisms,  including  photon  lifetime  [1], 
carrier  capture  and  escape  [2],  [3],  carrier  diffusion  [4],  carrier 
heating  [5],  [6],  spectral  hole  burning  [7],  and  circuit  parasitics 
[8].  The  photons  generated  by  the  stimulated  recombination 
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of  electrons  and  holes  travel  in  the  cavity.  On  average,  they 
exist  for  a  certain  amount  of  time  and  disappear  from  the 
cavity  due  to  the  intrinsic  absorption  or  transmission  from 
the  mirror  facets.  Therefore,  the  photon  lifetime  constitutes 
an  intrinsic  limitation  on  the  modulation  response.  The  major 
sources  of  carrier  heating  in  semiconductor  lasers  are  injection 
heating,  stimulated  recombination  heating,  and  free  carrier 
heating.  These  changes  in  the  carrier  temperature  will  be 
reflected  in  the  changes  in  gain.  An  increase  of  a  few  degrees 
will  result  in  a  decrease  in  gain  of  several  percent,  which 
appears  in  the  rate  equations  as  a  nonlinear  gain  suppression 
coefficient.  Spectral  hole  burning  theory  also  predicts  an 
increase  in  the  nonlinear  gain  suppression  coefficient  with 
increasing  carrier  temperature.  Parasitic  effects  come  from  the 
bias  circuit  and  the  shunting  of  modulation  current  around 
the  active  layer,  which  will  cause  a  low-frequency  roll-off 
of  modulation  response.  The  carrier  diffusion,  capture,  and 
escape  times  are  usually  defined  to  characterize  the  carrier 
transport  processes  [9],  which  give  a  parasitic-like  roll-off 
and  are  indistinguishable  from  parasitic  effects.  In  QW  lasers, 
the  carrier  transport  effect  is  an  important  limit  for  multiple- 
quantum-well  (MQW)  laser  modulation  bandwidth.  In  general 
when  the  number  of  wells  increases,  the  modulation  bandwidth 
of  the  device  initially  increases,  but  it  is  ultimately  saturated  by 
the  carrier  transport  effect.  This  saturation  effect  was  observed 
in  gain-coupled  InGaAsP  distributed  feedback  (DFB)  lasers 
with  more  than  eight  QW’s  [10]. 

A  schematic  of  high-speed  modulation  by  electrical  injec¬ 
tion  and  optical  injection  is  shown  in  Fig.  1.  For  electrical 
modulation,  the  electrons  are  injected  from  the  outer  edge 
of  the  left  separate-confinement-heterostructure  (SCH)  region 
and  the  holes  from  the  outer  edge  of  the  right  SCH  region. 
The  injected  carriers  diffuse  through  the  SCH  region  and  are 
captured  into  the  QW’s  before  recombining  by  stimulated 
emission  processes.  Compared  with  electrical  modulation, 
optical  modulation  with  an  optical  energy  in  the  absorption 
range  of  the  QW  directly  produces  photon-generating  carriers 
inside  the  test  laser’s  active  region  via  injection  of  a  modulated 
laser  beam  through  one  of  the  test  laser’s  mirror  facets. 
Therefore,  the  majority  of  carriers  transporting  through  the 
SCH  region  is  not  required  for  lasing  action  although  the 
coupling  between  SCH  and  QW  states  still  exists  for  optical 
modulation.  In  this  way,  optical  modulaton  removes  the  severe 
low-frequency  roll-off  due  to  the  transport  and  parasitic  effects 
and  helps  to  clarify  the  intrinsic  response. 

In  this  paper,  the  high-speed  modulation  response  of  a  QW 
DFB  laser  operating  at  1.55  fjm  is  investigated  for  two  types 
of  modulation,  electrical  and  optical.  Electrical  modulation 
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Fig.  1.  Energy-band  QW  potential  diagram  for  high-speed  modulation  by 
optical  or  electrical  injections.  For  electrical  injection,  the  carriers  have  to 
transport  and  get  captured  by  the  QW’s.  For  optical  injection,  we  choose  the 
optical  pump  energy  such  that  carriers  are  generated  in  the  wells. 


is  performed  by  applying  a  microwave  signal  directly  on 
the  injection  current  of  the  test  laser.  Optical  modulation 
is  performed  by  injecting  light  from  a  second  pump  laser 
operating  at  1.3  /im  which  itself  is  being  directly  modulated 
electrically,  so  that  the  injected  light  is  absorbed  in  the  QW’s 
of  the  1.55-^m  test  laser.  Both  sets  of  responses  are  compared 
with  theory  in  order  to  extract  the  values  of  the  fundamental 
parameters  which  determine  as  well  as  limit  the  modulation 
bandwidth.  Optical  injection  at  this  wavelength  eliminates 
circuit  parasitic  and  carrier  transport  effects,  and  the  intrinsic 
laser  response  can  be  directly  measured  [11],  [12].  Our  exper¬ 
iment  differs  from  the  experiment  performed  in  [11],  where 
the  optical  injection  is  at  a  wavelength  close  to  that  of  the  test 
laser.  Therefore,  competition  of  carriers  and  gain  saturation 
have  to  be  considered  for  both  the  pump  and  test  wavelengths. 
We  inject  pump  signal  into  the  center  of  the  active  waveguide 
region  of  the  test  laser  to  avoid  the  spatial  effects  on  the 
optical  modulation  and  optimize  the  modulation  response.  It 
has  been  shown  [14]  that  the  optical  modulation  response 
can  be  influenced  by  the  spatial  characteristics  (the  spot  size 
and  its  position)  of  the  pump  light.  Our  setup  is  similar 
to  that  in  [15],  where  the  intrinsic  absorption  of  the  pump 
light  modulates  the  test  laser  output.  However,  little  work  has 
been  performed  so  far  on  directly  comparing  the  electrical 
and  optical  modulation  response  at  different  temperatures.  We 
show  data  for  temperature  dependence  of  both  optical  and 
electrical  injection.  To  analyze  the  basic  parameters  controlling 
the  modulation  response,  the  simplest  approach  is  to  use  rate 
equations  for  the  cavity  photon  density  and  modulated  carrier 
density.  Rate  equations  provide  a  reasonable  description  of 
most  of  the  observed  phenomena. 

This  paper  is  organized  as  follows.  In  Section  n,  the  optical 
modulation  theory  is  presented  and  compared  with  that  of 
electrical  modulation  including  transport  effects.  In  Section  HI, 


our  experimental  setup  is  described  and  comparisons  between 
optical  and  electrical  modulation  responses  are  made,  and 
the  extracted  high-speed  parameters  are  shown.  Finally,  the 
temperature  dependence  of  important  laser  parameters  are 
discussed  in  Section  IV. 


II.  Rate  Equations  for  Optical 
and  Electrical  Modulation 


Fig.  1  shows  a  typical  SCH  QW  laser.  We  will  model  the 
carrier  transport  (diffusion,  capture,  and  emission)  in  such 
a  laser  structure  and  study  its  effects  on  the  small-signal 
modulation  with  electrical  injection  and  optical  injection.  Our 
objective  is  to  calculate  the  small-signal  frequency  response 
of  the  test  laser  for  optical  modulation  from  the  pump  laser 
and  use  the  same  rate  equations  with  a  different  source  term  to 
derive  the  electrical  modulation  response  of  the  test  laser.  In 
general,  transport  effects  for  electrical  microwave  modulation 
can  be  modeled  by  taking  into  account  the  carrier  density  in 
the  SCH  region,  the  carrier  density  in  the  well  region,  and  the 
photon  density  separately  [9].  Coupling  of  the  carrier  density 
in  the  barrier  states  above  the  QW’s  to  the  carrier  density 
in  the  QW’s  is  modeled  by  two  terms  representing  carrier 
capture  and  escape  into  or  from  the  wells,  respectively.  In 
this  case,  three  rate  equations  are  needed.  The  source  term 
enters  through  the  injection  current  in  electrical  modulation. 
The  model  considers  carrier  injection  from  the  outer  edges 
of  the  SCH  region,  diffusion  across  the  SCH  region,  and 
the  subsequent  capture  and  emission  of  carriers  by  the  QW. 
For  optical  injection  using  an  external  pump  laser,  the  pump 
photon  density  Sp  acts  as  the  source  term.  Since  the  optical 
energy  of  the  pump  laser  determines  whether  the  photons  are 
absorbed  in  the  well  or  in  the  barriers,  we  choose  the  pump 
wavelength  to  be  longer  than  the  bandgap  wavelength  of  the 
barriers  and  shorter  than  the  bandedge  wavelength  of  the  wells 
so  that  absorption  occurs  only  in  the  wells.  The  rate  equations 
for  both  cases  are  written  as 


dNb  -Vj)  Nf)  ^  NwVw 

dt  ^7^6  7 ~b  7" bu;  7 

dNw  _  NbVb  Nw  Nw  vgG(Nw)S 
dt  nwVw  rw  rwi  14 -  eS  V90tp  p 

dS  _  rvgG(Nw)S  S 
dt  1 4-  eS  rp 


(1) 

(2) 

(3) 


where  q  is  the  electron  unit  charge,  G(NW )  is  the  optical  gain 
at  the  carrier  concentration  Nw  in  the  bound  states  of  the 
wells,  Nb  is  the  carrier  density  in  the  barrier  (continuum)  states 
including  the  SCH  and  active  layers,  5  is  the  photon  density  of 
the  test  laser,  e  is  the  nonlinear  gain  suppression  coefficient,  T 
is  the  optical  confinement  factor,  rp  is  the  photon  lifetime,  7]^ 
is  the  injection  efficiency,  I  is  the  injection  electrical  current, 
Vb  is  the  volume  of  the  SCH  and  active  region,  Vw  is  the 
volume  of  the  well  region,  n  is  the  carrier  recombination 
lifetime  in  the  barrier  region,  rw  is  the  carrier  recombination 
lifetime  in  the  well,  Tbw  is  the  effective  carrier  diffusion  across 
the  SCH  region  and  capture  time  by  the  wells,  and  rwb  is  the 
thermionic  emission  and  carrier  diffusion  time  from  the  well  to 
the  barrier  states.  The  final  term  in  (2)  is  a  source  term  which 
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represents  the  photon  generation  due  to  optical  pumping  at  the 
wavelength  (1.3  fmi)  where  the  optically  injected  photons  are 
absorbed  in  the  well  regions.  Sp  is  the  pump  photon  density, 
and  ap  is  the  absorption  at  the  pump  wavelength.  The  above 
model  is  often  referred  to  as  the  reservoir  model  [9]  and  is 
equivalent  to  models  which  incorporate  additional  effects  such 
as  diffusive  transport  [14].  A  spontaneous  emission  term  has 
been  ignored  in  (3)  for  above  threshold  operation. 

Assuming  a  small-signal  optical  injection  due  to  the  external 
pump  laser 

Sp(t)  =  Sp0  +  8p(w)e?ut  (4) 


where  the  test  laser  is  biased  at  a  dc  current  I(t)  =  Jo  above 
threshold.  The  responses  can  be  solved  by  assuming 


S(t)=S0  +  8(w)e**  (5) 

Nb(t)=Nb0  +  nb(u)e>ut  (6) 

Nw(t)=Nw0  +  nw{u)e>ut  (7) 

and 

G(NW )  =  G(Nw0)  +  g\Nw  -  Nw0)  (8) 


where  gf  is  the  differential  gain  of  the  test  laser,  and  Nwo  is 
steady-state  carrier  density  in  the  QW. 

Substituting  (4) — (8)  into  (l)-(3),  the  steady-state  quantities 
can  be  obtained  by  setting  the  time-dependent  terms  to  zero. 
Then,  the  small-signal  equations  for  the  time-dependent  terms 
give 


n  Tib  Tib  TlwVw 

JUJUb  =0 - I - — 

Tfc  ^~btu  1~wb  *  b 


(9) 


junw 


ribVb 
Tbxv^w 

g  Sqtiw 


Tlw 

TVn 


Tin 


—  Vt 


G0s 


GoSqcs 


+ 


jus  -  Tvg 


1  +  eSo 

G0S0es 


T% vb 

4“  T)  g  Ctp  Sp 

cfS0n 


ffll  +  e5o  (l  +  eS0)2 
(10) 


+  : 


L  (14-  e50)2  1  +  cSq 


(11) 


where  the  steady-state  gain-loss  relations  rV5G0/(l-f  eSo)  = 
l/rp,  and  the  Taylor  expansion  1/(1  +  eS)  =  1/(1  4-  eSo)  — 
e5/(l-h  e5o)2  are  used.  This  set  of  small-signal  equations  can 
be  reduced  by  eliminating  rib  and  nw  to  give  a  relationship 


shown  in  (12),  at  the  bottom  of  the  page,  between  the  test 
laser  signal  s{u)  and  the  pump  signal  sp{uS). 

The  optical  modulation  response  can  be  written  in  a  nor¬ 
malized  form 


where 


M{u)  ul  .  .  1X 

777 XT  =  5 - TT~- —  (optical) 

M(  0)  u*  joj'y 

2  _  M/x)Sq  f1  e  \ 
r  rp(l  4-  eSo)  \  vggfTw) 

_(vg9'/x)So  .  ^0  1 

;  1  +  eS0  (14-  eS0)rp  xrw 


(13) 


(14) 

(15) 

(16) 


For  a  small-signal  electrical  injection,  I(t)  =  Io+i^e^, 
with  no  external  pump,  Sp(t)  =  0.  The  resulting  modulation 
response  is  given  by  [9]  and  [21]  in  (17),  shown  at  the  bottom 
of  the  page,  which  can  be  written  in  the  form 

M(u)  __ 

M( 0)  ”  (14-  j<jJTbw)(u)r  —(jS2^ 


(electrical)  (18) 


where  Tbw  is  a  constant  which  contributes  to  low-frequency 
roll-off.  The  interpretation  of  this  effect  is  not  straightforward. 
The  roll-off  in  the  modulation  response  derived  above  is 
equivalent  to  the  roll-off  caused  by  a  simple  electrical  parasitic, 
such  as  a  capacitance  in  parallel  to  the  device.  However,  in 
real  devices,  the  effective  capacitance  may  be  bias-dependent 
since  the  device  capacitance  is  related  to  the  storage  of 
charge  in  the  forward-biased  junction.  The  charge  storage  itself 
results  from  a  combination  of  effects  including  transport  and 
carrier  heating  and  is  difficult  to  model  in  terms  of  a  small 
number  of  rate  equations  [15]-[17].  Despite  these  difficulties, 
the  modulation  response  can  be  matched  very  well  for  the 
experiments  here  with  a  single  total  roll-off  frequency  fp  for 
a  given  temperature.  If  the  transport  effect  dominates,  the  roll¬ 
off  frequency  of  the  transport  effect  is  fp  =  1/(2^^).  In 
this  case,  we  can  measure  the  carrier  diffusion  and  capture 
time  of  the  QW. 

Because  the  stimulated  emission  occurs  in  the  QW’s,  the 
carrier  lifetime  in  the  SCH  or  barrier  region  is  usually  very 


M(u)  = 


jOjO 

Sp(u) 

_ v^ryso _ 

_ juUw  +  Thy/n  \  1  S0g'vg  1  /  eSp  \  Vg^Sp 

/  \  jUTwb(jUTbw  +  l+Tbw/Tb) )  Tw  1  +  eS0J  V  Tj,(l  +  <zS0) )  Tp(l  +  eS0)2 


(12) 


s(u ) 

M(u)  =  -7-f 
t{u) 


Vinj 

qVw 


VgVg1 


So 


(1  +  €5q) 


(1  +  jUTbw) 


|  JW^1  + 


Tbw 


Twb(jvTbw  +  1) 


)  +  ^ 
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Spg'vg 


1  +  eSo 
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eS0 
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vgg'Sp 
rp(l  +  eS0)2 
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TABLE  I 

Summary  of  Electrical  and  Optical  Modulation  Responses 


TABLE  n 

Structure  of  the  Test  Laser 


Electrical  Pumping  and  1.3  Optical  Pumping 


Electrical _ Optical _ 

M(cj)  _  _ <4  _  M(uS)  _  ^ _ 

M(0)  “  (l+ju>r^ul)(u>2  — cj2-f-iu?7)  M(0)  w2—  u>*-f 

t  ~  ™  +  r^l+lso'J  (Tp  +  &) 
_ *=1  +  ^ _ 


long,  so  we  set  «C  1  in  (12).  The  modulation  responses 

for  both  injections  are  summarized  in  Table  I.  Thus  the  com¬ 
parison  between  the  optical  and  electrical  modulations  is  clear. 
The  relaxation  frequency  ur  —  2nfr  and  damping  factor  7  for 
both  cases  are  the  same.  The  transport  factor  \  —  ^~\-Uw/Twb 
has  also  been  introduced  to  the  equations.  The  only  difference 
between  the  two  responses  is  that  optical  injection  has  no 
low  frequency  roll-off,  which  corresponds  to  transport  effects. 
There  would  also  be  a  difference  when  the  number  of  wells  is 
varied.  Thus,  optical  modulation  gives  an  intrinsic  modulation 
response  of  a  laser  since  the  injected  pump  light  is  chosen 
such  that  the  pump  light  is  absorbed  only  in  the  well  region. 

For  optical  injection,  the  carriers  are  injected  directly  into 
the  QW  region,  which  removes  the  requirement  of  carrier 
transport  from  the  SCH  region  to  the  QW  region  before 
lasing.  In  other  words,  optical  injection  removes  the  low- 
frequency  roll-off  caused  by  carriers  transporting  through  SCH 
region.  However,  the  SCH  and  QW  regions  are  still  coupled 
together.  There  are  always  carriers  being  captured  by  and 
escaping  from  the  QW’s.  The  rate  of  capture  and  escape 
will  influence  the  modulation  response  intrinsically.  Therefore, 
because  of  the  SCH  structure,  carrier  diffusion  capture-escape 
still  affects  the  optical  modulation  response,  which  comes  in 
through  the  x  factor.  Thus,  the  effective  differential  gain  is 
reduced  by  a  factor  of  x  f°r  the  same  photon  density.  This 
reduction  is  present  even  in  the  absence  of  low-frequency  roll¬ 
off  in  the  optical  injection,  which  results  in  the  reduction  of 
the  resonance  frequency.  The  effective  carrier  recombination 
lifetime  in  the  well  is  also  increased  by  a  factor  of  and 
the  nonlinear  gain  suppression  coefficient  remains  unchanged. 
Another  important  factor  in  high-speed  modulation  is  the 
K  factor,  which  is  the  slope  of  the  damping  factor  versus 
relaxation  frequency  squared 


(19) 

K~4n  {Tp+vgg'/X)' 

(20) 

The  K  factor  is  usually  used  to  determine  the  maximum 
possible  modulation  bandwidth  /maLX  =  2n  \f2jK .  A  small  K 
factor  means  a  large  laser  bandwidth.  Even  if  we  use  optical 
injection  modulation,  the  carrier  capture  and  escape  processes 
still  increase  the  K  factor  and  limit  the  maximum  modulation. 
The  above  parameters  characterize  the  intrinsic  modulation 
response  of  semiconductor  lasers. 


Well: 


Barrier: 


Number  of  wells  7  SCH  width:  705  A 

Material  InGaAsP  Stripe  width:  1.2  fxm 

Strain  —1.6%  comp.  Cavity  length:  400  //m 

Width  70  A 

PL  wavelength  1.5564  /on _ 
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Strain  lattice-matched 

Width  100  A 

PL  wavelength  1 .255  fim 
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Fig.  2.  Schematic  diagram  of  electrical  and  optical  modulation  experiment. 
Path  A  is  for  electrical  injection  experiment.  Path  B  is  for  optical  injection 
experiment  with  a  pump. 


in.  Experimental  Results 

The  high-speed  modulation  response  of  a  buried  heterostruc¬ 
ture  MQW  DFB  laser  was  measured.  The  composition  of  the 
undoped  active  region  is  described  in  Table  H.  The  barrier 
photoluminescence  wavelength  peak  occurs  at  1.255  jum. 
Therefore,  pump  light  at  1.3  [im  will  be  absorbed  in  the  QW’s. 

A  schematic  diagram  of  the  modulation  response  exper¬ 
iments  is  shown  in  Fig.  2.  The  laser  is  held  at  constant 
temperature  and  measurements  are  made  at  20°C,  25°C,  30°C, 
35°C,  and  40  °C.  In  the  electrical  modulation  experiments, 
shown  as  path  A  in  Fig.  2,  the  HP  8510  network  analyzer 
provides  a  small  microwave  signal  (0  dBm)  at  frequencies 
swept  from  45  MHz  to  10  GHz,  which  is  coupled  to  the  test 
laser  electrodes  through  a  bias-T  and  a  high-speed  probe.  The 
laser  light  is  coupled  to  a  fiber  using  a  lensed  fiber-optical 
interface  and  travels  through  an  optical  isolator  before  being 
measured  by  a  high-speed  (29  GHz)  photodetector.  The  small- 
signal  response  is  increased  by  an  18-dB-gain  RF  amplifier 
before  entering  the  network  analyzer,  which  measures  the 
magnitude  of  the  modulation  response  |M(cu)/M(0)|2.  The 
data  are  averaged  to  reduce  noise.  The  optical  response  mea¬ 
surements  shown  as  path  B  in  Fig.  2  are  similar  to  electrical 
modulation  experiments,  except  that  a  1 .3-^m  pump  laser  with 
a  wide  bandwidth  is  modulated  electrically,  and  its  modulated 
optical  output  is  injected  into  the  test  DFB  laser  facet.  The 
modulation  of  the  pump  laser  is  recorded  at  a  fixed  bias  and 
stored  for  the  calibration  of  the  optical  response  of  the  1.55-^m 
test  laser.  The  pump  light  is  completely  absorbed  in  the  test 
laser,  so  the  light  coupled  out  of  the  test  laser  does  not  require 
additional  filtering  of  the  pump  light. 
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Fig.  3.  Normalized  modulation  response,  101og(|M(w)/M(0)|2),  due  to 
electrical  injection  and  optical  injection  at  1.3  /i m.  The  solid  and  dashed  lines 
are  theoretical  modulation  responses  using  (13)  and  (18),  and  symbols  are 
measurement  results. 


Current  (mA) 


Fig.  4.  The  total  roll-off  frequency  (/p  =  1/27 tt^)  versus  the  bias  current 
at  different  temperature.  The  solid  line  is  the  average  value  of  roll-off 
frequency. 


The  electrical  and  optical  modulation  responses  are  first 
compared  at  a  fixed  temperature  of  20  °C.  Fig.  3  shows 
the  least-square  fits  of  the  measured  modulation  responses 
using  the  theory  derived  above  for  electrical  and  optical 
injections.  The  fits  are  used  to  extract  the  differential  gain 
and  nonlinear  gain  suppression  factor.  In  general,  the  optical 
injection  measurements  exhibit  more  noise  than  the  electrical 
injection  measurements,  since  the  system  losses  are  higher 
for  the  optical  injection  experiments.  However,  the  optical 
modulation  responses  have  less  parameters  to  fit,  which  means 
less  freedom  of  choosing  a  number.  On  the  other  hand,  the 
electrical  modulation  responses  usually  have  a  higher  signal- 
to-noise  ratio  and  more  parameters  to  be  handled  than  the 
optical  modulation.  Thus,  simultaneously  fitting  electrical  and 
optical  modulation  responses  can  compensate  each  method’s 
shortage  and  obtain  more  accurate  laser  parameters.  Also,  by 
comparing  the  two  responses,  the  roll-off  frequency  fp  can  be 
obtained  more  accurately.  Both  modulation  responses  show 
clearly  that  the  relaxation  peaks  increase  in  frequency  as  the 
dc  bias  is  increased.  The  modulation  responses  due  to  optical 
injection  exhibit  a  slower  roll-off  at  high  frequencies,  which 
indicates  the  absence  of  low-frequency  pole  on  the  optical 
response.  Also,  the  peak  of  the  optical  pump  response  is 
generally  higher  than  that  of  the  electrical,  and  the  differences 
between  the  two  responses  increase  with  increasing  current. 
For  35-mA  current  injection,  electrical  modulation  has  very 
obvious  roll-off  because  the  roll-off  frequency  is  3.52  GHz 
at  this  temperature  (20  °C).  Fig.  4  shows  that  the  roll-off  fre¬ 
quency  fp  at  various  temperatures  from  20  °C  to  40  °C  is  inde¬ 
pendent  of  the  current  bias.  As  we  discussed  before,  the  shunt 
capacitance  of  parasitic  effects  may  be  bias-dependent,  and  the 
roll-off  frequency  of  the  transport  effect  is  only  determined  by 
carrier  diffuse-capture  time  nw.  Although  we  cannot  exclude 
the  parasitic  effects,  we  still  can  see  that  the  transport  effects 
are  more  important  than  the  parasitic  effects  for  this  DFB  laser. 

In  Fig.  5(a),  the  relaxation  frequency  squared  is  plotted 
versus  optical  power.  The  slope  of  the  linear  fit  will  be  used 
later  to  extract  the  differential  gain.  The  damping  factor  is 
calculated  using  extracted  parameters  from  the  modulation 
responses  and  is  plotted  in  Fig.  5(b).  The  slope  of  the  linear  fit 
at  large  relaxation  frequency  is  the  K  factor.  The  deviation  at 
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Fig.  5.  (a)  The  relaxation  frequency  squared  versus  optical  power  at  20  °C. 
The  line  is  the  least-square  fit  to  the  data  using  (14).  (b)  The  damping  factor 
versus  the  relaxation  frequency  squared  at  20  °C.  The  dash  line  is  the  linear 
fit  for  large  relaxation  frequency,  with  a  slope  equal  to  the  K  factor. 


the  low-frequency  end  between  the  data  and  the  linear  dashed 
line  is  due  to  the  negligence  of  the  spontaneous  emission 
term.  If  we  include  the  spontaneous  emission  factor  RsP ,  the 
damping  factor  should  have  an  additional  term  R^p/S  which 
can  be  neglected  at  high  photon  density  S  (corresponding  to 
a  large  relaxation  frequency).  The  complete  set  of  extracted 
parameters  at  20  °C  is  listed  in  Table  IK.  The  differential  gain 
and  nonlinear  gain  suppression  coefficient  are  well  within  the 
values  typically  found  in  the  literature. 

IV.  Temperature  Dependence 

It  is  important  to  determine  which  factors  affecting  the 
modulation  responses  are  most  sensitive  to  temperature.  The 
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TABLE  III 

High-Speed  Laser  Parameters  and  Extracted  Values  at  20  °C,  Using  Simultaneous  Fitting  of  the  Electrical  and  Optical  Modulation  Responses 


Parameter 

Symbol 

Value 

Error 

Carrier  lifetime 

Tw 

0.37  ns 

±0.05  ns 

Effective  differential  gain 

a'/x 

1.32  X  10“15  cm2 

±0.01  x  10-15  cm2 

K  factor 

K 

0.20  ns 

±0.06  ns 

Nonlinear  gain  suppression 
coefficient 

€ 

3.11  x  10~17  cm3 

±0.5  X  10~17  cm3 

Photon  lifetime 

TP 

2.4  ps 

±0.2  ps 

Intrinsic  loss 

20  cm  -1 

±2  cm-1 

Mirror  loss 

a 

30  cm  _1 

±3  cm-1 

Group  velocity 

V9 

8.5  X  109  cm/s 

±0.1  x  109  cm/s 

Injection  efficiency 

^inj 

0.32 

±0.1 

Width 

W 

1.25  ftm 

Well  width 

Lw 

7  nm 

Barrier  width 

Lh 

10  nm 

SCH  width 

•f'sch 

70.5  nm 

Number  of  wells 

Nw 

7 

Active  volume 

V 

1.2  x  10-10  cm3 

Optical  coupling  constant 

Vc 

0.10 

±0.05 

Optical  confinement  factor 

r 

0.06 

±0.01 

Fig.  6.  The  injection  efficiency  and  threshold  current  as  a  function  of 
temperature.  The  efficiency  degrades  severely  with  temperature.  The  solid 
line  shows  the  trend.  The  dashed  line  is  fitting  of  threshold  current  using  the 
formula  Ith  =  /0eT/T°,  where  T0  =  24  K. 

fitting  procedure  described  above  is  applied  to  each  set  of  mod¬ 
ulation  responses  at  each  temperature.  In  Fig.  6,  the  injection 
efficiency  and  threshold  current  of  the  test  laser  are  plotted 
versus  temperature.  This  plot  is  obtained  from  measured  L-I 
curves.  The  injection  efficiency  is  defined  as  the  fraction 
of  current  above  threshold  which  results  in  the  radiative 
recombination.  A  fraction  of  the  injected  current  in  InGaAsP 
lasers  is  lost  (by  nonradiative  recombination  or  carrier  leakage) 
at  high  temperatures,  and  this  fraction  does  not  contribute 
to  radiative  recombination  or  optical  gain.  The  injection  ef¬ 
ficiency  degrades  relatively  severely  with  temperature,  by  a 
factor  of  more  than  two  over  the  30°  temperature  increase 
(about  10%  in  absolute  temperature).  The  injection  efficiency 
acts  to  change  the  modulation  response  primarily  by  reducing 
the  overall  magnitude  of  the  response.  Indirectly,  however,  the 
wasted  carriers  may  contribute  to  the  forward  bias  capacitance, 
further  influencing  the  modulation  response.  An  exponential 
temperature  dependence  of  the  threshold  current  is  generally 
used,  and  the  overall  characteristic  temperature  T0  is  24  K  for 
this  laser,  which  indicates  a  large  dependence  on  temperature. 


Fig.  7.  The  effective  differential  gain  (gf /x)  as  a  function  of  temperature, 
showing  strong  degradation  with  temperature,  x  =  1  +  Tbw/Twb  accounts 
for  carrier  transport/capture  effect. 

A  higher  carrier  density  is  necessary  to  achieve  the  threshold 
gain  condition  with  increasing  temperature.  This  results  in  the 
nonradiative  recombination  current  being  a  large  fraction  of  the 
total  current.  Several  mechanisms  have  been  proposed  to  ex¬ 
plain  the  observed  high-temperature  sensitivity  of  the  threshold 
current  of  InGaAsP  lasers.  These  are  carrier  leakage  over  the 
heterojunction,  Auger  recombination,  and  intervalance  band 
absorption. 

The  effective  differential  gain  (g'/x)  al so  shows  a  strong 
temperature  dependence  and  is  plotted  in  Fig.  7.  Over  the 
temperature  range  plotted,  the  differential  gain  decreases  by 
a  factor  of  approximately  two,  which  shows  a  linear  function 
of  the  temperature  [27].  This  effect  is  due  to  the  temperature 
dependence  of  the  Fermi  distribution,  the  carrier  density,  and 
the  intrinsic  transition  linewidth  broadening.  However,  it  is 
difficult  to  isolate  the  significance  of  each  quantity  on  the 
temperature  dependence  of  a  laser  because  so  many  physical 
phenomena  are  involved.  It  can  be  interpreted  as  follows. 
The  broadening  of  the  Fermi  occupation  probability  function 
with  increasing  temperature  spreads  the  carriers  over  a  larger 
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energy  range  for  a  given  overall  carrier  density.  The  result 
is  a  lower  spectral  concentration  of  inverted  carriers,  which 
leads  to  a  broadening  and  flattening  of  the  gain  spectrum. 
Thus,  the  gain  is  higher  at  lower  temperatures  for  the  same 
carrier  density,  and  the  carrier  density  required  to  achieve 
a  particular  gain  increases  with  temperature.  In  addition, 
thermionic  emission  of  electrons  from  the  QW  to  the  barrier 
region  is  also  enhanced  at  high  temperatures.  These  effects 
contribute  to  strong  recombination  outside  the  QW  and  carrier 
leakage.  Since  the  increased  carrier  population  in  the  barrier 
region  does  not  contribute  to  the  optical  gain,  the  differential 
gain  decreases  rapidly  at  high  temperatures.  This  decrease 
in  gain  strongly  reduces  the  modulation  bandwidth  [13]. 
The  differential  gain  plays  a  central  role  in  determining  the 
fundamental  frequency  response  of  semiconductor  lasers  since 
the  intrinsic  direct  modulation  speed  varies  as  the  square  root 
of  the  differential  gain  at  the  operating  carrier  density  and 
wavelength. 

In  Fig.  8(a),  the  carrier  lifetime  rw  in  the  QW’s  is  plot¬ 
ted  and  shows  a  decrease  with  increasing  temperature.  The 
carrier  lifetime  is  expected  to  shorten  with  increasing  temper¬ 
ature,  since  in  1.55-^m  devices  the  carrier  recombination  is 
dominated  by  Auger  recombination  [18],  which  is  strongly 
temperature-dependent.  In  Fig.  8(b),  the  total  roll-off  fre¬ 
quency  is  plotted  versus  temperature,  showing  an  increase  with 
increasing  temperature.  Usually,  the  low-frequency  roll-off  is 
a  combination  of  the  effects  of  the  shunt  capacitance,  external 
circuit  parasitic,  and  internal  transport  effect.  As  for  the 
transport  effect,  the  carrier  diffusion  and  capture  time  across 
the  SCH  region  (tiw)  is  usually  limited  by  the  hole  diffusion 
time,  which  is  insensitive  (for  MQW’s)  or  slightly  increasing 
(for  single-QW’s)  according  to  temperature  increments  from 
200  to  350  K  [9],  [21].  If  we  assume  the  carrier  transport  effect 
is  dominant,  we  can  calculate  the  carrier  diffusion-capture 
time  Tbw>  which  is  also  shown  in  Fig.  8(b).  However,  our 
data  shows  a  decrease  of  nw  with  increasing  temperature. 
This  can  be  explained  as  follows.  First,  in  simple  view,  the 
carrier  diffusion-capture  time  nw  is  inversely  proportional 
to  the  hole  diffusion  constant  D  =  ksT^x/q  where  the 
mobility  fx  is  also  temperature-dependent,  k&  is  the  Boltzmann 
constant.  Actually  the  mobility  \i  decreases  with  increasing 
temperature.  The  total  changes  in  diffusion  constant  D  are 
determined  by  changes  in  fxT.  To  be  more  exact,  the  carrier 
diffusion  is  usually  referred  to  the  diffusion  capacitance  of 
the  test  laser,  which  has  a  separate  temperature  dependence 
not  limited  to  constant  D  in  obvious  ways  [16].  Second, 
the  transport  effect  may  not  play  the  only  role  in  the  roll¬ 
off  frequency.  The  electrical  parasitic  effect  can  also  account 
for  the  low-frequency  roll-off.  However,  the  external  circuit 
parasitic  will  not  change  significantly  with  laser  temperature. 
Thus,  the  capacitance  of  the  test  laser  should  respond  to 
the  increase  of  roll-off  frequency.  This  trend  indicates  that 
the  capacitance  decreases  with  increasing  temperature.  This 
is  reasonable  since  the  charge  storage  capacitance  decreases 
with  decreasing  carrier  lifetime.  Generally,  because  of  this 
parasitic  frequency,  the  roll-off  frequency  cannot  be  used  to 
directly  extract  the  exact  value  of  carrier  diffusion-capture 
time.  Typical  values  for  the  well  capture  time  are  on  the 
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Fig.  8.  (a)  The  carrier  lifetime  is  plotted  as  a  function  of  temperature, 
showing  a  large  decrease  as  the  temperature  increases,  (b)  The  roll-off 
frequency  is  plotted  as  a  function  of  temperature,  which  increases  with 
temperature.  The  carrier  diffusion  and  capture  time  is  also  plotted  in  (b), 
which  decreases  with  temperature,  (c)  The  K  factor  is  plotted  as  a  function 
of  temperature,  which  is  insensitive  to  temperature,  (d)  The  nonlinear  gain 
suppression  coefficient  is  plotted.  The  line  in  each  figure  shows  the  trend 
only. 


order  of  tens  of  picoseconds  [9],  [14],  [21].  Fig.  8(c)  shows 
that  the  K  factor  is  insensitive  to  temperature,  which  agrees 
with  the  previous  literature  [21].  The  more  QW’s  there  are, 
the  more  independent  of  the  K  factor  on  temperature  [21]. 
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For  single-QW  lasers,  the  K  factor  can  vary  by  two  times 
between  200-350  K  [4].  The  temperature  insensitivity  of  the  K 
factor  implies  that  the  differential  gain  and  the  nonlinear  gain 
suppression  coefficient  have  a  similar  temperature  dependence 
and  cancel  each  other  to  maintain  a  temperature-insensitive 
e/{gf /x)  ratio  in  InGaAsP-InP  QW  materials.  This  means 
the  nonlinear  gain  suppression  coefficients  should  decrease 
with  increasing  temperature.  The  nonlinear  gain  suppression 
coefficients  of  different  temperatures  are  shown  in  Fig.  8(d), 
which  are  particularly  sensitive  to  the  extracted  value  of  the 
damping  factor  and  have  large  error  variations.  However, 
the  decrease  in  nonlinear  gain  suppression  coefficient  with 
increasing  temperature  still  can  be  seen.  The  nonlinear  gain 
saturation  results  from  a  variety  of  factors,  but  it  is  considered 
to  arise  primarily  through  dynamic  carrier  heating  [5],  [6] 
and  perhaps  spectral  hole  burning  [7].  If  the  K  factor  is 
constant,  it  means  the  maximum  bandwidth  /max  is  insensitive 
to  the  temperature.  The  lower  temperature  will  not  improve 
the  potential  bandwidth  of  QW  lasers.  Indeed,  it  gets  worse 
because  of  the  lower  roll-off  frequency.  On  the  other  hand, 
the  injection  efficiency  and  high  threshold  will  deteriorate  the 
optical  signal  at  higher  temperature.  Therefore,  designing  a 
high-speed  laser  is  a  tradeoff  among  all  aspects. 

V.  Conclusions 

The  temperature  dependence  of  the  electrical  and  opti¬ 
cal  modulation  responses  have  been  investigated,  and  rate- 
equation  theoretical  models  show  excellent  agreement  with  the 
data.  The  major  difference  between  the  two  responses  is  the 
low-frequency  roll-off,  which  is  constant  with  respect  to  bias 
current  at  a  given  temperature.  The  carrier  diffusion-capture- 
escape  time  influences  the  intrinsic  modulation  response  of 
QW  lasers  through  the  transport  factor  x-  We  fit  electrical  and 
optical  modulation  responses  at  the  same  time  at  each  current 
and  each  temperature.  The  combination  of  two  measurements 
allows  us  to  obtain  a  set  of  more  accurate  laser  parameters. 
Temperature-dependent  laser  parameters  are  shown  in  this 
paper.  When  temperature  increases,  the  most  dominant  effect 
is  the  linear  degradation  of  differential  gain,  and  reduction 
of  injection  efficiency  and  carrier  lifetime  in  the  QW’s.  The 
temperature-dependent  response  of  all  parameters  closely  fol¬ 
lows  physically  reasonable  trends. 
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We  present  a  comprehensive  theoretical  model  with  experimental  results 

feedback  laser  with  two  pump  modes.  Here  the  twin  lasing  modes  serve  as  pump  waves  for  the  tom  wave 
Sig  wWch  takes  placeP  within  the  laser  cavity.  We  show  what  we  ^^^tT^orfer  ^ 
ments  of  the  two-pump  four-wave  mixing  conversion  efficiency  spectrum  and  show  that,  to  first 
no  fimdamenKpendence  of  the  conversion  efficiency  on  pump  separation.  This  result  mdicates  s^mfic^t 
improvement  in  the  probe  and  conjugate  wavelength  separation  and  in  the  conversion  efficiency.  Our 'the 
retical  model  agrees  quite  well  with  the  measured  data  and  predicts  the  key  experimentally  observed  phenom¬ 
ena  winning  therelative  invariance  of  the  conversion  efficiency  on  pump  separation,  and  the  ^°nge  an 
ment  ofSn  conjugate  waves  when  the  probe  detuning  is  close  to  the  pump  separation.  ©  2000  Optical 
Society  of  America  [S0740-3224(00)01101-2] 
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1.  INTRODUCTION 

Much  of  the  recent  research  on  four-wave  mixing  (FWM) 
has  been  driven  by  its  potential  as  an  all-optical  wave¬ 
length  conversion  technique.1  Many  groups  have  pub¬ 
lished  studies  of  FWM  in  semiconductor  optical  amplifiers 
(SOA’s),2-4  where  an  external  pump  wave  and  an  external 
probe  wave  are  injected  into  the  amplifier.  Also,  other 
groups  have  studied  FWM  in  single-mode  distributed 
feedback  (DFB)  lasers,6-9  where  the  lasing  mode  serves 
as  the  pump  wave  and  an  external  probe  is  injected  into 
the  Inning  cavity.  One  of  the  largest  challenges  facing 
the  practical  use  of  this  technology  has  been  the  poor  con¬ 
version  efficiency  for  useful  wavelength  shifts.  Only  very 
recently  has  significant  improvement  been  achieved,  with 
the  use  of  extremely  long  SOA  s.4 

A  clever  solution  to  improve  the  conversion  efficiency 
for  larger  wavelength  shifts  was  proposed  several  years 
ago10  by  introduction  of  a  second  pump  wave  with  a  de¬ 
tuning  from  the  first  pump  wavelength.  The  dynamic 
grating  formed  by  the  beating  of  the  first  pump  with  the 
probe  wave  causes  the  second  pump  to  scatter  into  a 
shifted  conjugate  wave.  The  extra  wavelength  shift 
gained  in  this  approach  is  exactly  the  wavelength  separa¬ 
tion  of  the  two  pump  waves.  Subsequent  reports,  includ¬ 
ing  experimental  data11  (although  few  points  of  data  were 
offered)  and  theory,12'13  indicated  that  the  conversion  ef¬ 
ficiency  is  relatively  independent  of  the  pump  wavelength 
separation,  being  indirectly  affected  only  through  the  dis¬ 
persion  of  the  gain  and  of  the  refractive  index.  Since 
then,  experiments  have  been  carried  out  both  with  paral¬ 
lel  pumps  in  optical  fibers14  and  SOA’s16  17  and  with  or¬ 
thogonal  pumps.18-24  The  use  of  orthogonal  pumps  is 
particularly  attractive,  since  it  can  lead  to  polarization- 
independent  wavelength  conversion  and  optical  spectral 
inversion.  However,  only  very  recently  have  detailed  ex¬ 
perimental  data  been  published  demonstrating  the  rela¬ 


tive  independence  of  conversion  efficiency  on  the  pump 
wavelength  separation.15  17,25 

In  this  paper,  for  the  first  time  to  our  knowledge,  is  pre¬ 
sented  a  new  experimental  configuration  that  uses  a  DFB 
laser  with  two  lasing  modes  that  serve  as  dual  pumps. 
In  contrast  to  the  previous  studies,  which  require  two  ex¬ 
ternal  pump  waves  and  an  external  probe  signal  to  be  in¬ 
jected  into  a  SOA  or  optical  fiber,  the  lasing  modes  of  the 
DFB  laser  serve  as  the  two  pump  waves,  and  only  an  ex¬ 
ternal  probe  signal  is  needed.  In  Section  2  we  describe  a 
comprehensive  theory  for  the  two-pump  FWMs  which 
yields  analytical  formulas  for  the  small  detuning  case. 
In  Section  3  we  present  detailed  data  showing  the  spec¬ 
tral  dependence  of  the  conversion  efficiency  for  each  of  the 
important  conjugate  waves  and  discuss  some  of  the  inter¬ 
esting  physical  effects.  In  Section  4  we  offer  a  summary 
of  our  measurements  and  draw  some  brief  conclusions. 


2.  THEORETICAL  MODEL 

In  the  physical  model  of  the  laser  operation  and  the 
FWM,  several  simplifying  assumptions  are  made. 

First,  a  perturbative  approach  is  taken  in  that  it  is  as¬ 
sumed  that  the  power  of  the  probe  and  of  all  resulting 
conjugate  waves  is  negligible  compared  with  the  power  of 
each  pump  wave  and  that  these  negligible  values  can  be 
ignored  for  the  modeling  of  the  pump  waves. 

Second,  a  mean-field  approximation  is  employed  in 
which  the  intensify  of  all  the  fields  is  assumed  to  have 
little  spatial  variation.6,7  This  is  probably  unrealistic, 
since  the  presence  of  two  lasing  modes  almost  certainly 
means  that  there  is  some  spatial  dependence  involved. 
To  include  this  spatial  variation  in  the  analysis  of  the  la¬ 
ser  would  require  a  numerical  model  such  as  the  transfer- 
matrix  model,26  which  would  become  especially  compli¬ 
cated  with  the  correct  treatment  of  two  lasing  modes. 
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This  simplification,  while  physically  inaccurate,  yields 
analytical  solutions  and  gives  simple  explanations  for 
most  of  the  interesting  phenomena  observed  in  the  data. 

Third,  dispersion  in  both  the  refractive  index  and  the 
gain  is  ignored. 

Fourth,  the  nonlinear  gain  suppression  coefficient  is  in¬ 
cluded  but  is  assumed  to  be  constant.  This  allows  for  the 
inclusion  of  intersubband  processes  such  as  carrier  heat¬ 
ing  and  spectral  hole  burning,  which  are  important  for 
FWM  at  larger  pump-probe  detunings.  While  the  fre¬ 
quency  dependence  of  these  processes  is  ignored  in  this 
treatment,  the  detunings  studied  in  our  experiments  are 
small  enough  that  this  is  a  good  approximation. 

Figure  1  is  a  schematic  diagram  of  the  electric  field 
lines  relevant  to  this  study  of  two-pump  FWM.  The  total 
field  can  be  written  as 

E(tyz)  =  {Ea  +  EBexp{iM)  +  JEACexp(-iA0 
+  EBCex$(i2&t)  +  £1exp(-iflt) 

+  E2exp(iilt)  +  2?51exp[t(A  -  ft)*] 

+  2£B2exp[i(A  +  ft)*]  +  2£*exp[-i(A  -  ft)* ] 
+  EBxexp[i(2A  -  ft)*]}exp(-iwA*).  (1) 

In  this  expression  EA  and  EB  are  each  lasing  modes  of  the 
DFB  laser  that  serve  as  pump  waves.  a>A  represents  the 
lasing  angular  frequency  of  mode  A,  and  the  frequency 
separation  between  the  two  pumps  is  given  by  A.  The  ex¬ 
ternal  probe  wave  is  represented  by  Elf  and  its  frequency 
separation  from  pump  wave  A  is  given  by  ft.  Note  that, 
throughout  this  paper,  pump  A  is  defined  to  be  the  closest 
one  to  the  probe  wave,  so  which  mode  is  represented  by  A 
and  B  depends  on  the  location  of  the  probe  wave.  The 
pump  waves’  mixing  with  each  other  leads  to  the  creation 
of  cory  ugates  EAC  and  E BC .  The  interaction  of  the  probe 
with  the  two  pump  waves  leads  to  the  creation  of  the 
other  conjugates.  E2,  EB1,  and  EB2  are  created  mainly 
as  a  result  of  third-order  processes.  Terms  due  to  higher 
order  processes  are,  in  general,  left  out  inasmuch  as  they 
are  small.  The  exception  to  this  approximation  are  the 
fifth-order  processes  Ex  and  EBx ,  which  have  a  resonant 
condition  when  A  =  ft.  They  are  included  here  to  ex¬ 
plain  the  experimentally  observed  conjugate  waves  in  this 
regime. 


A.  Model  for  the  Ttao  Pump  Modes 

To  model  the  lasing  modes,  we  implement  the  following 

rate  equations28: 

d N  I  N 

—  =  — - vga(N  -  JV0)[1  -  e(SA  +  SB )] 

d*  qV  t 

X  ( SA  +  SB)t  (2) 

d SA 

—  =  Tvga(N  -  N0)[  1  -  e(SA  +  SB)]SA 
at 

SA 

- +  0R, p,  (3) 

Ta 

d SB 

—  =  Tvga(N  -  N0)[l  -  e(SA  +  SB)]SB 
at 

SB 

- +  PRb?>  (4) 

TB 

where  I  is  the  injection  current;  q  is  the  charge  of  an  elec¬ 
tron;  V  is  the  active-region  volume;  N  is  the  carrier  den¬ 
sity;  N0  is  the  transparency  carrier  density;  r  is  the  car¬ 
rier  lifetime;  vg  is  the  group  velocity;  a  is  the  differential 
gain;  T  is  the  optical  confinement  factor;  e  is  the  nonlinear 
gain  suppression  coefficient;  ^(-1  X  10~5)  is  the  ratio  of 
spontaneous  emission  photons  that  couple  into  the  lasing 
modes;  i?ap  is  the  spontaneous  emission  rate;  SA  and  SB 
are  the  density  of  photons  in  the  lasing  modes;  and  rA  and 
tb  are  the  photon  lifetimes  for  each  of  the  lasing  modes  A 
and  Bt  respectively.  In  general,  since  the  two  lasing 
modes  are  at  separate  wavelengths,  there  would  be  differ¬ 
ent  values  of  a  and  N0  for  each  mode.  In  DFB  lasers  the 
modal  behavior  is  dominated  by  the  relative  values  of  the 
cavity  loss  (i.e.,  photon  lifetime)  for  each  mode.  While 
the  dispersion  of  the  gain  will  have  a  minor  effect  on  the 
strength  of  the  modes,  these  effects  are  overwhelmed  by 
the  relative  values  of  rA  and  rB .  Therefore  we  make  the 
simplifying  assumption  that  a  and  Nq  are  constant. 

In  the  steady  state  the  left-hand  side  of  these  equations 
can  be  set  to  zero.  It  is  relatively  straightforward  to 
solve  for  the  photon  densities  of  each  mode  as  a  function 
of  the  carrier  density: 

y  ~  rBx 

SA  =  ta - ,  (5) 

TA  ”  tb 


Pump  B  Pump  A 

Eb  EX 


o  rAx-  y 
SB  -  tb  > 
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Rewriting  Eq.  (3)  in  terms  of  the  carrier  density  gives 


N.= 


(SA/rpl)  -  pR 


sp 


Tvga[l  -  e(SA  +  SB)]SA 


+  N0. 


(12) 


Using  Eqs.  (5>— (11)  in  Eq.  (12)  yields  a  single  nonlinear 
equation  for  the  carrier  density.  Given  a  value  for  the  in¬ 
jection  current,  one  may  solve  this  equation  numerically. 
Note  that  the  sign  in  Eq.  (7)  should  be  chosen  to  yield  a 
physically  allowable  solution;  that  is,  SA  and  SB  are  posi¬ 
tive.  pRe p  is  included  so  that  the  light  power  output  be¬ 
low  and  above  threshold  operation  can  be  modeled  with 
the  same  rate  equation. 


The  components  of  the  polarization  can  be  found  from 

P  =  €0XE,  (22) 

where29 

nc 

X  = - vga(N  -  N0)[a  +  i{l  -  e|E|2)].  (23) 

(Oq 

These  lead  to  the  following  components  of  the  polariza¬ 
tion: 


Pac  ”  ~ 


c0nca 

"A 


({SNA[a  +  i(l  -  e\E,\2)] 


B.  No  Probe  Case 

When  there  is  no  probe  wave  present,  FWM  still  exists 
owing  to  the  presence  of  both  lasing  modes  in  the  DFB 
cavity.  To  model  the  carrier  density  beating  that  is  due 
to  the  lasing  modes,  it  is  necessary  to  rewrite  Eq.  (2)  in 
terms  of  the  electric  field  (the  fields  are  assumed  to  be 
normalized  so  that  their  intensities  are  equal  to  the  pho¬ 
ton  density). 


6N 

dT 


I 

qV 


-  -  Vgd(N  -  N0)[  1  -  e|E(«)|2]|£(Ol2. 

(13) 


Assuming  that  the  carrier  density  has  the  form29 


N(t)  =  N,  +  «52V4exp<-iAi)  +  SNA*exp(ibt)  (14) 

and  using  the  rotating-wave  approximation,  we  can 
straightforwardly  solve  for  both  the  steady-state  carrier 
density  and  the  sinusoidally  varying  component29 


—  +  VgarN0(l  -  e|£,|2)|£U2 

1  +  B<ar(l  -  e|f?,|2)|£.|2  ’ 

Vgar{N,  -  N0)d  ~  2e\E.\2)EAEB* 

1  +  vgaT{\  -  e\Et\2)\E,\2  -  iAr  ’ 


(15) 

(16) 


where 

\E.\2  =  \Ea\2  +  \Eb\2  =  SA  +  SB. 

To  treat  the  electric  fields,  the  wave  equation 

n 2  c?E  1 
_  = 


(17) 


(18) 


V2E-  2  „  -  2  „ 

c2  dt 2  e0c2  St2 

is  used  with  the  electric  fields  and  the  polarizations  bro- 
ken  into  Fourier  components  as 

E(x,y,z,t)  =  t/(*,y)2  Ej(z)exp(-ia)jt),  (19) 


P(x,y,z,t)  =  U(x,y)2t  Pj(z)exp(-ia>jt).  (20) 
j 


Using  Eqs.  (19)  and  (20)  in  Eq.  (18),  multiplying  by 
U*(xty)7  and  integrating  over  the  x  and  y  dimensions 
yields  the  one-dimensional  wave  equation29 


d% 

d?" 


2,.  2 


it  to 


Ttoj 


C*  J  €qC2  j 


(21) 


-  ie{Ns-N0)EAEB*}EA 
+  (N,  -  N0)[a  +  i(l  -  €\E,\2)]Eac),  (24) 


60  nca 

PBC  = - ({SNA*[a  +  i(l  -  e|£,|2)] 

*>A 


-  i e(N,  -  N0)Ea*Eb}Eb 

+  ( N ,  -  N0)[a  +  id~  e\E,\2)]EBC).  (25) 

In  the  case  of  no  probe  wave  the  components  of  the  elec¬ 
tric  field  are  given  by 


Ej  -  Aj  exp(iAyz),  (26) 

where  kj  =  n  a>j  /c .  By  inserting  these  fields  into  the  one¬ 
dimensional  wave  equation  (21)  and  applying  the  slowly 
varying  envelope  approximation  [{d2Aj)!dz2  =  0],  we  ob¬ 
tain  the  following  coupled-mode  equations: 


dAj 

dz 


io>jT 
2  U€qC 


Pj 


exp  {-ikjz). 


(27) 


In  general,  the  coupled-mode  equations  are  solved  for  the 
spatial  variation  of  each  of  the  pump,  probe,  and  conju¬ 
gate  waves.  Here  we  make  the  simplifying  assumption 
that  the  waves  behave  as  if  they  were  in  an  optical  am¬ 
plifier.  Under  these  conditions  the  spatial  dependence 
for  each  pump  wave  is  solved  for.  These  solutions  are 
then  inserted  into  the  coupled-mode  rate  equations  for 
the  conjugate.  This  leads  to  the  following  expressions  for 
the  conjugate  waves: 


1 

Eac  =  2(1  -  e|2?,|2) 

hjard  ~  2c|g,|2)[g  +  i(l  ~  g|£«l2)3  . 

X{  1  +  vgr(l.-  e|E,|2)|Et|2  -  iAr 

X  Ea2Eb*,  (28) 

1 

Ebc  =  2(1  -  e|£,|2) 

«yxT(l  -  2e\E.\2)[a  +  i(  1  -  6|£,|2)]  . 

X  j  1  +  vgrd  -  e|JE,|2)|B,|2  +  iAr  +  “ 

x  Eb2Ea*.  (29) 


C.  External  Probe  with  a  Detuning  ft  <  A 

In  the  situation  in  which  an  injected  probe  wave  with  a 

pump-probe  detuning  ft  is  much  smaller  than  the  pump 
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separation  A,  the  copjugates  £2,  &bi*  Eb2  become 
significant.  The  carrier  density  in  this  case  is  expressed 


for  the  pump  waves.  Similarly,  for  the  probe  and  the 
conjugate  waves,  but  using  Eq.  (39)  to  eliminate  k A  and 
kB ,  we  obtain 


N(t)  =  Ns  +  <Wftexp(-ita)  +  SNn*  expiitlt).  (30) 

The  carrier  beating  due  to  the  two  pump  waves  has  been 
neglected,  since  it  has  no  interactions  with  the  slower 
exp(-ifU)  terms.  If  we  insert  Eq.  (30)  into  Eq.  (13)  and 
use  the  rotating-wave  approximation,  the  carrier  beating 

.  •  29 

term  is  given  as 

(N,  -  AT0)(1  -  2e\E,\2)Sa/P, 

SNa~  1  +  (1  -  c\E.\2)Po  ~  Mr  ’  (31) 

where  P,  =  1 /vgar,  Pq  =  \E,\2IP, ,  and 

Sa  =  ExEa*  +  EAE2*  +  EB1EB*  +  EBEB2 *.  (32) 
The  components  of  the  polarization  can  be  written  as 

PA  =  — - - (Af#  -  A^o)ta  +  i(l  "  4®*I2)1®a>  (33) 

WA 


<oaT 
2  n2 


—  (N,-N0)[a  +  i(l~€\Ea\2)] 
<*>A 


-  131 1  jp 

X  E ini - + 


>131 


e0 


(40) 


-he 


2B2 


waT  nca 

JL,  _ {Nt-N0) 
2  n2  o)A 


X  [a  +  t(l  -  €\Es\2)]E2#2  “ 


P  232 
€0 


(41) 


where  is  the  incident  probe  wave,  *i  is  a  coupling  co¬ 
efficient,  and  kb  1  is  0  in  the  equation  for  EB 1  • 

Using  the  rotating-wave  approximation,  we  obtain 
from  Eqs.  (39)-(41)  the  following  matrix  equation: 


pB  m  -  N0)[a  +  i(l  -  e\E.\2)]EB,  (34) 

*>A 


Pt  =  -  N0)[a  +  i(l  -  e\E,\2)]Ei 

WA 

+  {SNn[a  +  i(  1  -  e|E,|2)] 

-  ie(N,  -  N0)Sn}EA),  (35) 


P  =  .((N,  -  No) O  +  i(l  -  e\E,\2)]E2 

o>A 

+  {SNa*[<*  +  *d  -  *|£.l2)] 

-  ie(N,  -  N0)Sn*}EA),  (36) 

PB1  =  _!2^((J V,  -  N0)[a  +  i(l  -  e|£.|2)]£Bi 
"A 

+  {SNn[a  +  i(l  -  e|E,|2)] 

-  ie(N,  -  N0)Sa}EB),  (37) 


pB2  =  -  N0)[a  +  i(l  -  e\E,\2)]EB2 

"A 

+  {8Na*[a  +  i(l  -  e\E,\2)) 


-  if(N,  -  N0)Sn*}EB).  (38) 


Since  fl  is  small,  the  probes  and  the  conjugates  are  as¬ 
sumed  to  have  the  same  spatial  dependence  as  the  pump 
wave.30  Thus  the  fields  can  be  expressed  as  EA ,  E  x ,  and 
E2j  equivalent  to  Ajexp(ikAz)t  and  EBt  EB i,  and  EB 2, 
equivalent  to  Aj  exp (ik&z).  Using  these  fields  in  the  one¬ 
dimensional  wave  equation  (21)  and  assuming  that  the 
first  and  the  second  derivatives  of  the  field  amplitudes 
(Aj)  vanish  with  respect  to  z  in  the  leading-order  approxi¬ 
mation,  we  obtain 


-kA,B2EA,B  + 


n2o>A,B2  ^  <»AtB2T  n 

—a~Ea*  ~  ~i^~Pa-b 


(39) 


Ei  1 

■  kEC 

E2 * 

0 

Eb  i 

0 

.  eB2*. 

.  0  . 

M 


where  the  matrix  elements  are  given  as 


i^EA*\ 


Mn  =  ft  -  Jf|  “  ie  I  I^aI2. 

/  B 

M12=  -K^-—-ie]EA2, 

M*  =  ie\EB*EA, 

Mu  =  ~~p^  ~  *e  I EBEA< 

I  B*  • 
M2'=-K(-jrX  +  l€ 

I  B* 

Mn  =  -H  -  K(  +  ie  |  \EA\2, 

I  B* 

M23=  -K\[-—+ie\EA*EB*, 
M24  =  ~K[~^  +  ie\  ea*Eb, 

M3i  =  -Jf(  “  i*\EA*EB  > 


=  -e\  ~  *cl eaeb, 


B 

Mk^ci-kI  -i€||£Bl2, 


PA 


(42) 
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M, u  = 


M4 1  = 


M* 2  = 


-AT  “  ie \Eb2, 


-7TT  +  “  MV  • 


M4,  = 


M«=  +  icj|£B|2, 


where 


A  =  1  +  (1  -  e|£\|2)P0  "  iO*i 
B  -  [a  +  i(l  -  «|P,|2)](1  -  2c|B,|2), 
Tca(N.  -  N0) 


Solving  this  system,  we  obtain 


—  =  (M-'h!  = 

xEi 


ea*ebu* 


=  (ifr-1)*  = 


ea*eb*u 


where 


Tea  (AT,  -  N0) 


{-[a  +  i(l  -  e|B.|2)] 


2nP,  11  ' 

X  (1  -  2e|£,|2)  +  i«P,(l  +  Po  ~  iftr)},  (48) 
D  =  n(-n(l  +  P0)  -  »(flji2  -  ft2)r  +  e{ftP,Po2 
-  in*MP,(l  +  Po  -  Mr)  ~  2P^0]}).  W9) 


ft  2  - - - - — .  (50) 

R  rTVga(N,  -  N0) 

Note  that  Eq.  (45)  reduces  to  the  equivalent  of  Eq.  (67) 
of  Ref.  6  in  the  limit  of  no  second  pump  (|EB|2  -  0)  and 
no  nonlinear  gain  (e  =  0).  Also,  this  formulation  is 
equivalent  to  a  formulation  in  which  the  van  der  Pol 
equation  is  used  as  in  Ref.  7  for  the  single-pump  case.  In 
the  case  of  two  pumps  a  second  van  der  Pol  equation  de¬ 
scribing  the  evolution  of  the  waves  near  the  second  pump 
would  be  used.  Insofar  as  the  definition  of  Et  in  Ref.  7 
differs  from  ours  by  a  factor  of  £,  in  the  limit  of  no  second 
pump  (|Eb|2  =  0)  our  matrix  M  [in  Eqs.  (42)— (44)]  re¬ 
duces  to  a  2  x  2  matrix  that  can  be  shown  to  be  identical 
to  Eqs.  (11M14)  of  Ref.  7. 

Finally,  it  is  interesting  that  Eqs.  (42)-(44)  contain  no 
dependence  on  the  pump  separation  (A).  With  dispersion 
effects  neglected,  this  indicates  that  the  conversion  effi¬ 


ciency  for  the  shifted  conjugates  EB\  and  E B2  should  be 
the  same  as  that  for  E2  no  matter  what  the  pump  sepa¬ 
ration  is.  In  practice,  the  conversion  range  is  limited  by 
the  window  of  appreciable  net  gain  for  the  semiconductor 
material. 

D.  Resonance  Case  when  the  Probe— Pump  Detuning  ft 
~  A 

If  the  injected  probe  wave  has  a  larger  pump-probe  de¬ 
tuning  such  that  it  is  of  the  order  of  the  pump  spacing,  we 
can  no  longer  neglect  the  carrier  density  beating  that  is 
due  to  the  two  pumps  [the  exp(-£A£)  term].  In  this  case 
we  express  the  carrier  density  as 

N(t)  =  N8  +  SNAexp(-iM) 

+  SNa*  exp(iAJ)  +  SN^ 

+  SN(i*  exp(£fl£)  +  5Nxexp[-i(A  -  fl)£] 

+  SNX*  exp[i(A  -  Cl)t],  (51) 

where  SNX  represents  the  carrier  density  beating  that  is 
due  to  the  mixing  of  the  SN n  and  the  SN A  terms.  Insert¬ 
ing  Eq.  (51)  into  Eq.  (13)  gives  the  following  system  of 
equations: 


SNn 

fQil 

[SNX* 

■  M 

where 

Pn  =  iSlr-  1  -  (1  -  e\E,\2)P0, 

P12  =  ~(1  -  2e\Es\2)EAEB\ 

P 21  =  ~(1  -  2e\E,\2)EA*EB , 

*  $ 

P-a  =  -  ft)r  -  1  -  (1  -  e\E,\2)P0,  (53) 

Qi  =  -  N0)(l  -  2e\E,\2)EA*Eu 

*1 

Q2  =  ^-[6NS(1  -  2e\E,\2) 

-  2e(N,  -  N0)EA*EB]EA*E1 .  (54) 

These  equations  yield  solutions  of 

sNa  i  i  r  Pa  -PuirQil  .... 

[SNX*\  PUP&- PuPzA-Pn  pii  JIQ2J 

When  A  -  Cl  becomes  small  there  is  a  large  resonance  in 
the  SNX  terms.  This  beating  causes  each  pump  to  scatter 
into  the  conjugates  detuned  at  this  frequency.  As  a  re¬ 
sult,  E2  and  Ebi  experience  a  large  enhancement,  and 
two  new  conjugates,  Ex  and  EBxf  become  observable. 
Relationships  for  the  electric  fields  can  be  written,  result¬ 
ing  in  a  large  set  of  linear  equations.  We  leave  this 
analysis  out  because  it  becomes  very  difficult  to  track  and 
because  it  yields  little  insight  into  the  physics  of  the  sys¬ 
tem.  However,  we  do  summarize  here  the  field  depen¬ 
dencies  of  the  resulting  conjugates  from  this  fifth-order 
process: 
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Ex  =  Xx5)Eb\Ea\2Ea*Eu 
Ebi  =  X^Eb'E^ES, 

E2  =  xiS)EB2EA*2Elt 
EBx  =  x^IEbWES.  (56) 

3.  EXPERIMENTS 

For  our  experiments,  we  use  an  index-coupled  DFB  laser 
with  two  lasing  modes  as  the  mixing  device.  Another 
DFB  laser  similar  to  the  mixing  laser,  but  with  a  single 
mode,  is  used  to  provide  the  probe  signal.  The  probe  sig¬ 
nal  is  coupled  into  a  fiber  and  is  passed  through  an  optical 
isolator  to  eliminate  injection  of  light  from  the  mixing 
DFB  laser  into  the  probe  DFB  laser.  The  probe  light  is 
then  passed  through  a  20-dB  fiber  coupler  to  permit  moni¬ 
toring  of  the  power  of  the  probe  signal  injected  into  the 
mixing  DFB  laser.  Light  emitted  from  the  other  facet  of 
the  mixing  DFB  laser  is  coupled  into  another  fiber  and  is 
monitored  with  a  Hewlett-Packard  optical  spectrum  ana¬ 
lyzer  to  yield  the  output  powers  of  the  pump,  probe,  and 
conjugate  signals.  The  operating  temperatures  of  both 
probe  and  mixing  DFB  lasers  are  maintained  by  means  of 
thermoelectric  controllers,  with  the  mixing  DFB  laser  be¬ 
ing  kept  at  25  °C  for  all  the  experiments.  Tuning  of  the 
probe  signal  is  accomplished  by  changing  of  the  operating 
bias  current  and  of  the  temperature  of  the  probe  DFB  la¬ 
ser.  A  larger  tuning  range  is  achieved  by  use  of  several 
different  probe  DFB  lasers. 

4.  RESULTS  AND  DISCUSSION 

Figure  2  shows  a  typical  spectrum  from  the  antireflection 
facet  of  the  two-mode  DFB  laser  at  a  bias  current  of  I 
=  90  mA.  Again,  EA  and  EB  are  the  two  lasing  modes, 
and  EAC  and  EBC  are  conjugates  caused  by  the  mixing  of 
the  lasing  modes.  Curves  showing  light  output  power 
versus  bias  current  for  both  the  lasing  modes  and  the  con¬ 
jugates  are  shown  by  means  of  symbols  in  Fig.  3(a).  The 
power  shown  corresponds  to  the  peak  power  of  each  mode 
measured.  The  dashed  curves  represent  a  fit  of  the  pho¬ 
ton  density  in  each  mode,  obtained  from  Eqs.  (5)— (12), 
(28),  and  (29),  with  the  parameters  listed  in  Table  1.  The 
relationship  between  the  photon  densities  S  of  each  mode 
and  the  measured  peak  power  PmeafJ  can  be  expressed  as 

Praess  =  KxhvVgAS,  (57) 

where  Kx  is  a  coupling  constant,  h  v  is  the  photon  energy 
(—0.8  eV),  vg  is  the  group  velocity,  A  is  the  cross-sectional 
area  of  the  waveguide  (— 1  ^m  x  0.25  /im),  and  S  is  the 
photon  density.  Our  fit  yields  a  value  of  —0.44  for  Ki, 
which  is  quite  reasonable.  Figure  3(b)  shows  the  mea¬ 
sured  conjugate  powers  versus  the  following  power  laws, 
which  they  should  obey: 

Sac  a  Sa2Sb  , 

SBc  a  SASB2.  (58) 

The  slope  of  1  on  the  log-log  plot  indicates  that  the  con¬ 
jugates  do  indeed  obey  these  relationships. 


Figure  4  shows  typical  spectra  of  the  FWM  that  results 
when  a  probe  wave  is  injected  at  a  pump-probe  detuning 
that  is  much  smaller  than  the  pump  spacing.  These 
spectra  are  for  three  different  probe  detunings.  Also,  it 
should  be  noted  that  these  spectra  are  taken  from  the 
highly  reflecting  facet  of  the  mixing  DFB  laser.  The  mix¬ 
ing  of  the  probe  wave  with  the  first  pump  creates  three 
main  conjugate  waves.  Here  one  of  the  advantages  of  us¬ 
ing  two  pumps  for  FWM  is  very  clear.  Each  of  the  con¬ 
jugate  waves  is  approximately  the  same  power,  even 
though  EB2  is  shifted  by  an  extra  2  nm  (the  separation  of 
the  two  pump  waves)  relative  to  E2 ,  which  is  the  result¬ 
ing  conjugate  if  just  EA  were  present. 

This  result  can  be  seen  more  explicitly  in  Fig.  5,  where 
the  conversion  efficiency,  defined  as  the  ratio  of  the  out¬ 
put  conjugate  power  to  the  input  probe  power,  is  plotted 

10° 

10*1 
C  10'2 

I  Iff3 

g  10'4 
£  10'5 
10*6 
10*7 

1528  1530  1532  1534  1536 

Wavelength  (nm) 

Fig.  2.  Typical  lasing  power  spectrum  of  the  mixing  two-mode 
DFB  laser.  EA  and  EB  are  the  two  lasing  modes,  and  EAC  and 
Ebc  are  the  conjugates  created  by  the  mixing  of  these  modes. 


Fig.  3.  (a)  Light  output  power  versus  the  bias  current  (symbols), 
for  the  two  lasing  modes  as  well  as  the  created  conjugates  shown 
in  Fig.  2.  The  dashed  curves  represent  theoretical  fits  of  the 
photon  density  in  each  mode,  (b)  Power  of  each  conjugate  wave 
plotted  versus  the  proper  power  dependence  on  the  lasing  field 
powers  as  described  in  the  text  [Eqs.  (58)]. 
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Table  1.  Parameters  Used  in  the  Theoretical 
Calculations 


Parameter 

Symbol 

(Unit) 

Value 

Electron  charge 

<7(0 

1.602  X  10“19 

Active-region  volume 

V  (cm3) 

1  x  HTU 

Carrier  lifetime 

t(b) 

0.8  X  10“9 

Group  velocity 

Vg  (cm  s"1) 

8.33  x  109 

Differential  gain 

a  (cm2) 

2  x  10‘16 

Transparency  carrier  density 

N0  (cm"3) 

1  x  1018 

Nonlinear  gain  coefficient 

e  (cm3) 

7  x  10“ 17 

Optical  confinement  factor 

r 

0.1 

Photon  lifetime 

Tpl  (8) 

1.33333  x  10““ 

Tp  2  (S) 

1.33332  x  10““ 

Linewidth  enhancement 
factor 

a 

5 

Spontaneous  emission 

pRtv  (cm“3  s"1) 

2800 

Wavelength  (nm) 


Pig.  4.  Typical  FWM  power  spectra  for  the  case  of  an  injected 
probe  wave  with  a  small  detuning  on  the  left-hand  side  of  pump 
Ea  .  For  each  probe  detuning  shown,  three  conjugate  waves  E2 , 
EB\ »  and  EB2  are  produced. 


as  symbols  for  each  of  the  corrugates.  The  wavelength  in 
this  plot  corresponds  to  the  wavelength  of  the  probe  wave. 
Also,  it  should  be  noted  that  the  definitions  of  each  conju¬ 
gate  depend  on  the  location  of  the  probe  wave;  the  pump 
closest  to  the  probe  wave  is  always  considered  to be  EA . 
We  can  see  that  there  is  only  a  small  variation  among 
these  three  conjugates,  especially  at  small  pump-probe 
detunings.  These  variations  can  be  accounted  for  by  dif¬ 
ferences  in  the  gain  and  cavity  enhancement  experienced 
by  each  of  the  conjugates. 

The  curves  in  Fig.  5  correspond  to  the  theoretical  re¬ 
sults  that  we  obtained  by  using  the  absolute  magnitude 
squared  of  Eqs.  (45)  and  (47).  The  solid  curves  represent 
the  conversion  efficiency  of  E2 ,  and  the  dashed  curve  rep¬ 
resents  the  conversion  efficiency  of  the  shifted  conjugates. 
Differences  in  the  two  theoretical  curves  are  due  to  the 
difference  in  power  of  the  two  pump  waves.  The  param¬ 
eters  for  this  fit  are  again  those  listed  in  Table  1.  Note 
that  the  theoretical  conversion  efficiency  spectra  do  not 
take  into  account  either  dispersion  of  the  gain  and  of  the 
refractive  index  or  the  spatial  dependence  of  the  two 
pump  waves.  These  effects  could  account  for  the  fact 
that  the  theory  predicts  a  larger  difference  than  is  mea¬ 


sured  between  the  conjugates  in  the  cases  in  which  the 
probe  wave  is  near  one  of  the  pumps.  The  theory  also  ig¬ 
nores  the  presence  of  the  DFB  grating  as  well  as  the  end 
facets,  leaving  out  the  detailed  effects  of  the  cavity  on  the 
conversion  efficiency  spectrum.  This  is  why  the  theory 
does  not  predict  much  of  the  behavior  for  the  case  in 
which  the  probe  is  away  from  the  pump  wavelengths. 

To  demonstrate  more  clearly  the  relative  independence 
of  the  conversion  efficiency  on  pump  separation,  we  re¬ 
peated  the  experiments  for  a  similar  DFB  laser  with  two 
lasing  modes  ~6  nm  apart.  Figure  6  shows  the  conver¬ 
sion  efficiency  for  the  three  conjugates  (denoted  by  sym¬ 
bols)  as  well  as  the  amplified  spontaneous  emission  or  the 
lasing  spectrum  (denoted  by  the  solid  curve).  Again, 
there  is  a  relatively  small  difference  between  the  conver¬ 
sion  efficiency  of  E2  and  the  shifted  conjugate  EB2f  even 
though  the  extra  shift  is  6  nm  (—750  GHz).  In  Fig.  7  we 
show  the  conversion  efficiency  of  EB2  relative  to  that  of 
the  conjugate  in  an  experiment  conducted  with  a  similar 
DFB  laser,  but  with  a  single  lasing  mode  (single  pump 
wave).9  The  x  axis  represents  the  pump-probe  detuning 
frequency  for  the  single-pump  case,  and,  for  the  two- 
pump  result,  zero  detuning  represents  halfway  between 
the  two  modes.  This  was  done  so  that  the  two  cases 


Probe  Wavelength  (nm) 


Fig.  5.  Conversion  efficiency  (symbols)  for  each  of  the  generated 
conjugate  waves,  plotted  versus  the  probe  wavelength.  The 
curves  are  theoretical  fits  of  the  conversion  efficiency,  with  the 
solid  curves  representing  E2  and  the  dashed  curves  representing 
the  shifted  conjugates  (EBi  and  EB2).  There  is  relatively  little 
difference  in  conversion  efficiency  for  the  first  conjugate  and  the 
shifted  conjugates. 
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Fig.  6.  Conversion  efficiency  (symbols)  for  the  conjugate  waves 
in  another  DFB  laser  with  two  pump  wavelengths  —6  nm  apart. 
The  solid  curve  represents  the  laser  power  spectrum  or  the  am¬ 
plified  spontaneous  emission  spectrum  without  the  injected 
probe  wave.  Again,  there  is  relatively  little  difference  in  conver¬ 
sion  efficiency  for  the  first  conjugate  and  the  shifted  conjugates. 


60  J.  Opt.  Soc.  Am.  B/Vol.  17,  No.  1/January  2000 


J.  Minch  and  S.  L.  Chuang 


Fig.  7.  Comparison  of  the  conversion  efficiency  for  two-pump 
FWM  and  single-pump  FWM9  in  DFB  lasers.  The  conversion  ef¬ 
ficiency  is  shown  versus  the  probe  detuning  for  the  single-pump 
case.  For  the  two-pump  case,  zero  detuning  is  chosen  to  be  mid¬ 
way  between  the  pump  waves  so  that  a  comparison  can  be  made 
on  the  basis  of  the  same  overall  frequency  shift.  For  small 
enough  probe  detunings,  the  two-pump  case  gives  conversion  ef¬ 
ficiencies  between  100  and  1000  times  larger. 
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Fig.  8.  (a)  Typical  FWM  power  spectra  for  the  case  in  which  the 
probe  detuning  is  close  to  the  probe  separation.  In  this  case  con¬ 
jugates  E2  and  EBi  experience  a  large  enhancement,  while  new 
conjugates  Ex  and  EBx  are  also  created,  (b)  Close-up  of  newly 
created  conjugate  EBx . 

could  be  compared  on  the  basis  of  the  same  frequency 
shift.  The  two-pump  case  can  give  a  100- 1000 X  increase 
in  conversion  efficiency  if  the  probe  is  close  enough  to  one 
of  the  pump  waves.  Note  that  there  could  be  a  difference 
in  coupling  coefficients  between  the  two  cases,  which 
would  change  the  relative  values  of  these  two  conversion 
efficiency  curves;  however,  we  estimate  that  this  should 
account  for  no  more  than  a  factor  of  5  between  the  two 
curves.  Also,  the  difference  in  pump  power  between  the 
two  cases  needs  to  be  considered  in  this  comparison.  The 
exact  values  of  the  pump  powers  are  difficult  to  obtain  be¬ 


cause  of  the  unknown  coupling  coefficients.  The  total, 
maximized  lasing  power  from  each  laser,  collected  from 
the  side  in  which  the  probe  wave  was  injected,  was  ~1 
mW  for  each  laser.  For  the  two-pump  case,  however,  this 
power  is  split  between  two  modes,  which  makes  the  con¬ 
version  efficiency  lower  than  if  the  same  pump  powers 
were  used. 

When  the  probe  detuning  becomes  comparable  with  the 
spacing  between  the  pump  waves,  some  interesting  phe¬ 
nomena  occur.  Figure  8(a)  shows  the  resulting  spectra 
for  several  probe  detunings  close  to  the  pump  separation 
frequency.  Two  of  the  conjugates  ( E2  and  EB{)  become 
strongly  enhanced  in  this  figure.  Also,  two  new  conju¬ 
gates  appear  (Ex  and  EBx)  and  become  as  strong  as  the 
others.  A  close-up  view  of  the  newly  created  EBx  is 
shown  in  Fig.  8(b). 

The  conversion  efficiencies  of  each  of  these  conjugates 
are  plotted  in  Fig.  9.  Again,  the  x  axis  represents  the  po¬ 
sition  of  the  probe  wave.  The  point  at  which  the  pump- 
probe  detuning  is  equivalent  to  the  pump  separation  is 
— 1535.3  nm.  As  the  probe  wave  approaches  this  location 
each  of  the  conjugates  experiences  a  large  enhancement. 
When  the  probe  wave  is  away  from  this  location,  Ex  and 
EBx  disappear,  while  E2  and  EBX  are  still  observed.  As 
pointed  out  in  section  2,  this  happens  because  Ex  and  EBx 
are  created  by  a  fifth-order  process,  while  E2  and  EBi 
have  contributions  from  both  the  third-order  FWM  pro¬ 
cesses  and  this  resonant  fifth-order  process.  Where  the 
fifth-order  process  no  longer  dominates,  there  are  still 
contributions  to  these  conjugates. 
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'  Fig.  9.  Conversion  efficiencies  for  the  four  conjugates  that  expe¬ 
rience  a  large  resonance  when  the  probe  detuning  is  near  the 
pump  separation. 
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Fig.  10.  Normalized  magnitude  of  the  carrier  beating  terms  im¬ 
portant  in  the  creation  of  the  observed  conjugates.  The  x  axis 
corresponds  to  the  position  of  the  injected  probe  wave.  When 
the  probe  detuning  is  near  the  pump  separation  frequency  the 
large  resonance  in  SNX  causes  the  observed  enhancement  of  the 
conjugates  shown  in  Figs.  8  and  9. 
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Figure  10  shows  a  plot  of  the  relative  strengths  of  the 
carrier  beating  processes  SN n  and  SNXf  from  Eq.  (55), 
normalized  to  the  power  of  the  probe  wave.  For  probe 
wavelengths  that  give  small  pump-probe  detunings,  SN 
is  large  and  causes  the  large  enhancement  in  the  conver¬ 
sion  efficiency  seen  for  this  regime,  as  discussed  above. 

As  the  detuning  becomes  larger,  SN^  drops  off  rapidly. 
However,  as  the  detuning  becomes  closer  to  the  pump 
separation,  SNX  has  a  large  resonance,  which  leads  to  the 
large  enhancement  of  each  of  the  conjugates  shown  in 
Figs.  8  and  9. 

5.  SUMMARY  AND  CONCLUSIONS 

We  have  presented  what  is  for  the  first  time,  to  the  best  of 
our  knowledge,  a  comprehensive  study  of  the  properties  of 
four-wave  mixing  (FWM)  in  a  distributed  feedback  (DFB) 
laser  with  two  lasing  modes.  We  show  detailed  experi¬ 
mental  measurements  of  the  two-pump  FWM  conversion 
efficiency  spectrum  and  show  that,  to  first  order,  there  is 
no  fundamental  dependence  of  the  conversion  efficiency 
on  pump  separation.  Any  dependence  comes  from 
smaller  effects  such  as  cavity  effects  and  dispersion  of  the 
gain  and  of  the  refractive  index.  We  have  also  presented 
a  comprehensive  model  for  the  lasing  modes  and  the 
FWM  based  on  a  mean-field  approximation.  This  model 
agrees  quite  well  with  the  measured  data,  showing  no  de¬ 
pendence  on  pump  separation  for  the  two-pump  FWM 
conversion  efficiency.  Also,  it  explains  how  the  carrier  t 
density  beating  induced  by  the  pump-pump  beating  and 
the  pump-probe  beating  can  interact  and  lead  to  a  new 
resonant  enhancement  of  some  conjugate  waves,  as  ob¬ 
served  in  the  experimental  data  for  probe  detunings  of 
the  order  of  the  pump  separation. 

The  practicality  of  using  two-mode  DFB  lasers  for 
wavelength  conversion  is  severely  limited  by  the  difficulty 
in  building  such  a  device,  especially  with  a  tunable  sepa¬ 
ration  between  the  lasing  modes.  Little  work  has  been 
done  in  this  area.31  In  practice,  one  can  use  either  a 
single-mode  DFB  laser  with  one  external  pump,  which 
provides  a  tunable  detuning,  or  two  external  pump  lasers 
injected  into  a  SOA,  as  was  done  in  many  of  the  studies 
cited  in  our  reference  list.  The  analytical  results  ob¬ 
tained  from  our  theory  can  be  applied  to  these  cases,  since 
a  mean-field  approach  (which  leaves  out  the  details  of  the 
laser  cavity)  has  been  used. 
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